Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,2190TG33.8% 104.8 / 12.8 72F86V (TTT→GTT) adeDadenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (13/34);  new base G (23/1);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 2.12e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCTATTAACCATAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTG  >  NC_000913/3843971‑3844468
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
atTCTATTAACCATAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCCATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCCACGGCGGGGCCCGGGttt                                                                                                                                                                                                                                                             >  2:42434/1‑251 (MQ=255)
         aCCATAATTTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTCTTCTCGACCTGATCAATGGCGGTGAAATTTCCGGGCCAATTGGTATTACAGGACGTTACCTTGCCGGTGTTGGCGCCGAATACACTGATGCTCCGGCTTTGCCGCGGATTGATGCTCGCGGCGCCACGGCGGTGGCGGGGttttgttatgtc                                                                                                                                                                                                                                                    >  2:518227/1‑242 (MQ=255)
            aTAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTGCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCAc                                                                                                                                                                                                                                                 <  2:321669/251‑1 (MQ=255)
              aaaTTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCt                                                                                                                                                                                                                                               <  2:380946/251‑1 (MQ=255)
                        acaTTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGGTGCTCACCCGCATATTGa                                                                                                                                                                                                                                     >  1:182341/1‑250 (MQ=255)
                             aGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGAGGAGAAATTTCCGGCCCCATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGTTTATTGATGTTCACCTGCATCTTGAATcc                                                                                                                                                                                                                                 >  1:368714/1‑249 (MQ=255)
                                               gAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCAc                                                                                                                                                                                                          <  1:114160/251‑1 (MQ=255)
                                               gAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCAc                                                                                                                                                                                                          <  2:288950/251‑1 (MQ=255)
                                                 aTTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTCTATTATTGCTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGCTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGGCAAGGTTTTTTGATGCTCACCTCCATATTTCATTCCACATGATGACGCCgggcgcct                                                                                                                                                                                                          >  1:219853/1‑245 (MQ=255)
                                                    gTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACGCTGATGGTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAGCGGCGGTGGCAGGGTTTATGGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTg                                                                                                                                                                                                     <  1:351088/251‑1 (MQ=255)
                                                                    ggCGAAGCTGTTGCCGCTTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTTCCAGGGGTTTTTGATGTTCACCCGCCAATTTAATCCCACATGGTGGCGCCGGTGCACtttttaaccggttcccttcc                                                                                                                                                                                      >  2:352971/1‑235 (MQ=255)
                                                                     gCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTATAGACCTGATCAATGGCGGAGAACTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGCATACACTGATGCTCCGGCTTTGCAGCGGATTGACGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGCTGCTCACCTGCATATTGAATCCAGCATGATGAAGCCGGCCACTTTTGAAACCGCCACCCTGccc                                                                                                                                                                                    >  1:392901/1‑250 (MQ=255)
                                                                          gCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGACATTTCCGGCCCACTTGCGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCAGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGGCACGGGTTTTTGATTCTCACCGTCATATTGTATCCAGCCATATGAAGCCGGCTCCCTTTTAAACCGCCACCCCGCCGCCCgg                                                                                                                                                                               >  2:492686/1‑251 (MQ=255)
                                                                          gCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATTCTCACCCGCATATTGAATCCAGCATGATGACGCCGGGCACTTTTTAAACCGCTACCCTGCCGCGCgg                                                                                                                                                                               >  1:15299/1‑251 (MQ=255)
                                                                                gCCGATTATATTATTGATAATGTCGCTATGCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAAGACACTGATGCTCCGGCTTGGCAGCGGATTGATGCTAGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTTACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTggc                                                                                                                                                                         <  2:439918/251‑3 (MQ=255)
                                                                                 ccGATTATATTATTGATAATGTCTCTATTCTCTACCTGAGCAATGGCGGAGAAATTTCCGGCCCAAGGGGGAGGAAAGGACGGGACAGGGCCGGGGGGGGCGCAGAAGACACGGAGGCGCCGGCGGGGCAGCGGAGGGAGGCGCGCGGCGCAACGGCGGGGCCAGGGGGGAGGGAGGCGCACCGGCAGAGGGAAGCCAGCAGGAGGACGCCGGGCACTTTTGAACCCGCTACCCTGCCGCGCGTCCTgacg                                                                                                                                                                        <  1:328880/251‑1 (MQ=255)
                                                                                    attatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATTCTCAACTGCCTATTTAATCCAGCATGATGACGCCGGTCCCTTTTTAACACGCTACCCTTCCGCGCgggcccacgaac                                                                                                                                                                     >  1:363161/1‑241 (MQ=255)
                                                                                      tatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTg                                                                                                                                                                           <  1:448993/243‑1 (MQ=255)
                                                                                          ttattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCTGCCTATTGGATCCAGCCTGATTACGCCGGTCACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGGCCgctttt                                                                                                                                                               >  1:533974/1‑246 (MQ=255)
                                                                                           tattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTAtt                                                                                                                                                               <  1:418845/250‑1 (MQ=255)
                                                                                                aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGAc                                                                                                                                                         <  1:518227/251‑1 (MQ=255)
                                                                                                    tGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATTCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCCACCCCGCCGCGCGGCCTGACGACCGTTATTTGCGGCCCTc                                                                                                                                                     >  1:185267/1‑251 (MQ=255)
                                                                                                     gTCTCTATTCTCGACCTGCTCAATGGCGGAGAAATTTCCGGCCCACTTGTGATTACAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGGTTATTTATGGTACCCCGCCAATTTAATCCCGCGTGATGACGCCCGTCACATTTTAAACCGCCACCCTCCCGCGCCGCCTGACGAcccgtatttttgcccctcc                                                                                                                                                    >  2:468914/1‑234 (MQ=255)
                                                                                                        tctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaaa                                                                                                                                                                                                  <  1:444037/202‑1 (MQ=255)
                                                                                                         ctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCCATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGAGCAGAATAAACTGATGCTCCGGCTTTGCCGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTTTTGATGCTCACCTGCCTATTGAATCCAGCATGAATACGCCGGTCCCTTTTTAAACCGCTCCCCCGCCCCGCGGCCCGACGACCGTTATTTGCGGCCCTCATgg                                                                                                                                                 >  1:261623/1‑249 (MQ=255)
                                                                                                               ccGACCTGATCAATGGCGGAGAAATTTCCGGACCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTc                                                                                                                                          <  2:323514/250‑1 (MQ=255)
                                                                                                                cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCAc                                                                                                                                                                                                                                                 <  1:405825/151‑1 (MQ=255)
                                                                                                                 gACCTGATCAATGGCGGAGAACTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCCGCGGATTGATGCTCGCGGCCCAACGGCGGTGCCAGGGTTTATTTATGCTCACCCGCATATTGTATCCAGATTGATGACGCCGGTCACTTTTTAAAACCCTACCCTTCCGCGCGGCCTGAAGACCGTTTTTTGCGGCCCTCCTGCAATCGTCa                                                                                                                                         >  1:487298/1‑250 (MQ=255)
                                                                                                                 gACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAAGGGTTATTGAAGCTCAACTGCATCTTGAATACAGCATGATTACGGCCGTCCCCTTTGGAAACGCTACCCCTCCGCGTGGGCTGACGACCGGTATTTggggacctcaaggaatcggccc                                                                                                                                        >  1:587721/1‑230 (MQ=255)
                                                                                                                      gATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCTCTACCCTTCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGa                                                                                                                                   <  1:63536/251‑1 (MQ=255)
                                                                                                                      gATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTg                                                                                                                                    <  1:352971/250‑1 (MQ=255)
                                                                                                                             ggCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCg                                                                                                                             <  2:185267/250‑1 (MQ=255)
                                                                                                                             ggCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGCGCCCGGGGTTCTTGATGCCCACCTGCATATTGAATCCAGCATGATGACGCCCGGCACTTTTGTAACCGCCACCCTGCCCCGCGGCCCTACGCCCCTTCTTTGCGACCCTCCTGAAACCCTACACGTGATGGGCg                                                                                                                             >  2:165013/1‑250 (MQ=255)
                                                                                                                              gCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACAATGCCGGTGTTGGCGCAGAATACACTGATGCGCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTGATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTGGAAACCGCTACCCTGCCGCTCGGGCGGACGACCTTGATTGGGGTCCCTCATGAAATCGTCAACGTGAGGTGCGaa                                                                                                                           <  2:239220/251‑1 (MQ=255)
                                                                                                                                       tttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCCGAATACACTGATGCTCCGGCTTTGCAGCCGATTGATGCTCGCGGCGCACCGGCGGAACCCGGGATTATTGATGCTCACCTGCATATTGAATCCAGCATGACGACGCCGGTCACTTTTTAAACCCCTCCACCGCCGCGCGGCCTGACGACCGCTATTTGCGACCCTCCTGAAATCGTCAACGTGATTGGGGAAGCCGGAttt                                                                                                                  >  1:281776/1‑250 (MQ=255)
                                                                                                                                           cGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGTTCACCTGCATATTGCATCCAGCATGATGACGCCGGTCACTTTTTAAACCGCTACCCTGCCCCCCGGCCTGACGCCCGTTATTTGCGCCCCTCCTGAAAACGGCAACGGGATGGGGGAAGCCgggtttgccc                                                                                                              >  2:90726/1‑243 (MQ=255)
                                                                                                                                               ccAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCAcc                                                                                                                                                                                                                                                <  2:337218/121‑1 (MQ=255)
                                                                                                                                               ccAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGCGCCCGGGTTTTTTTATGCTCAcc                                                                                                                                                                                                                                                >  1:337218/1‑121 (MQ=255)
                                                                                                                                                  aTTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAAGCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCTGCCTGACGACCGTTATTTGCGACCCTCATGGAATCTTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGc                                                                                                       <  1:280709/251‑1 (MQ=255)
                                                                                                                                                     gtgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGccc                                                                                                     <  1:42434/250‑1 (MQ=255)
                                                                                                                                                      tgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGAGCCAGGGTTTATTGGTGCTCACCTGCATATTGAATCCAGCATGATGAAGCCGGACACTTTTTAAACCGCCACCCCGCCGCGCGGCCCGACGACCGGTATTTGCGACCCTCATGGAATCGTCAACGTGATGGGCGCAACCGGATTCGCCTGGTTTGccccc                                                                                                   >  1:325865/1‑249 (MQ=255)
                                                                                                                                                       gATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCGCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTGCCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGGCTGGTTGGCCCGCt                                                                                                  <  2:215689/251‑1 (MQ=255)
                                                                                                                                                         ttAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCCCCTGGATACTGAATACAGAATGACGACCCCGGGCCCTTTTGAAACCGCCAACCTGACGCGCGGCCTAACGGCCGTTATTTGCGAACCCCATGAACACGTCCACGCGATGGGAGAAGCCGGATTCGCCTGGTTTGCCCGCtgt                                                                                                >  1:42594/1‑251 (MQ=255)
                                                                                                                                                            aaGGACGTTACATTGACGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGGCAGGGGTTATTGATGCTCACCTGCATATTGAAACCAGCACGATGACGCCGGTCACTTTTGTAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCCTGAAATCGTCCCAGTGATGGGCGGAGCCGGCTTCGCCCGGTTTGCCCCCTGTGcc                                                                                             >  2:254798/1‑251 (MQ=255)
                                                                                                                                                               gACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGTTTATTGATGCTCACCTGCATATTTAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCCACCCTGCCGCGCGGCCCGACGGCCGCTATTTGCGGACCTCCTGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCAGGGTTGCCCCCTGTGCCGa                                                                                           >  2:202898/1‑250 (MQ=255)
                                                                                                                                                                         tGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGTTGCTCACCTGCATATTGAATCCAGCCTGATGACGCCGGTCACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTTTTTTCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTTCCCGCTGTGCCGAACaggcaaggg                                                                                >  1:233632/1‑250 (MQ=255)
                                                                                                                                                                          gCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCTAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaaggc                                                                                <  1:429473/250‑1 (MQ=255)
                                                                                                                                                                             ggTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATTAAATCGTCAACGTGATTGGCGGAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCaaa                                                                             >  1:90884/1‑250 (MQ=255)
                                                                                                                                                                               tgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAAcc                                                                          >  1:407874/1‑251 (MQ=255)
                                                                                                                                                                                      gcAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACtt                                                                   >  1:19707/1‑251 (MQ=255)
                                                                                                                                                                                               acTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGGGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTc                                                            >  1:489496/1‑249 (MQ=255)
                                                                                                                                                                                                  gATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCCACGGCGGTGCCAGGGGTTTTTGATGCTCACCTGCATATTGAATCCAGCATGATGCCCGCGGTCACTTTTGAAACCGCTACCCTGCCCCGCGGCCTGCCGACCGTTATTTGCGCCCCTCCTGAAATCGGCAACGTGATGGGGGGAGCCGGGTTCGCCTTGGTTGCCCGCTTTGCCGAAAAGGCGAGGCAAAACCAGTACTTACAGGGCAggcc                                                       >  2:95414/1‑248 (MQ=255)
                                                                                                                                                                                                  gATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGct                                                        <  1:118591/250‑1 (MQ=255)
                                                                                                                                                                                                          ggCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTAc                                               <  1:22783/251‑1 (MQ=255)
                                                                                                                                                                                                            cTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGGCCGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCCCTTTTGAAACCGCTACCCTGACGCGCGGCCTGACGACCGTTATTTTCCACCCTCATGAAATCGTTAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCCGTGCCGAACAGGCCAGGCAAAACCAGCACTTACAGGTCAGCTCTTGCGTAccc                                             >  1:426723/1‑251 (MQ=255)
                                                                                                                                                                                                                gCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCCCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCCTGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGGTTGCCCGCTGTGCCGAACAGGCAAGGCCAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGccc                                         >  1:238365/1‑251 (MQ=255)
                                                                                                                                                                                                                 cAGCGGATTGATGCTCGCGGCGCAACGGGGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGAGGACGCCGGTCACGTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTGATTTGCGACCCTCATGAAATCGGCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCt                                        <  2:533974/251‑1 (MQ=255)
                                                                                                                                                                                                                  aGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCCTGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGGGCCGAACAGGCAAGGCAAAATCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTg                                       >  1:437162/1‑251 (MQ=255)
                                                                                                                                                                                                                   gtggATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTTAAACCGCTTCCCTGCCGCGCGGCCGGACGACCGTGATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGGTGGCCCTCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGGCAGCTCTTGCGGTCCCGCCCTgg                                      <  2:249106/249‑1 (MQ=255)
                                                                                                                                                                                                                   gCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTgg                                      <  1:428357/251‑1 (MQ=255)
                                                                                                                                                                                                                         tGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCTGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACGAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGc                                 <  1:246697/250‑1 (MQ=255)
                                                                                                                                                                                                                          gATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGGTATTTGCGAACCTCATGAAAACGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTTCGTAACCGCCCTGGGAGGCTg                               >  1:270085/1‑251 (MQ=255)
                                                                                                                                                                                                                           atgtgcgcgGCGCAACGGCGGTGCCAGGGTTCATTGATTCTCTCCTGCAGATTGAAGCCAGCATGACGTCGCCGGTCACTTTTGAAATCGCTCCCCTGCCGCGCGGCCCGACGACCGTTATTTGCGACCCTCATGAAATCGTAAACGTGATGGGCCAAGCCGGATTCGCCGGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGGCAGCTCGTGCGTACCCGCCCTGGATGGCTGc                              <  2:499159/246‑1 (MQ=255)
                                                                                                                                                                                                                                   ggcgCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGtt                        <  1:502927/249‑1 (MQ=255)
                                                                                                                                                                                                                                   ggcgCAACGGCGGTGCCAGGGGTTATTTATGCTCACCTGCATATTGAATCCAGCCTGATGACGCCGGGCACTTTTTAAACCGCCACCCTGCCGCGCGGCCTGACGACCGTTATTTTCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGGCCGCTCTTGCGTACCCCCCCTGGGAGGCTGCGAGGtt                        >  2:502927/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cgcAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACTCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCTTGTTTTGCCCTCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACg                    <  2:182938/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaCGTCGGTGCCAGGGTTTATTGATGCTCGCCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGTATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTAGTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTg                 <  2:487298/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cggcggTGCCAGGGTTTATTGATGCTCACCTGCCTATTGAATCCAGCATGATGACGCCGGTCCCTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGGCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGGGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCCCTTGGGTACCCGCCCTGGAAGGCTGCGGTGTTCACGGTGcc               >  2:432944/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cggcggTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGcc               <  1:492686/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cggTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGc                                                                                                                                                                                      <  1:388817/81‑1 (MQ=255)
                                                                                                                                                                                                                                             cggTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGc                                                                                                                                                                                      >  2:388817/1‑81 (MQ=255)
                                                                                                                                                                                                                                                    aGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGGTCGTGCGTTCCCTCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTAccc     <  1:581069/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTAccc     >  2:483057/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        gTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGAGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTg  >  2:103029/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
ATTCTATTAACCATAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTG  >  NC_000913/3843971‑3844468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: