Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,469,456 C→G 8.7% T611S (ACC→AGC) ‡ mgtA → Mg(2(+)) importing P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,469,4560CG 8.7% 486.4 / 30.0 149T611S (ACC→AGC) ‡mgtAMg(2(+)) importing P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (65/71);  new base G (8/5);  total (73/76)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 151 total aligned reads displayed.

GGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGGCG  >  NC_000913/4469308‑4469605
                                                                                                                                                    |                                                                                                                                                     
 gAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCa                                                                                                                                                     >  2:1148205/1‑150 (MQ=255)
 gAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCa                                                                                                                                                     <  1:1312456/150‑1 (MQ=255)
     tACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCt                                                                                                                                                 <  1:1187099/150‑1 (MQ=255)
       ccGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGt                                                                                                                                               <  1:878974/150‑1 (MQ=255)
       ccGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGt                                                                                                                                               >  2:1339445/1‑150 (MQ=255)
       ccGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGt                                                                                                                                               >  1:1318151/1‑150 (MQ=255)
        cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGtt                                                                                                                                              >  1:1419570/1‑150 (MQ=255)
        cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGtt                                                                                                                                              >  2:1397159/1‑150 (MQ=255)
        cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGtt                                                                                                                                              <  1:1172138/150‑1 (MQ=255)
         gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                             <  1:269688/150‑3 (MQ=255)
         gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                             <  1:909187/150‑3 (MQ=255)
         gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGACATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATATTGCACAGCGTACCACGCTGttt                                                                                                                                             >  1:1040803/1‑150 (MQ=255)
          tAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgca                                                                                                                                            <  2:906944/150‑4 (MQ=255)
           aaaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaa                                                                                                                                           <  1:1325698/150‑5 (MQ=255)
            aaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGcc                                                                                                                                          <  1:276143/150‑1 (MQ=255)
              aaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCg                                                                                                                                        >  1:416479/1‑150 (MQ=255)
                tCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGcc                                                                                                                                      <  1:874673/150‑1 (MQ=255)
                        ggCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCg                                                                                                                              >  1:1259032/1‑150 (MQ=255)
                           taCAGTGAGTTAGTTGCAGCTTAAGTGTGCCATGTAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTTGTGATTTTGAAACGCTATCTGACGACGAACTGGCAATTCTTGCAGTGCGTACCTCTGTGTTTGCCGGCCTGACGCCGATg                                                                                                                           <  2:976612/149‑1 (MQ=255)
                           gACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCATACCACGCTGTTTGCCCGCCTGACGCCGATg                                                                                                                           <  1:1406493/150‑1 (MQ=255)
                                 gAGTTAGTTGCATCGAAAGTGTGCCATTAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTTATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATaaa                                                                                                                     >  2:674926/1‑150 (MQ=255)
                                 gAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATaaa                                                                                                                     >  1:341687/1‑150 (MQ=255)
                                   gttagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTTGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGa                                                                                                                   <  2:1179047/150‑1 (MQ=255)
                                   gttagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGa                                                                                                                   <  2:533007/150‑1 (MQ=255)
                                     tagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCGGCCTGACGCCGATGCATAAAGAAc                                                                                                                 <  1:666715/150‑1 (MQ=255)
                                     tagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAAc                                                                                                                 <  1:76001/150‑1 (MQ=255)
                                      agttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACg                                                                                                                <  2:636884/150‑1 (MQ=255)
                                       gttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGc                                                                                                               >  2:1018773/1‑150 (MQ=255)
                                       gttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGc                                                                                                               >  1:217377/1‑150 (MQ=255)
                                         tGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                             >  2:37673/1‑116 (MQ=255)
                                         tGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                             <  1:37673/118‑3 (MQ=255)
                                         tGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACCACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAt                                                                                                             >  2:638619/1‑150 (MQ=255)
                                          gCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAtt                                                                                                            <  1:914699/150‑1 (MQ=255)
                                           cAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtg                                                                                                           <  2:579524/150‑1 (MQ=255)
                                            aGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGcc                                                                                                                                      <  1:693067/122‑1 (MQ=255)
                                            aGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgt                                                                                                          >  1:1053504/1‑150 (MQ=255)
                                            aGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgt                                                                                                          >  1:1113434/1‑150 (MQ=255)
                                               gAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAc                                                                                                       >  1:146854/1‑150 (MQ=255)
                                                     gtgCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACt                                                                                                 <  1:1126814/150‑1 (MQ=255)
                                                       gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGa                                                                                               >  2:1049457/1‑150 (MQ=255)
                                                       gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGa                                                                                               >  1:347796/1‑150 (MQ=255)
                                                       gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGa                                                                                               >  2:700830/1‑150 (MQ=255)
                                                         cATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAg                                                                                             <  1:364610/150‑1 (MQ=255)
                                                         cATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAg                                                                                             >  1:90477/1‑150 (MQ=255)
                                                              aGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGa                                                                                        >  2:1104127/1‑150 (MQ=255)
                                                                tGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAg                                                                                      <  1:1005352/150‑1 (MQ=255)
                                                                    gCTGGACGCGGGAGAGGTGGTAATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCa                                                                                  >  2:534927/1‑150 (MQ=255)
                                                                       ggACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCAtgt                                                                               >  2:431502/1‑150 (MQ=255)
                                                                         aCGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTgg                                                                             >  2:948307/1‑150 (MQ=255)
                                                                            cgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTg                                                                          >  2:575808/1‑150 (MQ=255)
                                                                              ggAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCa                                                                                  >  1:611209/1‑140 (MQ=255)
                                                                              ggAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGc                                                                        <  1:64487/150‑1 (MQ=255)
                                                                              ggAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGc                                                                        <  2:251249/150‑1 (MQ=255)
                                                                                agagGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCtt                                                                      <  2:931060/150‑1 (MQ=255)
                                                                                   ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                                   >  2:1281048/1‑150 (MQ=255)
                                                                                   ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                                   >  2:1309009/1‑150 (MQ=255)
                                                                                   ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                                   >  1:963055/1‑150 (MQ=255)
                                                                                   ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGACGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                                   >  1:1013589/1‑150 (MQ=255)
                                                                                   ggtggtCATTCGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                                   >  2:367270/1‑150 (MQ=255)
                                                                                    gtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGTACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATg                                                                  <  1:1076632/150‑1 (MQ=255)
                                                                                    gtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATg                                                                  <  1:246787/150‑1 (MQ=255)
                                                                                      ggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                                <  1:223222/150‑1 (MQ=255)
                                                                                      ggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                                <  2:885071/150‑1 (MQ=255)
                                                                                       ctCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                               <  1:1147465/149‑1 (MQ=255)
                                                                                             ggTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGt                                                         >  2:325910/1‑150 (MQ=255)
                                                                                              gTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTa                                                        >  1:760417/1‑150 (MQ=255)
                                                                                                 gTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTa                                                     >  2:235789/1‑150 (MQ=255)
                                                                                                    atatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatg                                                  <  2:774012/150‑1 (MQ=255)
                                                                                                     tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatga                                                 >  2:917854/1‑150 (MQ=255)
                                                                                                     tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatga                                                 >  1:501790/1‑150 (MQ=255)
                                                                                                     tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatga                                                 >  2:103683/1‑150 (MQ=255)
                                                                                                      atTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgat                                                >  2:1413683/1‑150 (MQ=255)
                                                                                                       ttGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgatg                                               >  1:1269597/1‑150 (MQ=255)
                                                                                                           aaCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                           >  1:403778/1‑150 (MQ=255)
                                                                                                           aaCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                           >  2:1366217/1‑150 (MQ=255)
                                                                                                            aCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgccc                                          <  2:1148973/150‑1 (MQ=255)
                                                                                                             cGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcccg                                         >  2:815247/1‑150 (MQ=255)
                                                                                                                   cTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTAc                                   <  1:1104127/150‑1 (MQ=255)
                                                                                                                       cgacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACgcgc                               >  2:757176/1‑150 (MQ=255)
                                                                                                                           gaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCg                           >  2:669329/1‑150 (MQ=255)
                                                                                                                            aaCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACACCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCgg                          >  1:947213/1‑150 (MQ=255)
                                                                                                                               tGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGata                       >  1:1432858/1‑150 (MQ=255)
                                                                                                                                ggCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGat                        >  2:866334/1‑148 (MQ=255)
                                                                                                                                ggCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGat                        <  1:866334/148‑1 (MQ=255)
                                                                                                                                 gCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATc                     <  2:1032431/150‑1 (MQ=255)
                                                                                                                                  cAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCg                    >  1:675587/1‑150 (MQ=255)
                                                                                                                                  cAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCg                    <  2:217377/150‑1 (MQ=255)
                                                                                                                                   aaaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCgg                   <  1:1366217/150‑1 (MQ=255)
                                                                                                                                   aaaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCgg                   >  2:796706/1‑150 (MQ=255)
                                                                                                                                       ggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAtt               >  1:289231/6‑150 (MQ=255)
                                                                                                                                       ggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAtt               >  1:55129/6‑150 (MQ=255)
                                                                                                                                       ggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAtt               >  2:308119/6‑150 (MQ=255)
                                                                                                                                         tGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc             <  2:723230/150‑1 (MQ=255)
                                                                                                                                         tGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc             >  2:472423/1‑150 (MQ=255)
                                                                                                                                         tGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc             <  1:795256/150‑1 (MQ=255)
                                                                                                                                         tGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc             <  1:103683/150‑1 (MQ=255)
                                                                                                                                             ccgCGTACCACGCTGTTTGCCCGCCTGCCGCCGATGAATAAACAACGCATTGTCACCTTACTTAAGCGCGACGGGATTGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCGGCGGATATCGGCATTTCtgtg         <  2:146854/148‑1 (MQ=255)
                                                                                                                                             caGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtgtg         <  1:575808/150‑1 (MQ=255)
                                                                                                                                                   aCCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGgc   <  2:1113434/150‑1 (MQ=255)
                                                                                                                                                    ccACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGgcg  >  2:1161376/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGGCG  >  NC_000913/4469308‑4469605

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