Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,907,726 T→G 5.6% Q214P (CAG→CCG)  pstB ← phosphate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,907,7260TG 5.6% 316.7 / 5.2 101Q214P (CAG→CCG) pstBphosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (53/42);  new base G (3/3);  total (56/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

AACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCG  >  NC_000913/3907579‑3907873
                                                                                                                                                   |                                                                                                                                                   
aaCGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCctg                                                                                                                                                    >  2:1505296/1‑149 (MQ=255)
aaCGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGGGTTGCTGAACTCAATCAATTCGCCCAGGGACATAAACGCCGTGTGGTCGGAACAACGCGCAGCctg                                                                                                                                                    <  2:259204/149‑1 (MQ=255)
 aCGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCctgc                                                                                                                                                   <  2:668009/149‑1 (MQ=255)
     aCTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTAGCCCAGGTACATAAACGCCGTGTGGTCGGAACAATGCGCAGCCTGCTGCa                                                                                                                                               >  1:1411456/1‑149 (MQ=255)
         cTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCAtgtt                                                                                                                                           >  2:411099/1‑149 (MQ=255)
           gAATCAACCGTAACCTCCGTTGATGTAGTCTTCAGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGAACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCAtgttgt                                                                                                                                         <  2:1086459/149‑1 (MQ=255)
            aaTCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTtctg                                                                                                                                        >  2:153201/1‑149 (MQ=255)
            aaTCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACATGCAGCAGGCTGCGcgttgtt                                                                                                                                        <  1:597072/149‑2 (MQ=255)
             agcacccGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCACTCAATTCGCCCAGGTCCCTACTCGCCGCGGGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTgg                                                                                                                                       <  2:1132891/144‑1 (MQ=255)
             aTCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTgg                                                                                                                                       <  1:838291/149‑1 (MQ=255)
             aTCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTgg                                                                                                                                       <  1:367676/149‑1 (MQ=255)
              tCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGCTCGTCCGTGTTGCTGCACTCAATCAATTCGCCCAGGTACTTAAACGCCGTGTGGTCGGAACACCATGCAGCAGGCTGCGcgttgttcc                                                                                                                                      >  1:1013083/1‑146 (MQ=255)
              tCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACATGCAGCAGGCTGCGcgttgttcc                                                                                                                                      <  2:18913/149‑4 (MQ=255)
               cAACCGTAACGACCTGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCTTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGt                                                                                                                                     <  1:965121/149‑1 (MQ=255)
               cAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGt                                                                                                                                     <  1:970368/149‑1 (MQ=255)
                  ccGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGAc                                                                                                                                  <  1:411099/149‑1 (MQ=255)
                  ccGTAACGACCGGTGATGTAGTCTTCTGCTTGTTTCTTCGCGGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGAc                                                                                                                                  <  2:1510185/149‑1 (MQ=255)
                    gTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGa                                                                                                                                <  1:2572/149‑1 (MQ=255)
                     tAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGAt                                                                                                                               <  2:1252183/149‑1 (MQ=255)
                     tAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACGCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGAt                                                                                                                               <  1:600454/149‑1 (MQ=255)
                     tAACGACCGGTGATGTAGTCTTCTGTGTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGAt                                                                                                                               <  2:829562/149‑1 (MQ=255)
                       aCGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATca                                                                                                                             >  2:1310424/1‑149 (MQ=255)
                          aCCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGCACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGAGGTGGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATcacca                                                                                                                          >  1:262675/1‑149 (MQ=255)
                            cGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACg                                                                                                                        >  2:893314/1‑149 (MQ=255)
                            cGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACg                                                                                                                        >  1:1016683/1‑149 (MQ=255)
                              gTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGgt                                                                                                                      <  2:695437/149‑1 (MQ=255)
                              gTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACACGCCGCCTGCTGCATGTTGTGGGTGACGATCAACACGgt                                                                                                                      >  2:220535/1‑149 (MQ=255)
                              gTGATGTAGTCTTCTGTATGTTTCTTCTCTGGCTTGGTGAACAGATCGTCCGTGTTGCGGAACTCAATCACTTCGCCCAGGTACATAAACGCCGTGTGGGCGGAACAACGCGCAGCCTGCTGCATGTTTTGGGTGACGATCACCTCGgt                                                                                                                      <  1:445870/149‑1 (MQ=255)
                               tGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGgtg                                                                                                                     <  2:1260324/149‑1 (MQ=255)
                                  tGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTaa                                                                                                                  >  1:457820/1‑149 (MQ=255)
                                    tAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAGTTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATc                                                                                                                <  1:709967/149‑1 (MQ=255)
                                     aGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATcc                                                                                                               <  1:513874/149‑1 (MQ=255)
                                      gTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCt                                                                                                              >  2:1267727/1‑149 (MQ=255)
                                       tcttctGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTg                                                                                                             <  1:1076029/149‑1 (MQ=255)
                                         ttctGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCt                                                                                                           >  1:1510952/1‑149 (MQ=255)
                                          tctGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCtt                                                                                                          <  2:293754/149‑1 (MQ=255)
                                            tgtttgttTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCa                                                                                                        <  1:1510465/149‑1 (MQ=255)
                                             gtttgtttCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAg                                                                                                       <  1:942994/149‑1 (MQ=255)
                                                tgtttCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTc                                                                                                    <  1:586338/149‑1 (MQ=255)
                                                ggtttCTTCGCTGGCTTGGTGAACAGATCGTCCGGGTTGCTCACCTCAATCAATGCGCGCAGGTCCATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCTTGGTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTc                                                                                                    <  1:78021/148‑1 (MQ=255)
                                                 gtttCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCg                                                                                                   <  1:1201012/149‑1 (MQ=255)
                                                    tcttcGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGt                                                                                                 >  1:378909/1‑148 (MQ=255)
                                                    tcttcGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTg                                                                                                >  2:604747/1‑149 (MQ=255)
                                                      ttcGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGAt                                                                                              >  2:373648/1‑149 (MQ=255)
                                                      ttcGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGAt                                                                                              <  2:558059/149‑1 (MQ=255)
                                                      ttcGCTGGCTTGGTGAACAGATCGTCCGTGTGGCTGAACTCAATCACTTCGCGCAGGTTCATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGGGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGAt                                                                                              <  1:199605/149‑1 (MQ=255)
                                                       tcGCTGGCGTGGTGACCAGAGCGTCCGTGTTGCTGATCTCAATCATAGCGCCCATGTACATAAACGCCGAGGGGTCGGAACAACGCGACGCCTGCTGCATGGTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATc                                                                                             <  1:259095/149‑1 (MQ=255)
                                                        cGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCa                                                                                            <  1:755020/149‑1 (MQ=255)
                                                              cTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCtt                                                                                      <  2:752160/149‑1 (MQ=255)
                                                                 ggTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCaa                                                                                   <  1:407006/149‑1 (MQ=255)
                                                                  gTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTAGGTGATCAGCTCATCAAt                                                                                  >  2:68928/1‑149 (MQ=255)
                                                                     aaCAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACg                                                                               >  2:287637/1‑149 (MQ=255)
                                                                        aGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCCGTTCGGTGATCAGCTCTTCAATACGccc                                                                            >  1:577072/1‑149 (MQ=255)
                                                                          aTCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCgg                                                                          >  2:437853/1‑149 (MQ=255)
                                                                          aTCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCgg                                                                          >  1:1103782/1‑149 (MQ=255)
                                                                                  tgtTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATa                                                                  >  1:935942/1‑149 (MQ=255)
                                                                                  tgtTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATa                                                                  >  1:1234865/1‑149 (MQ=255)
                                                                                     tGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAggg                                                               >  2:1310116/1‑149 (MQ=255)
                                                                                        tGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCg                                                            >  2:1512719/1‑149 (MQ=255)
                                                                                            cTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAgcg                                                        >  2:328334/1‑149 (MQ=255)
                                                                                            cTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAgcg                                                        >  2:50910/1‑149 (MQ=255)
                                                                                            cTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAgcg                                                        >  2:1414492/1‑149 (MQ=255)
                                                                                             tcaatcaaTTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAgcgc                                                       >  1:1587281/1‑149 (MQ=255)
                                                                                              caatcaatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGcc                                                      <  2:1016683/149‑1 (MQ=255)
                                                                                               aatcaatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCg                                                     >  1:75800/1‑149 (MQ=255)
                                                                                                   aatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAaca                                                 >  1:1288377/1‑149 (MQ=255)
                                                                                                   aatTCGCCCAGGTACATAAACGCAGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGATGTGGGAGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCACTTAAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAaca                                                 >  1:829779/1‑149 (MQ=255)
                                                                                                    atTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAacac                                                >  2:790553/1‑149 (MQ=255)
                                                                                                    atTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAacac                                                <  1:1485382/149‑1 (MQ=255)
                                                                                                     ttCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACg                                               >  1:543779/1‑149 (MQ=255)
                                                                                                      tCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACgg                                              <  1:1310424/149‑1 (MQ=255)
                                                                                                        gCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCCCCACGGTGTAAGCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGCTCGAGCGCCGCACCCGGtt                                            >  1:237517/1‑149 (MQ=255)
                                                                                                        gCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGtt                                            >  2:571725/1‑149 (MQ=255)
                                                                                                         cccAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCGGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtc                                           <  1:524414/149‑1 (MQ=255)
                                                                                                            aGGTACATAAACGCCGTGTGTTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGTCGATCACCACGGTGTAATCCTGCGTCAGTGCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgtc                                        <  2:1262632/149‑1 (MQ=255)
                                                                                                            aGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgtc                                        <  1:681070/149‑1 (MQ=255)
                                                                                                              gTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGa                                      >  1:359063/1‑149 (MQ=255)
                                                                                                                     aaCGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagcagcag                               >  2:1571163/1‑149 (MQ=255)
                                                                                                                      aCGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagcagcagc                              >  2:1536755/1‑149 (MQ=255)
                                                                                                                        gCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCAc                            >  1:394342/1‑149 (MQ=255)
                                                                                                                          cGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACtt                          >  1:996211/1‑149 (MQ=255)
                                                                                                                              tgGTCGGAACAACGCGCAGTCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCgg                      >  2:1080138/1‑149 (MQ=255)
                                                                                                                                tcacccacaacaTGCAGCAGGCTGCGCGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGc                    >  1:1172234/7‑149 (MQ=255)
                                                                                                                                gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGc                    >  1:1014298/1‑149 (MQ=255)
                                                                                                                                gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGc                    >  2:1594277/1‑149 (MQ=255)
                                                                                                                                gTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGATCGTCGAGCAGCAGCACTTCCGGGc                    >  1:1267013/1‑149 (MQ=255)
                                                                                                                                 tCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCg                   >  2:401462/1‑149 (MQ=255)
                                                                                                                                 tCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCg                   >  1:1162026/1‑149 (MQ=255)
                                                                                                                                  acccacaacaTGCAGCAGGCTGCGCGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGa                  >  1:439349/5‑149 (MQ=255)
                                                                                                                                   ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCCCTTCCGGGCGaa                 >  2:1000739/1‑149 (MQ=255)
                                                                                                                                   ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGaa                 >  2:1549336/1‑149 (MQ=255)
                                                                                                                                    gAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAAt                <  1:893314/149‑1 (MQ=255)
                                                                                                                                    gAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCATCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAAt                >  1:205061/1‑149 (MQ=255)
                                                                                                                                     aacaacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATg               >  1:904438/1‑149 (MQ=255)
                                                                                                                                      acaacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACTCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATgg              >  2:169550/1‑149 (MQ=255)
                                                                                                                                      acaacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATgg              <  2:1288377/149‑1 (MQ=255)
                                                                                                                                        aacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCg            >  1:1264290/1‑149 (MQ=255)
                                                                                                                                                  ctgctgCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGTGATACCAcgcg  <  2:1014298/149‑1 (MQ=255)
                                                                                                                                                   tgctgcATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGa           <  2:289062/139‑1 (MQ=255)
                                                                                                                                                   tgctgcATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGa           >  1:289062/1‑139 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCG  >  NC_000913/3907579‑3907873

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