Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,185,098 A→G Y1285C (TAC→TGC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,0980AG99.3% 538.4 / ‑4.0 140Y1285C (TAC→TGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (68/71);  total (68/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 140 total aligned reads displayed.

CCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGG  >  NC_000913/4184952‑4185244
                                                                                                                                                  |                                                                                                                                                  
ccaccTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCAcc                                                                                                                                                <  2:1526307/151‑1 (MQ=255)
  accTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCt                                                                                                                                              <  1:1214866/151‑1 (MQ=255)
   ccTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTg                                                                                                                                             <  1:586515/151‑1 (MQ=255)
      ggTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCa                                                                                                                                           <  2:412912/150‑1 (MQ=255)
       gTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGc                                                                                                                                            <  1:1621027/148‑1 (MQ=255)
         cgacgaCAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGa                                                                                                                                        <  1:1573775/150‑1 (MQ=255)
         cgacgaCAAGATGCACGCGCGTTCCACCGGTTCTTACAGACTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGCGGGAGCTGGAAGCATACGGCGCAGCATGCACCCTGCAgg                                                                                                                                         >  1:1138025/1‑149 (MQ=255)
          gacgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGaa                                                                                                                                       >  2:1083437/1‑150 (MQ=255)
          gacgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGaa                                                                                                                                       >  2:688687/1‑150 (MQ=255)
           acgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAAt                                                                                                                                     <  1:89315/151‑1 (MQ=255)
            cgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATg                                                                                                                                    <  2:829425/151‑1 (MQ=255)
            cgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATg                                                                                                                                    <  2:1127197/151‑1 (MQ=255)
             gacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGc                                                                                                                                   >  2:185975/1‑151 (MQ=255)
             gacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGc                                                                                                                                   >  1:295545/1‑151 (MQ=255)
              acAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCt                                                                                                                                  <  1:627403/151‑1 (MQ=255)
              acAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCt                                                                                                                                  <  2:406670/151‑1 (MQ=255)
               cAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTc                                                                                                                                 >  1:1491385/1‑151 (MQ=255)
               cAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTc                                                                                                                                 >  1:503223/1‑151 (MQ=255)
                aaGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTc                                                                                                                                 >  2:885202/1‑150 (MQ=255)
                aaGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCa                                                                                                                                <  2:1302779/151‑1 (MQ=255)
                  gATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCAcc                                                                                                                              >  1:1704935/1‑151 (MQ=255)
                   aTGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCg                                                                                                                             >  2:640487/1‑151 (MQ=255)
                      cACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGtt                                                                                                                           >  1:754588/1‑150 (MQ=255)
                      cACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGtt                                                                                                                           >  2:1387085/1‑150 (MQ=255)
                        cgcgcgTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTa                                                                                                                          >  2:944648/1‑149 (MQ=255)
                        cgcgcgTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAg                                                                                                                        <  2:367440/151‑1 (MQ=255)
                         gcgcgTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGt                                                                                                                       <  1:547686/151‑1 (MQ=255)
                            cgTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtg                                                                                                                    <  2:497818/151‑1 (MQ=255)
                             gTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtga                                                                                                                   <  2:772931/151‑1 (MQ=255)
                               tCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtgat                                                                                                                  <  2:437191/150‑1 (MQ=255)
                                ccACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtgatga                                                                                                                <  2:1136497/151‑1 (MQ=255)
                                ccACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtgatg                                                                                                                 >  1:438644/1‑150 (MQ=255)
                                 cACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGAc                                                                                                               >  2:273643/1‑151 (MQ=255)
                                 cACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGAc                                                                                                               >  1:120291/1‑151 (MQ=255)
                                   ccGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACg                                                                                                              <  2:1162770/150‑1 (MQ=255)
                                     ggTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGa                                                                                                           <  1:1130233/151‑1 (MQ=255)
                                       ttCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAAc                                                                                                         >  2:1182330/1‑151 (MQ=255)
                                         cTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACgg                                                                                                       >  2:16000/1‑151 (MQ=255)
                                           tACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGt                                                                                                      <  2:759101/150‑1 (MQ=255)
                                           tACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTc                                                                                                     >  2:1337556/1‑151 (MQ=255)
                                             cAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGt                                                                                                   >  1:1215778/1‑151 (MQ=255)
                                             cAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGt                                                                                                   <  1:609763/151‑1 (MQ=255)
                                              aGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTa                                                                                                  <  1:994673/151‑1 (MQ=255)
                                               gCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTa                                                                                                  <  1:1446799/150‑1 (MQ=255)
                                               gCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTAc                                                                                                 <  2:547299/151‑1 (MQ=255)
                                                  tGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCa                                                                                               >  1:1632462/1‑150 (MQ=255)
                                                    gTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGa                                                                                            <  2:1599726/151‑1 (MQ=255)
                                                        cTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGt                                                                                         >  2:1376353/1‑150 (MQ=255)
                                                           agcagcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGtata                                                                                      >  1:956627/1‑150 (MQ=255)
                                                               gcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAAc                                                                                 <  1:1602246/151‑1 (MQ=255)
                                                                 cGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACAt                                                                               <  1:180172/151‑1 (MQ=255)
                                                                 cGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACAt                                                                               <  2:1016751/151‑1 (MQ=255)
                                                                    tGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGt                                                                            >  1:383929/1‑151 (MQ=255)
                                                                     gggTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCg                                                                             <  1:686048/149‑1 (MQ=255)
                                                                      ggtggtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTgg                                                                          <  1:152505/151‑1 (MQ=255)
                                                                       gtggtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGa                                                                         <  1:78334/151‑1 (MQ=255)
                                                                        tggtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGa                                                                         >  1:1565440/1‑150 (MQ=255)
                                                                        tggtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGAc                                                                        >  2:923808/1‑151 (MQ=255)
                                                                          gtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACgg                                                                      <  1:441739/151‑1 (MQ=255)
                                                                            aaGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGc                                                                     >  2:1115453/1‑150 (MQ=255)
                                                                            aaGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGc                                                                     >  1:110812/1‑150 (MQ=255)
                                                                            aaGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCa                                                                    <  2:1631098/151‑1 (MQ=255)
                                                                                  caGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATc                                                              >  1:1427883/1‑151 (MQ=255)
                                                                                      tCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATc                                                              >  1:1360811/1‑147 (MQ=255)
                                                                                         gtggtCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGa                                                       <  1:1007438/151‑1 (MQ=255)
                                                                                           ggtCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAg                                                      <  1:498830/150‑1 (MQ=255)
                                                                                            gtCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGcc                                                    >  2:315926/1‑151 (MQ=255)
                                                                                             tCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACAGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCg                                                   >  1:1380107/1‑151 (MQ=255)
                                                                                                gCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGc                                                <  2:615963/151‑1 (MQ=255)
                                                                                                 cGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGc                                                >  1:1194244/1‑150 (MQ=255)
                                                                                                 cGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGc                                                >  2:1627593/1‑150 (MQ=255)
                                                                                                  gTTTCGGGGAGATGGAAGTGTGGGCTCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCAt                                              <  1:120507/151‑1 (MQ=255)
                                                                                                  gTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCAt                                              <  1:598409/151‑1 (MQ=255)
                                                                                                     tCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGcc                                           >  1:1333002/1‑151 (MQ=255)
                                                                                                      cGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGcc                                           >  1:601840/1‑150 (MQ=255)
                                                                                                      cGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGcc                                           >  1:1678725/1‑150 (MQ=255)
                                                                                                         ggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGaa                                       <  1:681742/151‑1 (MQ=255)
                                                                                                           agaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATc                                     >  1:298932/1‑151 (MQ=255)
                                                                                                            gaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATc                                     <  1:372342/150‑1 (MQ=255)
                                                                                                             aTGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCt                                   >  2:1684520/1‑151 (MQ=255)
                                                                                                               ggAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTc                                 >  1:634711/1‑151 (MQ=255)
                                                                                                                 aaGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCaa                               <  2:878271/151‑1 (MQ=255)
                                                                                                                    tgtgGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGt                            >  1:712554/1‑151 (MQ=255)
                                                                                                                    tgtgGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGt                            >  2:1074937/1‑151 (MQ=255)
                                                                                                                     gtgGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGt                            >  1:22615/1‑150 (MQ=255)
                                                                                                                          cgcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgtt                      <  2:1483230/151‑1 (MQ=255)
                                                                                                                          cgcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgtt                      >  2:1320889/1‑151 (MQ=255)
                                                                                                                              ggAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGaaa                  <  2:209511/151‑1 (MQ=255)
                                                                                                                              ggAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGaa                   >  1:621786/1‑150 (MQ=255)
                                                                                                                                aaGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAaga                <  1:944575/151‑1 (MQ=255)
                                                                                                                                aaGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAaga                <  1:1632188/151‑1 (MQ=255)
                                                                                                                                   cATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAag                 >  1:159100/1‑147 (MQ=255)
                                                                                                                                   cATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAt             <  2:237989/151‑1 (MQ=255)
                                                                                                                                           gcagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGc     >  1:1060559/1‑151 (MQ=255)
                                                                                                                                            cagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttcgttcg      >  1:1472070/1‑149 (MQ=255)
                                                                                                                                            cagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCt    >  1:441530/1‑151 (MQ=255)
                                                                                                                                            cagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCt    <  2:1632564/151‑1 (MQ=255)
                                                                                                                                             agcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTg   >  1:1687212/1‑151 (MQ=255)
                                                                                                                                              gcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTgg  >  2:148020/1‑151 (MQ=255)
                                                                                                                                              gcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTgg  <  1:1675523/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGG  >  NC_000913/4184952‑4185244

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