Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,375,584 T→G D25A (GAT→GCT)  epmB ← lysine 2,3‑aminomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,375,5840TG98.8% 641.6 / ‑2.0 169D25A (GAT→GCT) epmBlysine 2,3‑aminomutase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (73/94);  total (74/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 170 total aligned reads displayed.

GATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACGTACTATAGCAACG  >  NC_000913/4375447‑4375720
                                                                                                                                         |                                                                                                                                        
   cGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTc                                                                                                                                     >  2:953076/1‑140 (MQ=255)
   cGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTc                                                                                                                                     >  1:1868075/1‑140 (MQ=255)
     tCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTcac                                                                                                                                   >  1:644683/1‑140 (MQ=255)
     tCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTcac                                                                                                                                   >  1:1958523/1‑140 (MQ=255)
     tCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTcac                                                                                                                                   >  1:1516861/1‑140 (MQ=255)
       cGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCacaa                                                                                                                                 <  1:304823/140‑1 (MQ=255)
            tGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCg                                                                                                                            >  1:1415782/1‑140 (MQ=255)
                 ttCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAg                                                                                                                       <  1:2453615/140‑1 (MQ=255)
                 ttCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAg                                                                                                                       <  2:1712061/140‑1 (MQ=255)
                   ctcCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGtt                                                                                                                     >  1:450080/1‑140 (MQ=255)
                   ctcCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGtt                                                                                                                     <  2:2511509/140‑1 (MQ=255)
                      cATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCg                                                                                                                  <  1:2524629/140‑1 (MQ=255)
                        tGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGtt                                                                                                                >  2:1459931/1‑140 (MQ=255)
                        tGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGtt                                                                                                                >  1:1212375/1‑140 (MQ=255)
                        tGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGtt                                                                                                                >  1:2220086/1‑140 (MQ=255)
                        tGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGtt                                                                                                                >  1:80952/1‑140 (MQ=255)
                        tGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGtt                                                                                                                >  1:1117886/1‑140 (MQ=255)
                          cgatcgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTaa                                                                                                              <  1:943717/140‑1 (MQ=255)
                           gatcgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAAc                                                                                                             <  1:991761/140‑1 (MQ=255)
                            atcgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAg                                                                                                                       >  2:684894/1‑129 (MQ=255)
                            atcgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAAcc                                                                                                            >  1:2400941/1‑140 (MQ=255)
                             tcgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCa                                                                                                           >  2:988951/1‑140 (MQ=255)
                              cgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCaa                                                                                                          >  1:1753366/1‑140 (MQ=255)
                               gatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         <  1:591795/140‑1 (MQ=255)
                               gatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         <  1:1726558/140‑1 (MQ=255)
                               gatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         <  2:1868075/140‑1 (MQ=255)
                                  aaaTGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGAAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAtctt                                                                                                      >  2:2494157/1‑140 (MQ=255)
                                   aaTGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAtcttc                                                                                                     <  2:2348554/140‑1 (MQ=255)
                                    aTGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAtcttct                                                                                                    <  1:1459931/140‑1 (MQ=255)
                                       aGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTg                                                                                                 >  1:2162879/1‑140 (MQ=255)
                                        gCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTgg                                                                                                <  2:1106992/140‑1 (MQ=255)
                                         cGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGa                                                                                               >  2:2294506/1‑140 (MQ=255)
                                         cGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         >  2:2113034/1‑130 (MQ=255)
                                         cGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         <  1:2113034/130‑1 (MQ=255)
                                         cGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGa                                                                                               >  2:1698899/1‑140 (MQ=255)
                                         cGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGa                                                                                               >  2:1063635/1‑140 (MQ=255)
                                         cGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGa                                                                                               <  2:2143508/140‑1 (MQ=255)
                                          ggggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTATCCAATCTTCTCTGGAt                                                                                              <  2:2001249/140‑1 (MQ=255)
                                           gggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAAc                                                                                                             <  2:1305214/124‑1 (MQ=255)
                                           gggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAAc                                                                                                             >  1:1305214/1‑124 (MQ=255)
                                           gggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATg                                                                                             >  1:1065944/1‑140 (MQ=255)
                                           gggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATg                                                                                             >  2:1292339/1‑140 (MQ=255)
                                           gggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATg                                                                                             >  2:2271076/1‑140 (MQ=255)
                                            ggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATgg                                                                                            <  1:953076/140‑1 (MQ=255)
                                            ggCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATgg                                                                                            <  1:1063635/140‑1 (MQ=255)
                                             gCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATggg                                                                                           <  1:1251029/140‑1 (MQ=255)
                                             gCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATggg                                                                                           <  1:429613/140‑1 (MQ=255)
                                             gCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATggg                                                                                           <  1:1313455/140‑1 (MQ=255)
                                              cacaCGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         >  1:1999911/1‑125 (MQ=255)
                                              cacaCGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAAt                                                                                                         <  2:1999911/125‑1 (MQ=255)
                                                  cGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTAtt                                                                                      <  2:683206/140‑1 (MQ=255)
                                                   gCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTAttt                                                                                     >  1:520614/1‑140 (MQ=255)
                                                    cAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTa                                                                                    <  1:628187/140‑1 (MQ=255)
                                                    cAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTa                                                                                    <  2:1547802/140‑1 (MQ=255)
                                                    cAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTa                                                                                    <  2:1582322/140‑1 (MQ=255)
                                                    cAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTa                                                                                    <  2:1117886/140‑1 (MQ=255)
                                                      gggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTAtt                                                                                      >  2:555359/1‑136 (MQ=255)
                                                      gggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAgg                                                                                  >  2:2436446/1‑140 (MQ=255)
                                                      gggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAg                                                                                   <  2:1190389/139‑1 (MQ=255)
                                                      gggCAAAAAGCCTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAgg                                                                                  >  1:623183/1‑140 (MQ=255)
                                                       ggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAggg                                                                                 >  2:1744844/1‑140 (MQ=255)
                                                       ggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCATTAACCAATCTTCTCTGGATGGGGTATTTAggg                                                                                 <  2:1205079/140‑1 (MQ=255)
                                                           aaaaGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTAc                                                                             >  2:226816/1‑140 (MQ=255)
                                                            aaaGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACa                                                                            <  1:2298004/140‑1 (MQ=255)
                                                              aGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAat                                                                          <  1:2271076/140‑1 (MQ=255)
                                                                   tttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATgcg                                                                     <  1:148370/140‑1 (MQ=255)
                                                                         gcTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCAttt                                                               <  2:974192/140‑1 (MQ=255)
                                                                         gcTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCAttt                                                               <  2:520614/140‑1 (MQ=255)
                                                                          cTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCAtttt                                                              <  1:1211870/140‑1 (MQ=255)
                                                                           tGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTtg                                                             <  2:1415782/140‑1 (MQ=255)
                                                                                ccGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCt                                                        <  1:1292339/140‑1 (MQ=255)
                                                                                      aaCAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTa                                                  <  1:2271379/140‑1 (MQ=255)
                                                                                           tttttCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAAt                                             <  2:450080/140‑1 (MQ=255)
                                                                                           tttttCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAAt                                             <  2:2162879/140‑1 (MQ=255)
                                                                                             tttCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGCTGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTa                                           >  2:1651313/1‑140 (MQ=255)
                                                                                                      cgTCTATATTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTc                                  >  1:373014/1‑140 (MQ=255)
                                                                                                          tataTTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCagag                              <  1:1194834/140‑1 (MQ=255)
                                                                                                          tataTTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCagag                              <  2:81986/140‑1 (MQ=255)
                                                                                                          tataTTCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCagag                              <  1:2472103/140‑1 (MQ=255)
                                                                                                               tCAAAAGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCt                         <  2:1956198/140‑1 (MQ=255)
                                                                                                                 aaaaGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTATGGTTACAATATGCGCCGTTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTa                       >  2:931140/1‑140 (MQ=255)
                                                                                                                 aaaaGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCagag                              <  1:293014/133‑1 (MQ=255)
                                                                                                                   aaGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATg                     >  2:362631/1‑140 (MQ=255)
                                                                                                                    aGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATgg                    <  2:1212375/140‑1 (MQ=255)
                                                                                                                    aGACGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATgg                    <  2:623183/140‑1 (MQ=255)
                                                                                                                      aCGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCa                  <  1:86277/140‑1 (MQ=255)
                                                                                                                      aCGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCa                  >  2:1308230/1‑140 (MQ=255)
                                                                                                                      aCGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCa                  <  1:114787/140‑1 (MQ=255)
                                                                                                                      aCGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCa                  >  2:86277/1‑140 (MQ=255)
                                                                                                                      aCGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCa                  <  2:373014/140‑1 (MQ=255)
                                                                                                                      aCGCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCa                  <  1:818741/140‑1 (MQ=255)
                                                                                                                        gCAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAAc                >  2:811843/1‑140 (MQ=255)
                                                                                                                         cAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACg               <  1:1153921/140‑1 (MQ=255)
                                                                                                                         cAGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACg               <  1:167261/140‑1 (MQ=255)
                                                                                                                          aGAAGTTCATCAGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACGt              <  2:2220086/140‑1 (MQ=255)
                                                                                                                                 catcaGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATgg                    >  2:413301/1‑127 (MQ=255)
                                                                                                                                 catcaGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACGTACTATAg       >  2:1818142/1‑140 (MQ=255)
                                                                                                                                 catcaGGAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACGTACTATAg       >  2:830214/1‑140 (MQ=255)
                                                                                                                                      ggAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATgg                    >  2:2088517/1‑122 (MQ=255)
                                                                                                                                      ggAGCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATgg                    <  1:2088517/122‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGTTAACCAATCTTCTCTGGATGGGGTATTTAGGGTTACAATATGCGCCATTTTGTGGCTTAGCTACCAATTAACAAATTTCAGAGGGCCTTATGGCAACGTACTATAGCAACG  >  NC_000913/4375447‑4375720

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