Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,687,5200AG49.1% 4.1 / 79.6 60G286G (GGT→GGC) fumAfumarase A
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base G (13/16);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAT  >  NC_000913/1687273‑1687749
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
ctTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGccc                                                                                                                                                                                                                                    >  1:308815/1‑251 (MQ=255)
    cATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAc                                                                                                                                                                                                                                >  1:164830/1‑251 (MQ=255)
    cATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAc                                                                                                                                                                                                                                >  1:165356/1‑251 (MQ=255)
     aTCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGGATTTACCACCAAACTGCGCGCCCAGAcc                                                                                                                                                                                                                               >  2:370120/1‑251 (MQ=255)
                aaGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTgcg                                                                                                                                                                                                                    >  1:455086/1‑251 (MQ=255)
                     cAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCg                                                                                                                                                                                                               <  1:379374/251‑1 (MQ=255)
                     cAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc                                                                                                                                                                                                                <  2:189707/250‑1 (MQ=255)
                       aCGCGCACCGCTTCGACTTCTATCGCTTTGCGAAGCTCTTCCGGGATATATTTGCCTGGCTTATGTTCCAGTTTTTCGATCCAGATCCCTTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGCGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGACTAACGCGGATGTCGCGCGCGAAGCATTTACCACCAAACCCCGCCCCCCGACCAAGCTTTTgtgcattcac                                                                                                                                                                                                             >  2:24447/1‑242 (MQ=21)
                                            cccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGAGCAGCAATTCTTTTTCCAg                                                                                                                                                                                           <  1:417966/248‑1 (MQ=255)
                                                   tGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                                                                                                                                                                 >  1:371123/1‑251 (MQ=255)
                                                      gcAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCgg                                                                                                                                                                              >  1:158384/1‑251 (MQ=255)
                                                        aGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATGCTTTTTCCAGTTCCACATCGCGGaa                                                                                                                                                                            <  1:383246/251‑1 (MQ=255)
                                                         gCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAAc                                                                                                                                                                           <  2:164830/251‑1 (MQ=255)
                                                                    atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTg                                                                                                                                                                >  1:174420/1‑251 (MQ=255)
                                                                         tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGacc                                                                                                                                                                           >  2:265326/1‑233 (MQ=255)
                                                                         tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAAc                                                                                                                                                                           <  1:265326/235‑1 (MQ=255)
                                                                          ttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                                                                                                                                                           <  1:363525/250‑1 (MQ=255)
                                                                           tGCCTGGCTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTTCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                                                                                                                                                         >  2:184264/1‑251 (MQ=255)
                                                                                 gagtATGTTCCAGTTTTTCTAGCCAGATCCTCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGGCGCCCATACCGACCGGTCAGGATGCGCCGGGACGTGGCAGGCGAATCACGCCGATGTCGTGAGCGAAGTATTTACCACCAAACGGCGCGCCCAGACCAATTTTTGGCGTTTCTATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                                                                                                                   <  1:483012/248‑1 (MQ=255)
                                                                                     aTGTTCCAGTTTTTCGATCCAGATCCCCTGTCGGTTGATCTTTGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATAGCGACCGGACAGGATGCCCCATGAAGTGGCAGGCGAATCACTCGGATGCCGTGAGCGAAGTATTTACCACCAAACTAAGCACCCAGACCAAGATTTTGCGCTTAGATCAGCAATTCTTTTTCCAGTTCCAAATCGCGGAACGCCTGACCGTGCTCATTTCCTTCCGTTg                                                                                                                                               >  1:320939/1‑251 (MQ=255)
                                                                                                tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatc                                                                                                                                     <  2:127013/250‑1 (MQ=255)
                                                                                                   cGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTcatcat                                                                                                                                   >  1:267657/1‑249 (MQ=255)
                                                                                                    gatcccgatcccATGAAGGTTGATCTTCGCTTTGATATTAAGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGCTGTCGTGAGAGAAGTCTTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCCATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatactcgt                                                                                                                                >  2:218805/1‑245 (MQ=37)
                                                                                                          gATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACTCCTGACCGTGCTCATTCCCTACCGTTGGCAGTTCATCATAGTATTTcg                                                                                                                          <  2:398444/251‑1 (MQ=255)
                                                                                                                cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGCTTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGa                                                                                                                      >  1:316348/1‑249 (MQ=255)
                                                                                                                     ggTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGt                                                                                                               >  1:248274/1‑251 (MQ=255)
                                                                                                                           tCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCCGACCAAGATTTTGCGCTTCGATAAGCAATTCTTTTTCCAGTTCCACATCGCGGCACGACCGACCCTGCTCATTCCCCTCCGTTGGCCTCTCCTCATAGTATTTTGCGGAAGCCCGTTTCAcc                                                                                                         >  2:169243/1‑251 (MQ=25)
                                                                                                                                  tttGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAg                                                                                                                                                                                           <  2:207456/162‑1 (MQ=255)
                                                                                                                                  tttGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAg                                                                                                                                                                                           >  1:207456/1‑162 (MQ=255)
                                                                                                                                        atTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCAcc                                                                                                         <  2:58348/238‑1 (MQ=255)
                                                                                                                                        atTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCAcc                                                                                                         >  1:58348/1‑238 (MQ=255)
                                                                                                                                        atTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCg                                                                                            >  1:18161/1‑251 (MQ=255)
                                                                                                                                          tACGGTCAGCAGAGCAGGAGCCGCCCATACCGACCGGCCAACATGCGCACTGACGTGGCAGACGAATCACGCGGATCTCGAGAGCCAAGAATTTACCACCAAACTGCGCACACCGACCAAGATTTTGCGCTTCGATCAGCCATTCTTTTTCCAGTTCCACCTCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGCTCGtt                                                                                          >  1:28778/1‑251 (MQ=255)
                                                                                                                                             ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCt                                                                                       >  1:379523/1‑251 (MQ=255)
                                                                                                                                             ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGAGGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATGTCGCTGAAGCCAGTTTCACCGTTTTAAGTTTCGTTTCt                                                                                       <  2:79245/251‑1 (MQ=255)
                                                                                                                                              gTCAGCAGAGCAGGAGACGCCCAGACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTg                                                                                      <  1:169243/251‑1 (MQ=255)
                                                                                                                                                 agcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTGCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATATGATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAg                                                                                   <  1:319167/251‑1 (MQ=255)
                                                                                                                                                   cagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTTCAGa                                                                                  >  1:96464/1‑250 (MQ=255)
                                                                                                                                                    agagcagGAGACGCCCATACCGACCGGGCATGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGTCCGTGCTCATTCCCTTCCTGTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAg                                                                                <  1:218805/251‑1 (MQ=255)
                                                                                                                                                       gcagGAGACGCCCATACCGACCGGGCAGTATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGATCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCGCCATCGCGGAACGCCTTTCCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTTTTTCGTGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTc                                                                             <  1:461218/251‑1 (MQ=255)
                                                                                                                                                              aCGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAAt                                                                      >  1:355595/1‑251 (MQ=255)
                                                                                                                                                                 cccATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATaac                                                                   <  1:228831/251‑1 (MQ=255)
                                                                                                                                                                  ccATACCGACCGGGCAGGATGCGCCGTGCCGTGGCAGGCGACTCACGCGGATGTAGCGAGCGAAGTATTTTCCACCAAACCGCGCACCCAGACCCAGATTTTGCACCTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATaacg                                                                  >  2:460183/1‑251 (MQ=255)
                                                                                                                                                                  ccATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGGCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATaacg                                                                  <  2:426214/251‑1 (MQ=255)
                                                                                                                                                                      accgaccgGGCAGCATGCGACGTGACGCGGCAGGCGAATCAAGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCCGACCACGATTTTGCGCTTCGATCCGCACTTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCCTAGTCTTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATaacgaacg                                                              >  2:474024/1‑251 (MQ=255)
                                                                                                                                                                         gaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCa                                                            >  1:435680/1‑250 (MQ=255)
                                                                                                                                                                            cgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAatat                                                        <  1:370120/251‑1 (MQ=255)
                                                                                                                                                                                cAGGATGCGCCGTGACGCTGCAGGCGAATCACGCGGAGGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATa                                                    <  2:23895/251‑1 (MQ=255)
                                                                                                                                                                                 aGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGt                                                                                               >  1:386401/1‑207 (MQ=255)
                                                                                                                                                                                 aGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGt                                                                                               <  2:386401/207‑1 (MQ=255)
                                                                                                                                                                                               cGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTGTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAATCCAGTTTCACCGTTTTAAGTTTCGTGTCTGCAGAAGTTCCACCAAGAACGAACGCAAGATGATACTGAGGAGAGgccgc                                     <  2:432028/251‑1 (MQ=255)
                                                                                                                                                                                                gTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCCCCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGgccccc                                    >  2:377927/1‑248 (MQ=255)
                                                                                                                                                                                                    cAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTAc                                <  2:316348/251‑1 (MQ=255)
                                                                                                                                                                                                      ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAgg                                         <  2:185428/240‑1 (MQ=255)
                                                                                                                                                                                                      ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAgg                                         >  1:185428/1‑240 (MQ=255)
                                                                                                                                                                                                         gAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGc                           >  1:466599/1‑251 (MQ=255)
                                                                                                                                                                                                               cgcgGATGGCGTGAGCGAAGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTgc                       <  2:374701/249‑1 (MQ=255)
                                                                                                                                                                                                                gcgGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTgcgc                     >  1:377632/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aGTATTTACCACCAAACTGCGCGCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTaa   <  2:248274/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGTATCTTCTCAACCAGGTAAt  >  2:66826/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
CTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAT  >  NC_000913/1687273‑1687749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑