Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,608,703	0	A→G	22.1%	12.3 ‑14.2 ‑0.0	72	intergenic (+54/+48)	zntA/tusA	zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/15):  ref base (37/16):  total (40/32)
Fisher's exact test strand distribution p-value = 6.77e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

CATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAG  >  NC_000913/3608505‑3608893                                                                                                                                                                                                           |                                                                                                                                                                                                  cATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGagag                                                                                                                                                                                                  >  4:190251/1‑200 (MQ=255)      cTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑a                                                                                                                                                                                            >  4:426126/1‑200 (MQ=255)        aTTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑gg                                                                                                                                                                                         >  4:403731/1‑200 (MQ=255)         ttGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  4:684412/200‑1 (MQ=255)               tGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAggg                                                                                                                                                                                  <  4:329713/200‑1 (MQ=255)                      gAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAgg                                                                                                                                                                          >  4:523345/1‑200 (MQ=255)                      gAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAgg                                                                                                                                                                          <  4:214109/200‑1 (MQ=255)                                                 gtgggggaggccggggTGGTGGCTGGAAGTGTGGGAAGAACGGGGGGGGCCGGTGCTGGGAACAGCAAAGGCGTAAGAATTTTG‑GCGCGGGAT‑AAAAGGAAACCCGTTCGCAACA‑TTGGGCGCGTTCGGCCCCCC‑CCCCCCCCGGGGGAGGGGGTT‑A‑GGGGGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                                                                                               <  4:215718/187‑1 (MQ=255)                                                 gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                                                                                               >  4:527175/1‑200 (MQ=255)                                                           acggggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGGGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGTAAACCGATCGCAACA‑TTGGTCGCGATCGGTCCCCT‑CGGCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  2:672542/147‑1 (MQ=255)                                                           aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  2:458194/150‑1 (MQ=255)                                                               ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCTCCT‑CGCCCCTCTGGGGAGAGGGt                                                                                                                                                                                               >  1:194192/1‑141 (MQ=255)                                                               ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGGGAGGGt                                                                                                                                                                                               >  1:83986/1‑141 (MQ=255)                                                                    tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  1:386058/141‑1 (MQ=255)                                                                    tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  1:393385/141‑1 (MQ=255)                                                                     ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                                                                           >  4:328308/1‑200 (MQ=255)                                                                             gtgcgtgcaaggaCGGGGGCGACGGTGGCGGGGACAGCAAATGCGTTAAGATTGTT‑GGGACGGAG‑ATCAGGGAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCGCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  3:400684/120‑1 (MQ=25)                                                                             gggCTGGCAGATACGGGAGCGACGGAGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  3:169783/130‑1 (MQ=255)                                                                             gTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑ggg                                                                                                                                                                                        <  3:735740/132‑1 (MQ=255)                                                                                 gagttgtgtaaaagggacgggTGCGGGTGACAGCGAATGCGTTAAGATTGTG‑CCGCAGGAG‑ATAAGGAAAACGGACCGAAACT‑TTGAGCGGGATCGGCCCCCC‑CCCCCCTCGGGGGAGGGGGTT‑G‑GGGTGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                               <  4:160037/181‑1 (MQ=255)                                                                                  ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                                              >  4:432880/1‑200 (MQ=255)                                                                                  cgCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCCCTGGGGAGCGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGac                                                                                                                                                                >  2:754500/2‑149 (MQ=255)                                                                                      gATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑C‑GGGTGAGGGG‑AAAAGGCGGCATCGACGCCa                                                                                                                                                            >  2:485981/1‑150 (MQ=255)                                                                                       aTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCCGGAG‑ATCAGGCAAACCGATCGCAACA‑TTGCGCCCGCTCGGTCCCCT‑CGCTCCTCTGCGGAGCGGGTT‑C‑GTGTGCGAGG‑aca                                                                                                                                                                             >  3:562353/1‑130 (MQ=255)                                                                                                cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGCGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCAt                                                                                                                                                                    >  3:331380/1‑132 (MQ=255)                                                                                                     gTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGACGCCAATCAGc                                                                                                                                                      >  1:359194/1‑141 (MQ=255)                                                                                                     gTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGc                                                                                                                                                      >  1:459623/1‑141 (MQ=255)                                                                                                     gTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCaa                                                                                                                                             >  2:160031/1‑150 (MQ=255)                                                                                                     gTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGCGAGGGG‑AAAAGGCGGCATCGAAg                                                                                                                                                               >  3:194206/1‑132 (MQ=255)                                                                                                          ggtgcaggggAATGCGTTAGGATTGTT‑GCGCAGAAG‑ATAAGGCAAACCGTGCGCAACA‑TGGGGCGCGATCGGCCCCCT‑CGCCCCTCTGGGGGGAGGTTT‑A‑GGGTGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                                                                        <  2:165389/141‑1 (MQ=255)                                                                                                               cAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCACCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                                                                     >  3:247642/1‑132 (MQ=255)                                                                                                                ggCGAATGCG‑TAAGATTGGT‑GGGGAGGGGAATAGGGCAAACGGGCCGAAACA‑TGGGGGGGCACGGGCCCCCC‑CCCCCCTCGGGGGGGGGGGTT‑G‑GGGGGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                                                    <  3:238826/131‑1 (MQ=255)                                                                                                                 gcaaaatttttAAGTTTGTG‑GCGCAGGAA‑AAAAGGCAAACCATTCGCAACA‑TGGGGCGCG‑TCGGGCCCCTCCGCCCATCGGGGGGGGGGTTG‑A‑GGGTGGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                                               <  4:666477/191‑1 (MQ=255)                                                                                                                             aagttgttgGCGAAGGAG‑AAAGGGCAAACCGTGCGAAAAT‑TTGGGCGCGATGGGTCCCCC‑CCCCCCTTGGGGGAGGGGGTT‑A‑GGGTGAGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                                    <  4:651592/197‑1 (MQ=255)                                                                                                                              gtttgtg‑gGGCAGGAG‑AGAAGGAAAACCGACCGCAACT‑TTAGGCGCGACAGGTCCCCC‑CGCCCCTGGGGGGAGGGGGTG‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                                                                    <  2:263338/148‑1 (MQ=255)                                                                                                                              gagtgtg‑gCGAAGGAG‑AGAAGACAAGCGGAT‑GTAACACTGGGGCGCGGTCGGTCCCGC‑CTCCCCT‑GGGGGAGGGGGTTGA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGAt                                                                                                                    <  2:426084/147‑1 (MQ=255)                                                                                                                                  gtt‑gCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGAGTTCCATAAAGGTACaaa                                                              >  4:390142/1‑200 (MQ=255)                                                                                                                                    t‑gCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑TGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                                            <  4:639220/200‑1 (MQ=255)                                                                                                                                     tggcgaggagaagaaGGAAAACGGACCGAAACA‑TGGGGGGGAATC‑GTCCCCC‑CGCCCCCTTGGGGGGGGGGTT‑A‑GGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                               <  2:403702/137‑1 (MQ=255)                                                                                                                                           ggAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATCAGCCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa                                                      >  4:133860/1‑200 (MQ=255)                                                                                                                                           ggAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                                                                 >  1:92625/1‑141 (MQ=255)                                                                                                                                           ggAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                                                                        >  2:432898/1‑150 (MQ=255)                                                                                                                                           ggAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGGCCAt                                                                                                        >  2:952446/1‑150 (MQ=255)                                                                                                                                            gag‑aTAAGGCAAACCGATTGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGTGGCATCGAAGCTAATCAGCCGCTCTTAAACGCCTTTACGAATCAAATa                                                                                                                         >  3:700545/1‑132 (MQ=37)                                                                                                                                            gag‑aAAGGGAAACCCGATCGCAACA‑TTGATCCCGCCCGTCCCCCT‑CCCCCCTCTGGGAAGAGGGTT‑G‑GGGTGGGGAA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                                                     <  4:215127/200‑1 (MQ=255)                                                                                                                                             gtaaaaaGGAAAAACGATCGAAACA‑TTGAGCGCGATGGGACCCCC‑CGCCCCTGTGGGGGGAGGGTT‑G‑GGGTGAGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                                                     <  4:381489/195‑1 (MQ=255)                                                                                                                                             ag‑aTAAGGCAAACTGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCTCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCATTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAATACCCAGGAat                                                    >  4:312279/1‑200 (MQ=255)                                                                                                                                             ag‑aTAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                               >  1:499972/1‑141 (MQ=255)                                                                                                                                                    ggCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                                                                  >  3:110292/1‑132 (MQ=255)                                                                                                                                                       aaccagctCGAAACA‑TTGAGCGGGATGGGTCCCTC‑CGCCCCTCGGGGGAGAGGGTT‑G‑GGAGGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                               <  3:239675/125‑1 (MQ=255)                                                                                                                                                             gtCGAAACA‑TGGAGCCCGATGGGCCCCCG‑CGCCCCCAGGGGGAG‑GGGTT‑G‑GGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                                                                       <  2:688392/149‑1 (MQ=255)                                                                                                                                                                   aCA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc                               >  4:862166/1‑200 (MQ=255)                                                                                                                                                                      aatggggggcagcGGCCCTCT‑CG‑CCCTTTGGGGGGGGGGTT‑AGGGGTGGGGGA‑AAAAGGAGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCGAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                                               <  2:455647/138‑1 (MQ=255)                                                                                                                                                                        tgcgggcgcccgccccccC‑CCCCCCTCTGGGG‑GGGGGTT‑A‑GGGTGAGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                                                                            <  2:49068/136‑1 (MQ=255)                                                                                                                                                                         gAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                                                                     >  1:70348/1‑141 (MQ=255)                                                                                                                                                                         gAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                                                                     >  1:192926/1‑141 (MQ=255)                                                                                                                                                                            ggagaccGGTCCCC‑‑CGCCCCCGGTGGGGGGGGGTT‑AGGGGTGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                                                                           <  3:589642/126‑1 (MQ=255)                                                                                                                                                                                aTCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                                                              >  1:470542/1‑141 (MQ=255)                                                                                                                                                                                aTCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAAt                   >  4:8905/1‑200 (MQ=255)                                                                                                                                                                                    gTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                                                                          >  1:204716/1‑141 (MQ=255)                                                                                                                                                                                    gTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGc               >  4:697320/1‑200 (MQ=255)                                                                                                                                                                                     cccccT‑CCCCCCTGGGGGGGGGGGTTT‑A‑GGGTGGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                                                                         <  1:110291/140‑1 (MQ=255)                                                                                                                                                                                      ccccT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                                    >  2:418320/1‑145 (MQ=255)                                                                                                                                                                                           tccccccTCAGGGAAGGGGTTT‑G‑GAGTAGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                                    <  1:194207/138‑1 (MQ=255)                                                                                                                                                                                           ctcccccCCCGGGGGGGGGGTT‑AGGGGTGAGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                                            <  2:158606/147‑1 (MQ=255)                                                                                                                                                                                           accccccTGGGGGGGGAGGGTG‑G‑GGGGGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                                    <  1:651882/138‑1 (MQ=255)                                                                                                                                                                                               ccgtttGGGGGGGGGGGA‑A‑GGGTGGAGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTGACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg                                                       <  2:476009/145‑1 (MQ=255)                                                                                                                                                                                                ccTCTGGGGCGCGGGCT‑A‑GGGTGAGGGG‑AACAGGCGGCATCGAAGCCAATCAGCCCCTATCGACCGCCTTTACGACTCAAATAACGATCAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCCTa                                                                        >  3:408577/1‑132 (MQ=255)                                                                                                                                                                                                ccTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                                               >  1:138016/1‑141 (MQ=255)                                                                                                                                                                                                ccTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGcc    >  4:430595/1‑200 (MQ=255)                                                                                                                                                                                                  tctGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGTTACAAAACCAAGAAATATCGCGGGAGGGGGCCGGATCGCCGGCGATATTAAGCAAGGTTTCGCCAg  >  4:163538/1‑200 (MQ=255)                                                                                                                                                                                                     ggggAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc                                                          >  1:377232/1‑141 (MQ=255)                                                                                                                                                                                                     ggggAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc                                                          >  1:128235/1‑141 (MQ=255)                                                                                                                                                                                                        gagaGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                                                >  3:147747/1‑132 (MQ=255)                                                                                                                                                                                                           agGGTT‑G‑GGGGGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                                                    <  1:305134/141‑1 (MQ=255)                                                                                                                                                                                                           |                                                                                                                                                                                                  CATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCGCAGGAG‑ATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCT‑CGCCCCTCTGGGGAGAGGGTT‑A‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAG  >  NC_000913/3608505‑3608893

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.