Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,2190TG31.0% 121.6 / 12.5 71F86V (TTT→GTT) adeDadenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (20/29);  new base G (22/0);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 6.05e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATTCTATTAACCATAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTG  >  NC_000913/3843970‑3844468
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
aatTCTATTAACCATAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGtt                                                                                                                                                                                                                                                             <  1:406961/251‑1 (MQ=255)
      aTTAACCATAAATTTCATCACAGTATCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCGGATTATATTATTGATAATGGCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAAGTGTGATTAAATGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAtt                                                                                                                                                                                                                                                        <  2:214215/249‑1 (MQ=255)
               aaaTTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCt                                                                                                                                                                                                                                             >  1:312628/1‑251 (MQ=255)
                 aTTTCCTCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGCTCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTACAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGCTGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCTGc                                                                                                                                                                                                                                           >  2:93486/1‑251 (MQ=255)
                  tttCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTTTTGGTGCTCACCTGc                                                                                                                                                                                                                                           >  1:14648/1‑250 (MQ=255)
                   ttCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGGTGCTCACCTGcct                                                                                                                                                                                                                                         >  1:390353/1‑249 (MQ=255)
                        cacaTTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTATTGATGCTCCCCTGCATATTg                                                                                                                                                                                                                                    >  1:506905/1‑250 (MQ=255)
                          caTTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATTCTCCCCCGCATATTTAAt                                                                                                                                                                                                                                 >  2:193110/1‑251 (MQ=255)
                          caTTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGTTTATTGTTGCTCCCCTGCATATTGAAt                                                                                                                                                                                                                                 >  1:452531/1‑251 (MQ=255)
                           aTTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGAGCCAGGGTTTATTGATGGTCACCTgacaaatgaat                                                                                                                                                                                                                                  >  1:386200/1‑240 (MQ=255)
                             tAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATcc                                                                                                                                                                                                                               <  1:328283/250‑1 (MQ=255)
                                       gAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatgatg                                                                                                                                                                                                                      <  2:390353/249‑1 (MQ=255)
                                           aCCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGGTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACGGATGCTCCGGCTGTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAAGCCAGCATGATGACGCCg                                                                                                                                                                                                                <  2:253736/251‑1 (MQ=255)
                                                        aGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaa                                                                                                                                                                                                   <  1:215099/250‑1 (MQ=255)
                                                         gCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAc                                                                                                                                                                                                 <  1:309077/251‑1 (MQ=255)
                                                                        gACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCGGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACGGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGAGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCGGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgca                                                                                                                                                                                  <  2:6204/251‑2 (MQ=255)
                                                                             tGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGTTTATTGGTGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCCCTTTTTAAACCGCTACCCTTCCGCGCGGcc                                                                                                                                                                             >  1:98084/1‑251 (MQ=255)
                                                                              tttGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTTATGATCACCTGCATATTTGATCCCGCATATTGTCGCCGGTCACTTTTGAAAACGCCACCCTTCcgctctgccc                                                                                                                                                                            >  2:206186/2‑244 (MQ=255)
                                                                                tGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTg                                                                                                                                                                           <  1:585074/250‑1 (MQ=255)
                                                                                  ccGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgacg                                                                                                                                                                        <  1:17045/251‑1 (MQ=255)
                                                                                     attatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGAcc                                                                                                                                                                     <  1:266542/251‑1 (MQ=255)
                                                                                            tattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTAtt                                                                                                                                                               <  1:243072/250‑1 (MQ=255)
                                                                                             attGATAATGTCTCTATTCTCGCCCTGATCAATGGCGGCGAAATTTCCGGCCCAGTTGTGATTAAAGGACGTTACCTTGAAGGTGTTGGCGCAGAATACACTGATGCTCCGGATTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCACGGGGTTTTTGGTGCTCACATGCATATTGAATTCCGCCCGGTGACGCCCGGCACTTTTTAAAACGCTACCCCGCCGCgccgccctcccgccggtatttt                                                                                                                                                             >  2:506359/1‑235 (MQ=255)
                                                                                                gATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCg                                                                                                                                                           <  2:77006/250‑1 (MQ=255)
                                                                                                  tactGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGCGCCCGGGGTTTTTGTTGTCCCCCTGCCTATTGAATTCCGAATGATGACGCCGGGCACCTTTTAAACCGCCACCCCTCCCCGCGGcccgcccgccgctattttcgatc                                                                                                                                                        >  2:4203/4‑230 (MQ=255)
                                                                                                   aaTGTCTCTATTCTCGACCTGATCAATTGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGAccc                                                                                                                                                       <  2:473790/251‑1 (MQ=255)
                                                                                                   aaTGTCTCTATTCTCGACCTGATCAATGGCCGAGAAATTTCCGGCCCAATTGTGAATAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGAcc                                                                                                                                                        <  2:452531/250‑1 (MQ=255)
                                                                                                      gTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTc                                                                                                                                                     <  2:65386/250‑1 (MQ=255)
                                                                                                         tctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATACAGCATGATGACGCCGGGCACTTTTGAAACCGCCACCCTGCCGCGCGGGCTGACGCCCGTTATTTGCGACCCTcctg                                                                                                                                                  >  1:355945/1‑247 (MQ=255)
                                                                                                                 cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCg                                                                                                                                            <  2:98084/248‑1 (MQ=255)
                                                                                                                      tGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCCCACCTGCATATTGAATTCAGCATGATGACGCCGGTCACTTTTGAAACCCCTAACCCGCCGCGCGGACTGGAGGCCGTTATTTGTGACCCTCATGaaaactgcaaagtg                                                                                                                                    >  1:353569/1‑240 (MQ=255)
                                                                                                                          cAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGGTTTTTGATGGTCACCTGAATTATGGAACCCGCATGATGGAGCCGGGCACTTTTTCAAACCCATCCCCGCCGCGCGGCCTGACGGCCGTTCTTTGCGGCCCCCATTGAATCGTTCACGTGATgg                                                                                                                                >  2:411234/1‑251 (MQ=255)
                                                                                                                                       aTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGGTGCTCCCCTGCATATTGCATCCCGCCTGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCCTGAAATCGTCCACGTGATGGGCGAAGCCGGAtt                                                                                                                   >  1:578182/1‑251 (MQ=255)
                                                                                                                                         ttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTc                                                                                                                  <  2:100373/250‑1 (MQ=255)
                                                                                                                                            cGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCt                                                                                                              <  1:421277/251‑1 (MQ=255)
                                                                                                                                             ggCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTgac                                                                                                                                                                         <  1:20410/191‑1 (MQ=255)
                                                                                                                                                     tgtgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTAGTTGTGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGACACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTTCGACCCTCATGAAATCGTCCACGTGATGGGAGAAGACGGATTCGCCTGGTTTgacc                                                                                                     >  1:449258/1‑248 (MQ=255)
                                                                                                                                                      gtgATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTTCCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGccc                                                                                                     <  2:578182/250‑1 (MQ=255)
                                                                                                                                                         aTTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCCCCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCtg                                                                                                 >  1:551696/1‑251 (MQ=255)
                                                                                                                                                         aTTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACCGCGGTGCCAGGGGTTATTGATGCTCACCTGCCTATTGAATCCAGCATGATGACGCCGGTCACTTTTTAAACCGCTCCCCTGCCGCGCGGCCTGGCGACCGTTATTTGCGACCCCCATGAAATCGTCCACGGGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGccg                                                                                                 >  1:370404/1‑249 (MQ=255)
                                                                                                                                                               ggACGTTACATTGCCGGTGTTGGCGCAGACTACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGAGCCAGGGGTTTTTGATGCTCACCTGCAAATTGAATCCAGCATGGTGACGCCGGTCAATTTTGAAACCGGTAACCTGCCGCCCGGCCTGGCGCACGGTATTTTCGACCCTCATGAAAATGTCCACGTGGAGGGGGAAGCCGGGTTCACCTGGGTTTCCCGCcgtggcc                                                                                            >  2:387834/1‑247 (MQ=255)
                                                                                                                                                               ggACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCTGCATATTGAATCCCGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCCGCCGCGCGGCCTGACGACCGTTATTTGCGCCCCTCATGAAATCGTCACCGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCgc                                                                                           >  2:126935/1‑250 (MQ=255)
                                                                                                                                                                     tACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGAAACGGAGGTCCCAGGGTTTATTGATGCTCACCTGCATATTGAAACCAGCCTCATGACGCCGGGCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTTTTTGCGACCCTCCTGGAATCGTCCACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAAAaggc                                                                                     >  1:538824/1‑251 (MQ=255)
                                                                                                                                                                       cATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCATGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaa                                                                                   <  2:227609/251‑1 (MQ=255)
                                                                                                                                                                          tGCCGGTGTTGGCGCAGAATACACTGATGCTCCGCCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCCTATTGAATCCAGCATGCTGGCGCCGGTCACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGGGATGGGCGAAGCCGGATTCGCCTGGGTTGCCCGCTGGGCCGAACaggcaaggc                                                                                >  2:383466/1‑251 (MQ=255)
                                                                                                                                                                             cGGTGGTGGAGCAGAATACACTGGTGCTCCCGCTTTGCAGCTGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGTTGCTCACCTGCATATTGACTCCAGCATGATGACGCCGGTCCCTGTTGAAGCCGCTACCCTGCCGCGCGGCGTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCGGCTGTGCCGACCAGGCAAGGCaaa                                                                             <  2:129296/251‑1 (MQ=255)
                                                                                                                                                                                         aGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGAGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCCTATAGAATCCAGAATGATGACGCCGGTCTCTTTTGAAAACGCTACCCTCTCGCGCGGCCCGCCGACCCTTATTTGCGACCCTCCTGAACTCGACCACGGCATGAGGGAAGCCGGATTCCCCTGGTTTCCCCCCTGTGCCGGAAAGGAAAAGCAAAAACAGTACTTAc                                                                 >  1:352636/1‑251 (MQ=255)
                                                                                                                                                                                         aGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTAGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCCCCTGCATAGTGACTCCAGCAAGCAGACTCCGGGCACTTTTTAAACCGCTTCCCTGCCGCGCTGCCCAACGACCGGTAATTGCGGACCTCATGAAATAGGCAACGGGATGGGCGACACCGGCTGCGCCTGGTTTGCCCGGTGAGCCGGACAGGCAAGGCCAAAACAGTACttcc                                                                 >  2:192800/1‑249 (MQ=255)
                                                                                                                                                                                              acacTGCTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCCTGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTc                                                            >  1:468368/1‑251 (MQ=255)
                                                                                                                                                                                                   gATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGTCTGTGCCAGGGTTTATTGAGGCTCACCTGCATATTGAATCCAGCACGAGGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCTCGGCCTGACGACCGTTATTTGCGACCCGCATGAAATCGTCAACGTGAGGGGCGAAGCCGGATTCGTCTGTTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGc                                                         <  2:567547/249‑1 (MQ=255)
                                                                                                                                                                                                   gATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGctc                                                       >  1:242497/1‑251 (MQ=255)
                                                                                                                                                                                                     tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAg                                                                        <  2:21535/232‑1 (MQ=255)
                                                                                                                                                                                                     tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGGGAAGCCGGGTTCGCCCGGTTTGCCCGCTGGGCCGGACAGGCAAGGCAAAAACAg                                                                        >  1:21535/1‑232 (MQ=255)
                                                                                                                                                                                                      gCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTTCATATTGCATACAGCATGATGAAGCCGGTCCCTTTTGAAAACGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATAGTCAACGTGATGGGCGGAGCCGGATTCGCCTGGTTTGGCCGCTGTGCCGAACAGGCAAGGCAAAACCCGTACTTACAGGTCAGCTCTTg                                                    >  1:37497/1‑251 (MQ=255)
                                                                                                                                                                                                      gCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTg                                                    <  2:506905/251‑1 (MQ=255)
                                                                                                                                                                                                      gCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACACTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCtt                                                     <  1:93486/250‑1 (MQ=255)
                                                                                                                                                                                                          tggCTTTGCAGCGGATTGATGCTAGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGGTGGCGCCGGTCCCTTTTTAAACCGCTACCCTTCCGCGCGGCCTGACGACCGTTATTTTCGACCCTCATGAAATCGTCAACGGGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGGTGTGCCGAACAGGCCAGGCAAAACCAGTACTTACAGGTAAGCTCTTGCGt                                                 >  2:366937/2‑250 (MQ=255)
                                                                                                                                                                                                           ggCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGGAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTAc                                               >  1:122273/1‑251 (MQ=255)
                                                                                                                                                                                                                   aGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGGCACTTTTGAAAACGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCCCCCTg                                       >  1:220327/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGa                                     <  1:383466/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCCTATTGAATCCAGCATGATGACGCCGGTCACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAAAGTGATGGGCGAAGCCGGATTTGCCTGGTTTGCCCGCTGTGCCGGACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCCGGAAGGCTGCgg                            >  2:7868/1‑250 (MQ=255)
                                                                                                                                                                                                                               cTCGCGGCGCAACGGCGGAGCCAGGGGTTTTTTACGCTCCCCCGCCTATTGAATCCAGCATGATGACGCCCGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCACCGGGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACCGGCAAGGCAAAACCAGTACTTACAGGTCAGCGCTTGCTTACCCGCCCCGGAAGGCTGCGa                            >  2:383469/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ggcggTGCCAGGGTTTATTGATGCTCACCTGCATATTGCATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCt                                                                                                              >  2:402419/1‑155 (MQ=255)
                                                                                                                                                                                                                                            ggcggTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCt                                                                                                              <  1:402419/155‑1 (MQ=255)
                                                                                                                                                                                                                                             gcggTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGcc               >  1:162047/1‑249 (MQ=255)
                                                                                                                                                                                                                                                gTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGttt          <  1:438108/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGttt          >  1:474029/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 tGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGCAGGCTGCGATGTTATCGGTGCCTGtttt         >  2:314088/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGCTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGtttt         >  1:230126/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   ccAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTAc       >  1:34893/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ccAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAAGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGTTGACAGTTTTAc       >  1:282663/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     cAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGACCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTAcc      >  2:73960/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          tttATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGTTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGGGCCAGTTTTACCCTTg  >  2:173607/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
AATTCTATTAACCATAAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTG  >  NC_000913/3843970‑3844468

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