Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009133,745,3710T→G16.9% 10.0 ‑1.7 0.0320W190G (TGG→GGG) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  new base (31/23):  ref base (125/141):  total (156/164)
Fisher's exact test strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 353 total aligned reads displayed.

CCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATT  >  NC_000913/3745121‑3745610
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
 cTTTCTGGTTCGGCCGTGGCGGATCCCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCtggtg                                                                                                                                                                                                                                                   <  2:731093/251‑1 (MQ=255)
  tttCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCCTGGTGATGGGCGCGACGCTGATGACTTCCtggtgg                                                                                                                                                                                                                                                  >  2:2439289/1‑251 (MQ=255)
                 ttGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGc                                                                                                                                                                                                                                 <  1:3212730/251‑1 (MQ=255)
                        tACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAAtt                                                                                                                                                                                                                          <  2:5554031/251‑1 (MQ=255)
                             ccgccgTTGCCGCTCTTCTGGTGCCGATGATGCGCAGCGCCAACTAGCCGTTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTAccc                                                                                                                                                                                                                     <  1:5244311/251‑1 (MQ=255)
                                ccgTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCtggtg                                                                                                                                                                                                                                                   <  2:1733155/220‑1 (MQ=255)
                                  gTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGGCTTTCTCCTCGGATGGTTGTTGGCGCAAAGCTGGTGGTTTTCCGGTGGGGGCCGGGCC‑CGTCCCCTACATTT‑CCCCGcc                                                                                                                                                                                                                 >  2:1206001/1‑250 (MQ=255)
                                      ccGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCGTTTCCCCCGGCATGATGATGGGGGCAACGCCGATGCGTTCCTGGTGGTGG‑CAGGCG‑CGTCGCCTTTATTTACCCCGCCAGcca                                                                                                                                                                                                            >  2:589293/1‑249 (MQ=255)
                                      ccGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCtggtggt                                                                                                                                                                                                                                                 <  2:4375825/216‑1 (MQ=255)
                                        gCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGCTTATCCCTCAGTAAGCTGTTTATGGCGGGCATTTCCCCCTGCACGACGTTTGGCGCAACGCTGATGCTTTCCTGGTGGTGG‑CACGCC‑CGTCGCCCTTATTTACCCCGCCAGCcaaa                                                                                                                                                                                                          >  1:5166995/1‑251 (MQ=255)
                                             aCTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCTGCATGATGATTGGCGCAACGCTGATGCTTTCCTGGTGGTGG‑CAGGCC‑CGGCGCCTTAAtttttcccccccgccaacaacaa                                                                                                                                                                                                     >  2:5416397/1‑231 (MQ=255)
                                                                 gcgcCAGCTATCTGGTCACCCGGGCGGCGGGGTTGATTCCCTCTGGCGGCATTAGCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGGGGATTATCCACCAGTAAGCGGTTTATGGCGGGCATTGCCCCGGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATc                                                                                                                                                                                       <  1:2331416/245‑1 (MQ=255)
                                                                 gcgcCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCCTTATGATGGGGGCAACGCTGATGCTTACCTGGTGGGGG‑GCGGCC‑GGTCGCCTTTATTTAACCCCCCAGCCAAAAGCAAcacagcaagaaaacttgcc                                                                                                                                                                                  >  2:2294223/1‑232 (MQ=255)
                                                                 gcgcCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGcc                                                                                                                                                                                                                                <  1:1366565/204‑1 (MQ=255)
                                                                    ccAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCGGGCTGATGGTGGGGGCCACCCTGGTGCTTTCCTGGGGGTGG‑CGgggc                                                                                                                                                                                                                                        >  2:4204964/1‑192 (MQ=255)
                                                                                   aCCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCGGCCTGGTGATGGGCGCCACCCTGCTGCTTTCCTGGGGGGGG‑CCGGCC‑CGTCGc                                                                                                                                                                                                                                 >  1:1016276/1‑185 (MQ=255)
                                                                                           gcggGGCTGATTGCCTTTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCGTCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGGTGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAAttt                                                                                                                                                                                                                         <  1:2577991/185‑1 (MQ=255)
                                                                                              gggCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGGTGATGGGCGCAACGCTGATGGTTTCCTGGGGGTGG‑GAGGCC‑CGTCGCCTTAATTTACCCCCcc                                                                                                                                                                                                                 >  2:1853310/1‑190 (MQ=255)
                                                                                              gggCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGcc                                                                                                                                                                                                                 <  1:1853648/190‑1 (MQ=255)
                                                                                                           ctGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGGATTGCCCCCGGCATGATGATGGGCGCAAAGCTGATGCTTACCtgggggt                                                                                                                                                                                                                                                 >  1:2696674/1‑147 (MQ=255)
                                                                                                           ctGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCtggtggt                                                                                                                                                                                                                                                 <  2:2696235/147‑1 (MQ=255)
                                                                                                               cggcATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCCACGCTGATGCTTACCTGGTGGTGG‑GAGGCC‑AGGCGCCCTAATTTTCCCCCCCCGCAAAAAACAACGAAGCCGGGAAACTGGCACTCTTTTTTCTCCGGGATTTGGGCGGTGTttttttccgtgcatattt                                                                                                                                   >  1:4372631/1‑238 (MQ=255)
                                                                                                                 gcATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTAccc                                                                                                                                                                                                                     <  2:651389/167‑1 (MQ=255)
                                                                                                                 gcATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTAGGGCGGGCATTGCCCCTGGCATGATGATGGGCGAAACGCTGATGCTTCCCTGGGGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCatta‑t                                                                                                                                 <  2:5202575/250‑1 (MQ=255)
                                                                                                                              aaTTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTTCCTGGGGGGGG‑GAGGCC‑AGTCGCCCTTATTTTCCCCCCCCGCCAAAAGCAACGATGCCGGGAAACTGGCACTCTTTTTTCTCCGGCATCt                                                                                                                                                              >  1:932232/1‑209 (MQ=255)
                                                                                                                              aaTTATCCCCCCTCCCATCCCTTTTATTATCTTCGGCTTTCCCAGCGGATTATCCATCAGAAGGCTTTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTCCCGGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCatta‑t                                                                                                                                 <  1:5184285/237‑1 (MQ=255)
                                                                                                                               attattCCCCCTTCCATTCCGTTTGTTATCTGGGGCGTTTCCAGCGGATTATCCAACAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCAAAGCGCTGATGCTTGCCTGGGGGTGG‑CAGGCC‑GGTCACCTTAATTTACCCCGCCAGCAAAAAGc                                                                                                                                                                                                       <  1:5886269/167‑1 (MQ=255)
                                                                                                                                ttattCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGGATGATGATGGGGGCAAACCTGGTGCTTTCCTGGTGGGGG‑CAGGcc                                                                                                                                                                                                                                        >  2:2776187/1‑134 (MQ=255)
                                                                                                                                ttattCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGGGG‑CAGGCC‑AGGCGCCCTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTTTCTCCGGGATCTGGGCGCTGTTTCTTCCCGTCATCATTT‑TTGGCggggt                                                                                                                        >  2:2979522/1‑242 (MQ=255)
                                                                                                                                         cTTCCATTCCGTTTATTATCTTCGGCTTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCTTTATTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTggg                                                                                                                                                           <  2:920343/201‑1 (MQ=255)
                                                                                                                                         cTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGGGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCg                                                                                                                                                                                                                                  >  2:3497144/1‑130 (MQ=255)
                                                                                                                                         cTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCg                                                                                                                                                                                                                                  <  1:2837600/130‑1 (MQ=255)
                                                                                                                                         cTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCg                                                                                                                                                                                                                                  <  1:3497711/130‑1 (MQ=255)
                                                                                                                                              aTTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCGGCATGATGATGGGGGCAAAGCTGATGCTTACCTTGTGGTGG‑GAGGCC‑AGTCGCCCTAATTTTCCCCCCCCGCCAAAAACAACGATGCCGGAAAACTGGCACCCTTTTTTCCCCGGGATCTGGGGGGCGTTTTTTTCAGGCATCATTA‑TTTGCGGGgttttgttctcgcccttttttc                                                                                                    >  2:3279007/1‑231 (MQ=255)
                                                                                                                                                                   gTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGGGG‑GAGGCC‑AGTCGCCTTAATTTTCCCCGCCCGCAAAAAACCACGATGCAGGGAATCTGGCACTCTTTTGTCTCCCGTATCTGGGCGCTGTTTCTTCCCGCCATCATTA‑TTGGCgggggtttttctcgccctctttt                                                                                                      >  1:1476862/1‑207 (MQ=255)
                                                                                                                                                                     ttCCAGCGGATTACCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGttt                                                                                                                       <  1:4314538/208‑1 (MQ=255)
                                                                                                                                                                       ccAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCCTTGCCCCCGGCATGATGATGGGCGCCACGCTGATGCTTACCTGGGGGGGG‑GAGGCC‑AGTCGCCCTAATTTACCCCGCCCGCCAAAAGCAACGGTGCCGGAAATCTGGCACTCTTTTTTCTCCCGTATCTGGgcg                                                                                                                                                         >  1:2414199/1‑173 (MQ=255)
                                                                                                                                                                            ggATTATCCATCAGTAAGCTGTTTATGGCGGGCCTTGCCCCCGGCATGATGATGGGCGCAACGCTGACCCCTACCTTGGGGGGG‑CCGGCC‑AGTCGCCTTACTTTACCCCCCCAGCCAAAAACCAAGATGCAGGGAATCTCGCCCTCCTTTTTCTCCGGGATCTggg                                                                                                                                                           >  1:2367132/1‑166 (MQ=255)
                                                                                                                                                                            ggATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGGGCAACGCTGTTGCTTACCTGGGGGGGG‑CAGGCC‑AGTCGGCTTTATTTACCCCGCCAGCCAAAAGCCACGATGCAGGGAATCTGGCACTCTTTTGTCTCCGGGATCTGGGCGCTGTTTCTTCCCGGCATAATTATTTGGGGGGGTTTGTTCCGCGCGTTTTTCCCCCCCCCAAACCGGGGcggctgtgggctttt                                                                      >  1:4546536/1‑237 (MQ=255)
                                                                                                                                                                            ggATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCCACGCTGATGCTTACCTTGGGGGGG‑CAGGCC‑AGTCGCCTTAATTTACCCCCCCAGCAAAAAGCAACGATGCAGGGAATCTTGCCCTCTTTTTTTTCCGGTATCTCGgcg                                                                                                                                                         >  1:5833090/1‑168 (MQ=255)
                                                                                                                                                                            ggATTACCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGACGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTggg                                                                                                                                                           <  2:2366725/166‑1 (MQ=255)
                                                                                                                                                                             gATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTCCCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCAc                                                                                                                                                                                 <  2:4342740/143‑1 (MQ=255)
                                                                                                                                                                               ttATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTTGCCTGATGATGGGGGCAACGCTGATGCTTTCCTGGGGGGGG‑CAGGCCAAGCCGCC‑TAATTTTCCCCCCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCttttttt                                                                                                                                                                        >  2:5177268/1‑147 (MQ=255)
                                                                                                                                                                                 aTCCATCAGTAGGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTGCCCCGCCAGCAAAAAGCAACGATGCAGGaaa                                                                                                                                                                                         <  1:5100279/131‑1 (MQ=255)
                                                                                                                                                                                 aTCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGG‑CAGGCC‑AGTCGCCTTAAtt                                                                                                                                                                                                                          >  2:975686/1‑98 (MQ=255)
                                                                                                                                                                                  tCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAg                                                                                                                                            <  1:3685735/175‑1 (MQ=255)
                                                                                                                                                                                  tCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGGGGAACGCCTATGGTTTCCCGGGGGGGG‑CCGGCC‑CGGCGCCTTTATTTAACCCGCCCGCCAAAAGAAAAGAAGCAGGGAATCTGGCACTCCTTTGTCTCCGGGATCTTGGCGCTGTTTTTTCCAg                                                                                                                                            >  2:3685133/1‑175 (MQ=255)
                                                                                                                                                                                      tCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGAGGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTcc                                                                                                                                              <  1:4066392/169‑1 (MQ=255)
                                                                                                                                                                                           aaGCTGTTTATGGCGGGCATTTCCCCCGGCATGATGATGGGCGCAACGCCGATGCTTACCTGGGGGGGG‑GAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCCGGGAACCTGGCAc                                                                                                                                                                                 >  1:5566909/1‑129 (MQ=255)
                                                                                                                                                                                             gCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGGGCAACGCTGATGCTTACCtggtgg                                                                                                                                                                                                                                                  >  1:3984002/1‑64 (MQ=255)
                                                                                                                                                                                             gCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCtggtgg                                                                                                                                                                                                                                                  <  2:3983343/64‑1 (MQ=255)
                                                                                                                                                                                                gTTTATGGCGGGCATTTCCCCCGGCATGGTGATGGGGGCAACGCCGATGCTTACCTGGGGGGGG‑CAGGCC‑CGTCGCCTTAATTTACCCCCCCAGCAAAAAACAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGCATCTGGgcg                                                                                                                                                         >  2:4605899/1‑148 (MQ=255)
                                                                                                                                                                                                 tttATGGCGGGCATTGCCCTTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAg                                                 <  1:2980008/250‑1 (MQ=255)
                                                                                                                                                                                                       gCGGGCATTGCCCCTGGCTTGATGATGGGCGCACCGCGTATTCTTACCTGGTGGGGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTAttt                                           <  1:1144097/250‑1 (MQ=255)
                                                                                                                                                                                                        cGGGCTTTGCCCCTGGCATGATGTTGGCCGCAACGCTGATGCTTCCCTGGTGGTGG‑CAGGCC‑AGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTAc                                         <  1:4422639/251‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGcc                                                                                                             <  2:3078381/175‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:4990964/1‑132 (MQ=255)
                                                                                                                                                                                                                gCCCCTGGCATGATGATGGGCGAAACGCTGATGCTTACTTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTc                                                                                                                  <  1:4210158/170‑1 (MQ=255)
                                                                                                                                                                                                                gCCCCGGGCATGATGAGGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CGGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGttt                                                           <  2:788106/225‑1 (MQ=255)
                                                                                                                                                                                                                 ccccTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:1295421/1‑131 (MQ=255)
                                                                                                                                                                                                                 ccccTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:2214945/1‑131 (MQ=255)
                                                                                                                                                                                                                 ccccTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:643453/1‑131 (MQ=255)
                                                                                                                                                                                                                 ccccGGGCATGATGATGGGCGCAACGCTGATGCTTCCCTGGTGGTGG‑AAGGCC‑AGCCGCCTTATTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTAc                                         <  1:3042637/242‑1 (MQ=255)
                                                                                                                                                                                                                     tGGCAGGATGATGGGCGCACCGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAg                                                                                                                                            <  1:2458907/140‑1 (MQ=255)
                                                                                                                                                                                                                      ggCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTtgc                                                                              <  1:4674592/200‑1 (MQ=255)
                                                                                                                                                                                                                      ggCATGATGAGGGGCGCAGCGCTGATGTTTACCTGGTGGTGG‑CGGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTAc                                         <  1:1828041/237‑1 (MQ=255)
                                                                                                                                                                                                                        cATGATGATGGGCGCACCGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTttcttc                                                                                                                                               <  2:4009682/134‑1 (MQ=255)
                                                                                                                                                                                                                        cATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTttcttc                                                                                                                                               >  1:4010341/1‑134 (MQ=255)
                                                                                                                                                                                                                           gatgatgGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtct                                                                                                                                                                       >  2:5642261/1‑107 (MQ=255)
                                                                                                                                                                                                                                   gcgcAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTTATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:152154/1‑113 (MQ=255)
                                                                                                                                                                                                                                   gcgcAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtct                                                                                                                                                                       >  1:769970/1‑99 (MQ=255)
                                                                                                                                                                                                                                   gcgcAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:5320624/1‑113 (MQ=255)
                                                                                                                                                                                                                                   gcgcAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:4128313/1‑113 (MQ=255)
                                                                                                                                                                                                                                   gcgcAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:4097507/1‑113 (MQ=255)
                                                                                                                                                                                                                                   gcgcAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTcc                                                                                                                                              >  1:1591174/1‑124 (MQ=255)
                                                                                                                                                                                                                                   gcccAACGCTGATGCTTCCCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCg                                                                                              <  2:5319726/168‑1 (MQ=255)
                                                                                                                                                                                                                                   ccgcACCGCTGATTCTTCCCTGGTGGTGG‑CAGGCC‑AGTCGCCTTATTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCacac                   <  2:4562022/245‑1 (MQ=255)
                                                                                                                                                                                                                                      caccGCGGAGGCTGACCGGGGGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCAt           <  2:5832094/248‑1 (MQ=255)
                                                                                                                                                                                                                                        aCGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTttcttc                                                                                                                                               >  1:4634590/1‑118 (MQ=255)
                                                                                                                                                                                                                                         cGCTGATGCTTCCCTGGGGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTAcccc                                                                                                   <  1:4575339/160‑1 (MQ=255)
                                                                                                                                                                                                                                         cGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:2500928/1‑107 (MQ=255)
                                                                                                                                                                                                                                         cGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:4735543/1‑107 (MQ=255)
                                                                                                                                                                                                                                         cGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:860099/1‑107 (MQ=255)
                                                                                                                                                                                                                                         cGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGtc                                                                                                                                          <  1:5511980/122‑1 (MQ=255)
                                                                                                                                                                                                                                             gATGCTTACCTGGTGGTGG‑CGGGCC‑GGTCGCCTAAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCg                                                                                              <  1:339812/161‑1 (MQ=255)
                                                                                                                                                                                                                                             gATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:339737/1‑103 (MQ=255)
                                                                                                                                                                                                                                             gATGCTTACCTGGTGGTGG‑CAGGCA‑AGTCCCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGttt                                                                                                        <  2:2404491/151‑1 (MQ=255)
                                                                                                                                                                                                                                              gtgtttcccgggtggggg‑CAGCAC‑AGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTAc                                         <  2:5117986/203‑1 (MQ=255)
                                                                                                                                                                                                                                              aTGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:976435/1‑102 (MQ=255)
                                                                                                                                                                                                                                               tGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:3323992/1‑101 (MQ=255)
                                                                                                                                                                                                                                               tGCTAACCGGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGtca                                                                                                                                         <  1:2257062/117‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttcccTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTc                                                                                                                  <  2:4221394/133‑1 (MQ=255)
                                                                                                                                                                                                                                                   tcccTGGTGGTGG‑CAGGCC‑AGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTg                                                                                                                               <  2:4638047/120‑1 (MQ=255)
                                                                                                                                                                                                                                                   tACCTGGTGGTGG‑CAGGCC‑GGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:149314/1‑97 (MQ=255)
                                                                                                                                                                                                                                                   tACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTg                                                                                                                               <  2:2105044/122‑1 (MQ=255)
                                                                                                                                                                                                                                                   tACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  1:2105418/1‑97 (MQ=255)
                                                                                                                                                                                                                                                     ccTGGTGGTGG‑CAGGCC‑AGTCCCCTTATTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCg                                                                                              <  1:4936383/153‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggtggtgg‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                         >  2:4577317/1‑92 (MQ=255)
                                                                                                                                                                                                                                                        ggtggtgg‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGttt                                                                                                                                                   >  2:1656109/1‑98 (MQ=255)
                                                                                                                                                                                                                                                        ggtggtgg‑AAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGttt                                                                                                                                                   <  1:1656413/98‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
CCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGG‑CAGGCC‑AGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTA‑TTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATT  >  NC_000913/3745121‑3745610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ 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.