Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,422,263 G→A G186S (GGC→AGC)  yfcH → epimerase family protein YfcH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,422,2630GA93.6% 138.5 / ‑1.2 48G186S (GGC→AGC) yfcHepimerase family protein YfcH
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (26/18);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTG  >  NC_000913/2422021‑2422499
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
cACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCg                                                                                                                                                                                                                                      >  1:324025/1‑251 (MQ=255)
   cGTCGGTACTCATTTCCGGCTCGGAAACGGGCTATTATGGGGAATAAGGTGATGTGGTGTTTCCCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCGGCGTGCCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCTAtt                                                                                                                                                                                                                                   <  2:356545/251‑1 (MQ=255)
      cGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGc                                                                                                                                                                                                                                >  1:52507/1‑251 (MQ=255)
                   cgcttCGGCACCGGTCTATTTTGGTGACGTGGGTAAAGTGTTGGTTCCCGAAGAGGACCCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCGGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCGGAGCGGGCCGATTGGCTCTGGTCGGCAgg                                                                                                                                                                                                                   <  2:192640/247‑2 (MQ=255)
                             ccgggcaataagggatattaaGGTGAAGTGGTGGTTACCAAAGGGAACCCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAATAAATTGCCTGCCGGGCGCAAAGTGACAAACCGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGTCGGTATTCGCGTGAAAAGGCTGCCGCCGTGTCGTCTTGGCCTGAGCGGGCCGAGTGGCTCTGGTCGGCAGTATCTGGCCTg                                                                                                                                                                                                         <  2:224250/231‑1 (MQ=255)
                                      tatGGAGACTTAGGAGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCACAGTGACAAAACGCGAGTGAGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGACTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCCGTATCTGGCCTGGATTCatat                                                                                                                                                                                                >  1:349277/1‑251 (MQ=255)
                                                                   agAAGGGGAACCGCCGCATAACGAATTTACCCGTAACCTCTGCGCTCGCTGGGAAGTAATTGCCTGCCGGGCGCAAAGTGACAAACCGCGAGTGCGCCTGCTGCGTACGGGTGTAGTGCTGGCACCGGAGGGCGGTATTCGCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGgc                                                                                                                                                                   <  2:197345/250‑1 (MQ=255)
                                                                         ggAACCGCCGCATAACGAATTTACCCATAAACTCTGAGCTCGCTGGGAAGACATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTgg                                                                                                                                                             >  1:36797/1‑251 (MQ=255)
                                                                                       caCGAATTTATCAATAAACTCTGCGCTCGCTGGGCAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGACTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAACTGCTGCCGCCGCTTCGTCTTGGCCTGAGCGGTCCGGTTGGCACTGGTCGGCAGTATCTGGCCTGGTTTCCTATTGATGGTATGGTCAAAGGTATTCTCTGGATGCAGGATCATGAGCTgcg                                                                                                                                                >  2:128970/2‑251 (MQ=255)
                                                                                         cGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAACATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGATCTGcgcg                                                                                                                                              >  1:145265/1‑251 (MQ=255)
                                                                                                 cccATAAACTCTGCGCTCGCGGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCAttt                                                                                                                                      <  1:291898/251‑1 (MQ=255)
                                                                                                     tAAACTCTGCGCTCGCTGGGAAGAAATTTCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGAGTCATATTGATGATATGGTCAACGGTATTCTATGGCTGGTGGTTAGCGCTCTGCGCGGGCCAtgaact                                                                                                                                   >  2:70044/1‑245 (MQ=255)
                                                                                                       aagtcggcgcTGGCTGGGAAGAAATGTCCTGCCGGGCGCAAAGTGACAAACCGCGAGTGTGCCTGCTGGGTACCGGGGTAGTGCTGGCCCCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATg                                                                                                                                <  2:102503/245‑1 (MQ=255)
                                                                                                         ctctGCGCTCGCTGGGAAGAAATTGCCTGCCAGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGAAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGAGGGCCATTTAATATGGt                                                                                                                              >  1:224079/1‑251 (MQ=255)
                                                                                                             gcgcTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCg                                                                                                                          <  1:354856/251‑1 (MQ=255)
                                                                                                              cgctcgctGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGc                                                                                                                         >  1:32846/1‑251 (MQ=255)
                                                                                                                  cgctGGTAAGAAACTGCCAGCCGCGCGCATAGTGACAAAACGCGAGTGTGCCTCCGGAGGACCGGTGTAGTGCGGGCACCGGAGTGCGGTACTGTCGATAAAATGCGGCCGCCGTTTTGGCTTGGCCTGAGCGGGCAGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGAGGAGATGGGCAACGTTATTCTCGGGCTGCTTGTTAACGATCTGCGCGGGCCATTTAATATGGTTTCGCTCTa                                                                                                                     <  2:211081/251‑1 (MQ=255)
                                                                                                                  cgctGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGccc                                                                                                                       <  1:163718/249‑1 (MQ=255)
                                                                                                                     tGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTAGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGGTAACGAGCTGCGCGGGACATTTCATATGGGTTCGCCCTcccc                                                                                                                  >  2:250788/1‑251 (MQ=255)
                                                                                                                      gggAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTAcc                                                                                                                   <  1:276078/249‑1 (MQ=255)
                                                                                                                           gaaATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACg                                                                                                            >  1:4045/1‑251 (MQ=255)
                                                                                                                            aaaTTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGc                                                                                                           >  1:447535/1‑251 (MQ=255)
                                                                                                                              aTTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGCTTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCa                                                                                                          >  2:382152/1‑250 (MQ=255)
                                                                                                                                   ctgccGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAAc                                                                                                    >  1:203829/1‑251 (MQ=255)
                                                                                                                                        gggCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGGTAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAAttt                                                                                               <  1:198517/251‑1 (MQ=255)
                                                                                                                                            gcAAAGTGACAAAACGCGAGTGTGCCTGCTGGGTACCGGTGTAGTGCTGGAACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCAGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCACTGAACAATTTGccc                                                                                           >  1:189917/1‑251 (MQ=255)
                                                                                                                                            ccAAAGTGACAAAACGCGAGGTGGCCGGTCGCGTCCCGGTGTAGTGCTGGCACCGGAGGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCGGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGccc                                                                                           <  2:276165/250‑1 (MQ=255)
                                                                                                                                             cAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGccc                                                                                           >  1:67986/1‑250 (MQ=255)
                                                                                                                                             cAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCa                                                                                          <  1:173843/251‑1 (MQ=255)
                                                                                                                                              aaaGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCAt                                                                                         >  1:321689/1‑251 (MQ=255)
                                                                                                                                               aaGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATg                                                                                        >  1:217835/1‑251 (MQ=255)
                                                                                                                                                                                 cGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAttt                                                       >  2:432240/1‑251 (MQ=255)
                                                                                                                                                                                   gtgtAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCg                                                                  >  1:134206/1‑238 (MQ=255)
                                                                                                                                                                                   gtgtAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCg                                                                  <  2:134206/238‑1 (MQ=255)
                                                                                                                                                                                          gCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTc                                               <  1:397793/250‑1 (MQ=255)
                                                                                                                                                                                                          ggTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCCTGGCCTGAGCGGGCCGATTGGCTCTGGACGGCAGTATCTGGCCTGGCTTCATATTGATGATATGGTCCACGGTATTCTCTGGCTGCTGGCTAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCCTGCGCTGCATCGCCCGGCCCTTTTGCCCGTCCCCGCGCCCGCCATTCg                              >  2:207218/1‑251 (MQ=255)
                                                                                                                                                                                                              ttCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGCGCGGGCCGCTTGGCTCTGGTAGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCCATGACCAATTTGCCCATGCGCTCTGTCATGCGCTGCCTCGCCCGGCCATTTTGCGCGTCCCCGCGCCCGCCCTTCGACTg                          >  2:271194/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCGGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAAt                     <  1:350649/251‑1 (MQ=255)
                                                                                                                                                                                                                       aaaTGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTCGCTCTGGTCGGCAGTATCTGGCCTGGCTTCCTATTGATGATATGGTCCACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATgggg                 >  2:200952/1‑250 (MQ=255)
                                                                                                                                                                                                                         aTGCTGACGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTAGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGCCCGCCATTCGACTGTTAATGGGCg                >  1:364113/1‑250 (MQ=255)
                                                                                                                                                                                                                               ccgccgTTTGTTCTGGGCCTGAGCGGGCCGATTGGCTCTGGCCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCAGGCGCTGCATCGCCCGGCCATTTTGCGCTTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTc         <  2:217835/251‑1 (MQ=255)
                                                                                                                                                                                                                               ccgccgTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCATGATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTc         <  1:227920/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ccgTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGt      >  1:401569/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTTGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCcgcg                                         <  1:124603/214‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGGCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTAc    >  1:207253/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACAAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCcgcg                                         >  2:124603/1‑214 (MQ=255)
                                                                                                                                                                                                                                      ttCGTATTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTTTTAATGGGCGAATCTTCAGTACTg  >  2:161738/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tCGTTTTGGCCTGAGGGGGCCGATTGGCTCTGGTCGGAATTATCTGGCCTGGATTCATATTGATTATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTg  <  2:321689/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
CACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTG  >  NC_000913/2422021‑2422499

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