Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,715,919 T→G 5.5% intergenic (‑30/+29) 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,9190TG 5.5% 1092.2 / 34.3 338intergenic (‑30/+29)pdxY/tyrSpyridoxamine kinase/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (145/173);  new base G (17/3);  total (162/176)
Fisher's exact test for biased strand distribution p-value = 7.78e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

TGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATG  >  NC_000913/1715771‑1716068
                                                                                                                                                    |                                                                                                                                                     
tGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAtt                                                                                                                                                      <  1:4564065/150‑1 (MQ=255)
   aCGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtct                                                                                                                                                   <  1:4704376/150‑1 (MQ=255)
   aCGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtct                                                                                                                                                   <  2:3672116/150‑1 (MQ=255)
     ggTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCa                                                                                                                                                 <  1:4931644/150‑1 (MQ=255)
     ggTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCa                                                                                                                                                 <  2:4277367/150‑1 (MQ=255)
       tgtTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACt                                                                                                                                               >  2:2047318/1‑150 (MQ=255)
        gtTCAGCGGCCAGCCGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTc                                                                                                                                              <  1:2911004/150‑1 (MQ=255)
        gtTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTc                                                                                                                                              <  1:119239/150‑1 (MQ=255)
         ttCAGCGGCCAGACGTTCGCGCCCAGGCGGCGGATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTGTCACTcc                                                                                                                                             <  2:4160115/150‑1 (MQ=255)
         ttCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcc                                                                                                                                             <  1:2324/150‑1 (MQ=255)
          tCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTccc                                                                                                                                            <  2:2442050/150‑1 (MQ=255)
           cAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                           <  1:2035931/150‑1 (MQ=255)
            aGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTccccc                                                                                                                                          <  1:635772/150‑1 (MQ=255)
            aGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTccccc                                                                                                                                          >  2:3720379/1‑150 (MQ=255)
             gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGACACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATCGCGAGAATATTTTTCATCATGTTCCCTGGATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCt                                                                                                                                         >  2:3099098/1‑150 (MQ=255)
             gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTCCCCCt                                                                                                                                         <  2:462928/150‑1 (MQ=255)
                 ccAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGTTCTCCCTCCCCCTTTcc                                                                                                                                     >  1:4361830/1‑150 (MQ=255)
                  cAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCa                                                                                                                                    <  2:4428009/150‑1 (MQ=255)
                   aGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTGACTCCCCCGTTCCAc                                                                                                                                   >  1:2427274/1‑150 (MQ=255)
                   aGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCAc                                                                                                                                   >  1:4277714/1‑150 (MQ=255)
                   aGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCAc                                                                                                                                   >  2:4743121/1‑150 (MQ=255)
                        ttCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATACACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAAt                                                                                                                              <  2:1131478/150‑1 (MQ=255)
                        ttCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAAt                                                                                                                              <  2:1150253/150‑1 (MQ=255)
                        ttCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGAGATTTCTCACTCCCCCTTTCCACTTAAt                                                                                                                              <  2:1431632/150‑1 (MQ=255)
                           gcgcCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCa                                                                                                                           <  2:4805060/150‑1 (MQ=255)
                            cgcCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCat                                                                                                                          >  1:599172/1‑150 (MQ=255)
                             gcCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCatt                                                                                                                         >  2:5189139/1‑150 (MQ=255)
                                cAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGGTTTCTCTCTCCCCCTTTCCACTTAATGCattatt                                                                                                                      >  1:5200686/1‑150 (MQ=255)
                                 aGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTAttt                                                                                                                     <  1:953109/150‑1 (MQ=255)
                                    cggcgCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCa                                                                                                                  <  2:2683994/150‑1 (MQ=255)
                                            cGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCa                                                                                                          <  1:3142443/150‑1 (MQ=255)
                                              gAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGGTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGa                                                                                                        >  1:1049610/1‑150 (MQ=255)
                                                   tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTa                                                                                                   >  2:481678/1‑150 (MQ=255)
                                                   tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTa                                                                                                   >  2:587942/1‑150 (MQ=255)
                                                       ccGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTc                                                                                               <  2:2316554/150‑1 (MQ=255)
                                                       ccGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACGCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTc                                                                                               >  1:3330324/1‑150 (MQ=255)
                                                           aCTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtttt                                                                                           >  1:362465/1‑150 (MQ=255)
                                                           aCTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtttt                                                                                           <  1:2817584/150‑1 (MQ=255)
                                                            cTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttttt                                                                                          <  2:1280821/150‑1 (MQ=255)
                                                                tACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                                      <  1:4131204/150‑1 (MQ=255)
                                                                tACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                                      <  2:316789/150‑1 (MQ=255)
                                                                 aCCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcgc                                                                                     >  2:1667045/1‑150 (MQ=255)
                                                                   ccGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGa                                                                                   <  2:4079283/150‑1 (MQ=255)
                                                                       aTGACCATAAACACCGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                                               <  2:3461292/150‑1 (MQ=255)
                                                                       aTGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                                               <  1:2790716/150‑1 (MQ=255)
                                                                       aTGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATGTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                                               <  2:3908880/150‑1 (MQ=255)
                                                                        tGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAg                                                                              <  2:5104283/150‑1 (MQ=255)
                                                                         gACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGt                                                                             >  2:1585303/1‑150 (MQ=255)
                                                                            cATAAACACCGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                          <  2:1583895/150‑1 (MQ=255)
                                                                            cATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                          <  1:2715475/150‑1 (MQ=255)
                                                                            cATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                          <  2:1868933/150‑1 (MQ=255)
                                                                             aTAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGTTTTGTCAGTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACAGCGACGCAGTAAgg                                                                         >  1:3268698/1‑150 (MQ=255)
                                                                              tAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                                        <  1:747230/150‑1 (MQ=255)
                                                                              tAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                                        >  1:2467335/1‑150 (MQ=255)
                                                                                 acaacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                                     >  1:3396034/1‑150 (MQ=255)
                                                                                 acaacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                                     <  1:5358615/150‑1 (MQ=255)
                                                                                 acaacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                                     >  2:4895300/1‑150 (MQ=255)
                                                                                           aCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                                           >  2:4432509/1‑150 (MQ=255)
                                                                                           aCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATGGCTCGCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACTGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                                           >  1:1183127/1‑150 (MQ=255)
                                                                                               gATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGGTTTCTCTCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGAc                                                       >  1:1609399/1‑150 (MQ=255)
                                                                                               gATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCAGTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGAc                                                       >  1:707645/1‑150 (MQ=255)
                                                                                                   gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                                   >  2:822864/1‑150 (MQ=255)
                                                                                                     gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCCCTCCCCCTTGCCACTGAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTCAAACGACCAAACAGACGAtctt                                                 >  1:2874912/1‑150 (MQ=255)
                                                                                                      ataaTATTTTTCATCATGTTCCATTTATAAAAACCAGGGGGAGTGTTTTCTCACTCCCCATTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGTTAAAACGACCAAACAGACGAtcttc                                                <  1:3198181/148‑1 (MQ=255)
                                                                                                      agaATATTTTTCATCATGTTCCCTGTATAAAAACGAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttc                                                <  2:2843391/150‑1 (MQ=255)
                                                                                                      agaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttc                                                <  2:3266094/150‑1 (MQ=255)
                                                                                                        aaTATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttctt                                              >  2:913765/1‑150 (MQ=255)
                                                                                                         atatTTTTCATCCTGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCATTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                                             <  2:5301455/150‑1 (MQ=255)
                                                                                                          tatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCTTTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttct                                            <  1:1762145/150‑1 (MQ=255)
                                                                                                              tttCATCATGTTCCCTTTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGAAAAACGACCAAACAGACGATCTTCTTCTTTaa                                        <  1:1441482/150‑1 (MQ=255)
                                                                                                               ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                                       >  2:352580/1‑150 (MQ=255)
                                                                                                               ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                                       >  2:400217/1‑150 (MQ=255)
                                                                                                               ttCATCATGTTCCCTGAATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTAAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACAGCGACGCAGTAAGGTAAAACGACCAAACAGACTATCTTCTTCTTTaaa                                       >  1:970725/1‑150 (MQ=255)
                                                                                                                  atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                                    <  2:2739006/150‑1 (MQ=255)
                                                                                                                  atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                                    <  1:4020444/150‑1 (MQ=255)
                                                                                                                  atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                                    <  2:4183877/150‑1 (MQ=255)
                                                                                                                   tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag                                   >  2:4567970/1‑150 (MQ=255)
                                                                                                                    catGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGt                                  <  1:1674419/150‑1 (MQ=255)
                                                                                                                    catGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGt                                  >  1:5128060/1‑150 (MQ=255)
                                                                                                                      tGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTAt                                >  1:3914593/1‑150 (MQ=255)
                                                                                                                      tGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCCCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTAt                                >  1:1160232/1‑150 (MQ=255)
                                                                                                                         tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa                             <  1:802824/150‑1 (MQ=255)
                                                                                                                         tCCCTGTAAAAAAACCAGGGGCAGTGATTACTCACTCCCCCTTTCCACTAAATGAAAAATTTCCAGCAAATCACACACTAATACTTTTTACCCCGACCCACTAAGCTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa                             <  1:4868415/150‑1 (MQ=255)
                                                                                                                            cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa                          >  2:2171522/1‑150 (MQ=255)
                                                                                                                            cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa                          >  2:3005334/1‑150 (MQ=255)
                                                                                                                               tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTACTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCg                       <  2:5256331/150‑1 (MQ=255)
                                                                                                                               tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCg                       <  1:1281003/150‑1 (MQ=255)
                                                                                                                                 taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGa                     >  2:686509/1‑150 (MQ=255)
                                                                                                                                   aaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACt                   <  2:1482267/150‑1 (MQ=255)
                                                                                                                                     aaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGt                 <  1:5189587/150‑1 (MQ=255)
                                                                                                                                         aGGGGGTGTGATTTCTCACTCCCCCTTTCCTCTTAATGCATTATTTCCCGCAAATCAGACCGTATTTCTTTTTACCGCGACGCAGTATGGTAAAACGACCACACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt             <  1:3099098/150‑1 (MQ=255)
                                                                                                                                         aGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt             >  2:3387740/1‑150 (MQ=255)
                                                                                                                                         aGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt             >  2:4055037/1‑150 (MQ=255)
                                                                                                                                         aGGGGGAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttttt             >  2:174108/1‑150 (MQ=255)
                                                                                                                                           ggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa           >  2:788731/1‑150 (MQ=255)
                                                                                                                                           ggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTAAGGATCGGACTGTTTTTCa           >  1:2668832/1‑150 (MQ=255)
                                                                                                                                           ggggAGTGAGAAATCACTCCCCCTTTCCAATTAATGCATTATTTCCAGCAAAACAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCa           >  2:2559911/1‑150 (MQ=255)
                                                                                                                                             ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCAcc         <  1:2920016/150‑1 (MQ=255)
                                                                                                                                              gAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCg        <  1:4997609/150‑1 (MQ=255)
                                                                                                                                                   aTTTCTCACTCCCCCTTTCCAATTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCTGGATCGGACTGTTTTTCACCGTTAAt   <  1:219771/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATG  >  NC_000913/1715771‑1716068

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: