Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,715,925 A→G 5.4% intergenic (‑36/+23) pdxY ← / ← tyrS pyridoxamine kinase/tyrosyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,715,9250AG 5.4% 514.0 / 4.2 164intergenic (‑36/+23)pdxY/tyrSpyridoxamine kinase/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (75/75);  new base G (7/2);  total (87/77)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

GGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATG  >  NC_000913/1715776‑1716074
                                                                                                                                                     |                                                                                                                                                     
ggTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCa                                                                                                                                                       <  1:1837408/150‑1 (MQ=255)
ggTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCa                                                                                                                                                       <  1:3612633/150‑1 (MQ=255)
      cAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGTAAACTCTGCCGCACTGTTACCCGCTTGACCATAAACAACGTGATACTGTATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCTGGGGGAGTGAGAAATCACTcccc                                                                                                                                                 <  1:1050602/150‑1 (MQ=255)
      cAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                                 <  2:4122381/150‑1 (MQ=255)
        gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGACTATTTTTCATCATGTTCCCTGGATAAAAACCAGGGGGAGTGATGTCTCACTcccccg                                                                                                                                               >  2:1118517/1‑149 (MQ=255)
        gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCt                                                                                                                                               >  2:538942/1‑150 (MQ=255)
        gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCt                                                                                                                                               >  2:3424850/1‑150 (MQ=255)
        gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCGCTCCCCCt                                                                                                                                               <  1:3546945/150‑1 (MQ=255)
         cGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGGTTTCTCACTCCCCCtt                                                                                                                                              >  1:3392815/1‑150 (MQ=255)
         cGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTCCCCCtg                                                                                                                                              <  1:3461261/150‑2 (MQ=255)
         cGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACACTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGGGTTCTCACGCCCCCtt                                                                                                                                              >  1:3392782/1‑150 (MQ=255)
             cAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAATCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCa                                                                                                                                          >  2:206703/1‑150 (MQ=255)
                  gTTCGCTACCAGGCGGCGCATCGGATACTCTGCAGCTATTTTTCCCTCATGAACATAAACAACGTGAGTCTTTATAGCTATAATATTTTTAATCATTTTCACTTTATAAAAACCTGGTGGTGTGATTTCTCACTCCCCCTTTCCACTTac                                                                                                                                     <  2:1557759/150‑2 (MQ=255)
                  gTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCCCTCCCCCTTTCCACTTaa                                                                                                                                     >  1:3068754/1‑150 (MQ=255)
                  gTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTaa                                                                                                                                     <  2:1557607/150‑1 (MQ=255)
                  gTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTaa                                                                                                                                     <  2:1993795/150‑1 (MQ=255)
                  gTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTaa                                                                                                                                     <  2:2474829/150‑1 (MQ=255)
                    tCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATg                                                                                                                                   <  1:2088286/150‑1 (MQ=255)
                     cgcgCCCAGGCGGCGCATCGGAAACTCGGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGGATACAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTGACTGc                                                                                                                                  <  2:868879/150‑1 (MQ=255)
                       cgcCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCat                                                                                                                                <  2:1475295/150‑1 (MQ=255)
                       cgcCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCat                                                                                                                                <  2:1475287/150‑1 (MQ=255)
                         cccAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCatta                                                                                                                              <  1:850330/150‑1 (MQ=255)
                              gcggcgCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACACCGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                                         <  1:1421945/150‑1 (MQ=255)
                              gcggcgCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                                         <  2:3759314/150‑1 (MQ=255)
                               cggcgCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCa                                                                                                                        <  1:1392206/150‑1 (MQ=255)
                                        ggAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGGTGTCTCGCTCACCCTTTCCACTTAATGCATTATTTCCAGCAAATCAg                                                                                                               >  1:2158951/1‑150 (MQ=255)
                                        ggAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCCCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAg                                                                                                               >  1:2717598/1‑150 (MQ=255)
                                        ggAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAg                                                                                                               >  2:1620999/1‑150 (MQ=255)
                                        ggAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAg                                                                                                               >  1:2710941/1‑150 (MQ=255)
                                         gAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGa                                                                                                              <  1:901519/150‑1 (MQ=255)
                                         gAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGa                                                                                                              <  1:3808245/150‑1 (MQ=255)
                                         gAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGa                                                                                                              <  1:3536070/150‑1 (MQ=255)
                                              tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTa                                                                                                         >  1:1949650/1‑150 (MQ=255)
                                              tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTa                                                                                                         >  2:1033106/1‑150 (MQ=255)
                                                 gCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAAtt                                                                                                      <  2:2260082/150‑1 (MQ=255)
                                                 gCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCCCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAAtt                                                                                                      >  2:3676016/1‑150 (MQ=255)
                                                    gCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtt                                                                                                   >  1:2855310/1‑150 (MQ=255)
                                                    gCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtt                                                                                                   <  2:4102540/150‑1 (MQ=255)
                                                     cACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttt                                                                                                  >  2:1429815/1‑150 (MQ=255)
                                                     cACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttt                                                                                                  <  1:1592838/150‑1 (MQ=255)
                                                      aCTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtttt                                                                                                 <  1:1865935/150‑1 (MQ=255)
                                                        tGTTACCCGCATGACTATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTa                                                                                               >  1:1834597/1‑150 (MQ=255)
                                                        tGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCTCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTa                                                                                               >  1:1834537/1‑150 (MQ=255)
                                                         gtgaCCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGGTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                                              <  2:2746819/147‑1 (MQ=255)
                                                         gTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                                              <  2:275540/150‑1 (MQ=255)
                                                           tACCCGCATGACCATAAACAACGTGAGACTGGATATCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGTGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                                            <  1:2907063/150‑1 (MQ=255)
                                                            aCCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcgc                                                                                           >  1:818574/1‑150 (MQ=255)
                                                                gCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACg                                                                                       >  2:2706297/1‑150 (MQ=255)
                                                                 cATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGc                                                                                      <  2:1647588/150‑1 (MQ=255)
                                                                       cATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCTCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                                >  1:494168/1‑150 (MQ=255)
                                                                        aTAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                                               >  1:1546919/1‑150 (MQ=255)
                                                                        aTAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                                               >  1:4125999/1‑150 (MQ=255)
                                                                        aTAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                                               >  1:2233723/1‑150 (MQ=255)
                                                                         tAAACAACCTGAGACTGGAGAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                                              <  2:3138419/150‑1 (MQ=255)
                                                                             caacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaaa                                                                          <  1:3580445/150‑1 (MQ=255)
                                                                             caacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaaa                                                                          <  1:2521736/150‑1 (MQ=255)
                                                                              aacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAAc                                                                         >  2:353325/1‑150 (MQ=255)
                                                                                 gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAATCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCTGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAATCGAc                                                                      >  2:322937/1‑150 (MQ=255)
                                                                                  tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                                                     >  2:2442262/1‑150 (MQ=255)
                                                                                     gaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaaa                                                                  >  2:4259233/1‑150 (MQ=255)
                                                                                          gATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGAc                                                             <  2:1582735/150‑1 (MQ=255)
                                                                                              gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTACTTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                                         >  2:3342135/1‑150 (MQ=255)
                                                                                                gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtctt                                                       >  1:617826/1‑150 (MQ=255)
                                                                                                 agaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttc                                                      <  2:2743956/150‑1 (MQ=255)
                                                                                                   aaTATTTTTCATCATGTTCCCTGTATAAAAACCTGGGGGAGTTATTTATCACTCCCCCTTTCCACTTAATGAATTATTTCCAGAAAATCAGACATTAATTCTTTTTACCGCGACGCAGTAAGTTAAAACGACCAAACAGACGAtcttctt                                                    <  1:1702032/150‑1 (MQ=255)
                                                                                                    atatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGTAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                                                   <  1:4066678/150‑1 (MQ=255)
                                                                                                      atTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                                                 <  2:451084/150‑1 (MQ=255)
                                                                                                      atTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                                                 <  2:451030/150‑1 (MQ=255)
                                                                                                       tttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGGTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCttt                                                >  1:1821154/1‑150 (MQ=255)
                                                                                                        ttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTa                                               >  2:1881012/1‑150 (MQ=255)
                                                                                                         tttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaa                                              >  1:3208144/1‑150 (MQ=255)
                                                                                                          ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTTATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                                             <  1:1189771/150‑1 (MQ=255)
                                                                                                          ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                                             >  2:2422638/1‑150 (MQ=255)
                                                                                                           tcatcaTGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaag                                            <  2:3393773/150‑1 (MQ=255)
                                                                                                                    tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTCATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa                                   <  1:2101111/150‑1 (MQ=255)
                                                                                                                        tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAACTTCTTCTTTAAAGAAGTATTCAGGAt                               >  1:2436461/1‑150 (MQ=255)
                                                                                                                            taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACATTAATTCTTTTTACCGCGACGCAGTAAGTTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGa                           <  1:4006474/150‑1 (MQ=255)
                                                                                                                              aaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACt                         <  2:1861825/150‑1 (MQ=255)
                                                                                                                                  ccAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttt                     <  1:2753037/150‑1 (MQ=255)
                                                                                                                                   cAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATACTTTTTAACGCGAAGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTAAGGATCGGACTGtttt                    >  1:1908758/1‑150 (MQ=255)
                                                                                                                                    aGGGGGAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt                   >  2:263972/1‑150 (MQ=255)
                                                                                                                                    aGGGGGAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt                   >  2:2431451/1‑150 (MQ=255)
                                                                                                                                     gggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTc                  >  2:1895438/1‑150 (MQ=255)
                                                                                                                                      ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa                 >  2:3061287/1‑150 (MQ=255)
                                                                                                                                      ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa                 <  2:2465538/150‑1 (MQ=255)
                                                                                                                                      ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa                 >  2:2105820/1‑150 (MQ=255)
                                                                                                                                      ggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa                 >  2:2147010/1‑150 (MQ=255)
                                                                                                                                        ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAATGAAGTATTCAGGATCGGACTGtttttatcc               >  2:2929457/1‑146 (MQ=255)
                                                                                                                                             gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTaa          <  1:840945/150‑1 (MQ=255)
                                                                                                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAt         <  1:2070597/150‑1 (MQ=255)
                                                                                                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAt         <  2:238844/150‑1 (MQ=255)
                                                                                                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGAATTCAGGATCGGACTGTTTTTCACCGTTAAt         >  2:3506835/1‑150 (MQ=255)
                                                                                                                                                ttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATgg       >  2:1949200/1‑150 (MQ=255)
                                                                                                                                                ttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATgg       <  1:1934897/150‑1 (MQ=255)
                                                                                                                                                ttCTCACTCCCCCTTTCCACTTAATGCATGATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATgg       <  1:3886644/150‑1 (MQ=255)
                                                                                                                                                  ctcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAATAAGTATACAGGATCGGACTGTTTTTCACCGTTAATGGTg     >  1:2680435/1‑150 (MQ=255)
                                                                                                                                                  ctcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTg     <  1:1049253/150‑1 (MQ=255)
                                                                                                                                                  cacACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTg     >  1:4159935/3‑150 (MQ=255)
                                                                                                                                                   tcACTCCCCCTTTCCACTTAATGCATTAGTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGa    <  1:4167550/150‑1 (MQ=255)
                                                                                                                                                    cACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGAt   >  1:3577972/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATG  >  NC_000913/1715776‑1716074

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: