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% 561.9 / ‑4.6 150Y1285C (TAC→TGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (63/86);  total (63/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 152 total aligned reads displayed.

AACCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGG  >  NC_000913/4184950‑4185245
                                                                                                                                                    |                                                                                                                                                   
aaCCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAACCATACGGCGCAGCATgca                                                                                                                                                   <  1:289390/151‑1 (MQ=255)
 accaccTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATgcac                                                                                                                                                  >  2:585670/1‑151 (MQ=255)
     ccTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTg                                                                                                                                              <  2:1078372/151‑1 (MQ=255)
         gTCGACGACAAGATGCACGCGCGTTCCCCCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAgg                                                                                                                                          <  4:260653/151‑1 (MQ=255)
         gTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAgg                                                                                                                                          >  4:322986/1‑151 (MQ=255)
             acgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAAt                                                                                                                                      <  1:470686/151‑1 (MQ=255)
             acgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAAt                                                                                                                                      <  2:523302/151‑1 (MQ=255)
              cgacAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAAt                                                                                                                                      <  4:556136/150‑1 (MQ=255)
                acAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGATGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGc                                                                                                                                    >  2:17478/1‑150 (MQ=255)
                  aaGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTc                                                                                                                                  >  2:1695395/1‑150 (MQ=255)
                    gATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCAcc                                                                                                                               >  1:1589550/1‑151 (MQ=255)
                     aTGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCg                                                                                                                              >  2:72117/1‑151 (MQ=255)
                     aTGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCg                                                                                                                              >  3:545244/1‑151 (MQ=255)
                     aTGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCg                                                                                                                              <  4:16495/151‑1 (MQ=255)
                       gCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGtt                                                                                                                            <  3:38082/151‑1 (MQ=255)
                       gCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGtt                                                                                                                            <  1:1357657/151‑1 (MQ=255)
                             gcgTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCt                                                                                                                      >  3:218939/1‑151 (MQ=255)
                             gcgTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCt                                                                                                                      <  3:275556/151‑1 (MQ=255)
                               gTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtg                                                                                                                     >  1:1372815/1‑150 (MQ=255)
                                ttCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGAGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtgat                                                                                                                   >  2:186967/1‑151 (MQ=255)
                                  ccACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtgatga                                                                                                                 <  2:662742/151‑1 (MQ=255)
                                  ccACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCtgatg                                                                                                                  >  2:1012221/1‑150 (MQ=255)
                                   cACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGAAGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGAc                                                                                                                <  3:422308/151‑1 (MQ=255)
                                     ccGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGt                                                                                                              <  2:1162989/151‑1 (MQ=255)
                                       ggTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTg                                                                                                             >  1:1691130/1‑150 (MQ=255)
                                        gTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGa                                                                                                            >  4:354474/1‑150 (MQ=255)
                                          tCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAAc                                                                                                          <  2:152135/150‑1 (MQ=255)
                                          tCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAAc                                                                                                          <  1:772512/150‑1 (MQ=255)
                                              aCAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTc                                                                                                      <  2:130226/150‑1 (MQ=255)
                                              aCAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCg                                                                                                     >  2:948842/1‑151 (MQ=255)
                                              aCAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCg                                                                                                     <  3:567768/151‑1 (MQ=255)
                                                aGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGCGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGt                                                                                                    >  1:35781/1‑150 (MQ=255)
                                                aGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGt                                                                                                    >  4:108724/1‑150 (MQ=255)
                                                aGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTa                                                                                                   <  2:713703/151‑1 (MQ=255)
                                                aGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTa                                                                                                   <  2:958442/151‑1 (MQ=255)
                                                        tACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATg                                                                                           >  2:1736999/1‑151 (MQ=255)
                                                         aCTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATg                                                                                           <  1:67308/150‑1 (MQ=255)
                                                            cagcagCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGtat                                                                                        <  1:211512/150‑1 (MQ=255)
                                                             agcagcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATaa                                                                                      <  2:545059/151‑1 (MQ=255)
                                                              gcagcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATaaa                                                                                     <  3:209266/151‑1 (MQ=255)
                                                               cagcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACg                                                                                                               >  4:289922/1‑124 (MQ=255)
                                                               cagcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACg                                                                                                               <  3:289940/124‑1 (MQ=255)
                                                               cagcCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATaaaa                                                                                    <  4:429906/151‑1 (MQ=255)
                                                                   cGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACAt                                                                                <  4:565542/151‑1 (MQ=255)
                                                                    gCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATc                                                                               >  1:374023/1‑151 (MQ=255)
                                                                       gggTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCg                                                                              <  2:1151803/149‑1 (MQ=255)
                                                                       gggTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGt                                                                             <  4:520751/150‑1 (MQ=255)
                                                                        ggtggtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTgg                                                                           >  1:528053/1‑151 (MQ=255)
                                                                         gtggtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGa                                                                          <  3:148779/151‑1 (MQ=255)
                                                                            gtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACgg                                                                       >  2:1714038/1‑151 (MQ=255)
                                                                            gtAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACgg                                                                       >  2:346091/1‑151 (MQ=255)
                                                                              aaGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCa                                                                     <  1:1214359/151‑1 (MQ=255)
                                                                               aGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCa                                                                     >  2:1271567/1‑150 (MQ=255)
                                                                                  cacaGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAAcc                                                                  >  2:1103660/1‑150 (MQ=255)
                                                                                   acaGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCAt                                                                >  2:172710/1‑151 (MQ=255)
                                                                                       ttCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATc                                                               >  2:1726315/1‑148 (MQ=255)
                                                                                        tCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATc                                                               >  2:560703/1‑147 (MQ=255)
                                                                                        tCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGa                                                            <  2:1714105/150‑1 (MQ=255)
                                                                                           gtggtCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGa                                                        <  1:160054/151‑1 (MQ=255)
                                                                                            tggtCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATgg                                                         >  1:1441714/1‑149 (MQ=255)
                                                                                               tCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCg                                                    >  2:1527585/1‑151 (MQ=255)
                                                                                               tCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCg                                                    >  1:472518/1‑151 (MQ=255)
                                                                                                   cGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCa                                                <  1:774052/151‑1 (MQ=255)
                                                                                                    gTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCAt                                               <  2:213162/151‑1 (MQ=255)
                                                                                                    gTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCAt                                               <  1:620428/151‑1 (MQ=255)
                                                                                                       tCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGc                                             <  2:961537/150‑1 (MQ=255)
                                                                                                         ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGcc                                            >  1:939230/1‑149 (MQ=255)
                                                                                                         ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAg                                          <  2:265776/151‑1 (MQ=255)
                                                                                                              gaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATcc                                     <  2:454152/151‑1 (MQ=255)
                                                                                                               aTGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCt                                    <  2:1137182/151‑1 (MQ=255)
                                                                                                                 ggAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCtt                                   <  1:811897/150‑1 (MQ=255)
                                                                                                                 ggAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCtgc                                  <  1:817930/151‑3 (MQ=255)
                                                                                                                 ggAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTc                                  >  3:169511/1‑151 (MQ=255)
                                                                                                                  gAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCa                                 <  2:819996/151‑1 (MQ=255)
                                                                                                                   aaGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCaa                                <  2:65328/151‑1 (MQ=255)
                                                                                                                    aGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAAc                               >  2:1426004/1‑151 (MQ=255)
                                                                                                                    aGTGTGGGCGCTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAAc                               <  2:1044773/151‑1 (MQ=255)
                                                                                                                           gcgcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgt                        <  1:948724/151‑1 (MQ=255)
                                                                                                                           gcgcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgt                        <  2:51081/151‑1 (MQ=255)
                                                                                                                           gcgcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgt                        <  1:844196/151‑1 (MQ=255)
                                                                                                                            cgcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgtt                       <  1:1406762/151‑1 (MQ=255)
                                                                                                                             gcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgttg                      >  1:1220635/1‑151 (MQ=255)
                                                                                                                             gcTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgttg                      <  2:1739765/151‑1 (MQ=255)
                                                                                                                              cTGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgttg                      >  1:895675/1‑150 (MQ=255)
                                                                                                                               tGGAAGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGa                     >  1:364004/1‑150 (MQ=255)
                                                                                                                                  aaGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAATAACATCGTGGACGGCAACCATCAGATTGAGCCGGGCATGACAGAATCCTTCAACGTATTGTTGAAaga                 <  1:1466774/151‑1 (MQ=255)
                                                                                                                                  aaGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAaga                 <  2:411955/151‑1 (MQ=255)
                                                                                                                                   aGCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAag                  >  1:1231130/1‑149 (MQ=255)
                                                                                                                                    gCATACGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAag                  >  2:138547/1‑148 (MQ=255)
                                                                                                                                        aCGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttcg           <  2:860385/151‑1 (MQ=255)
                                                                                                                                        aCGGCGCAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttc            <  4:62221/150‑1 (MQ=255)
                                                                                                                                           gcgcAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttcgttc        <  4:258175/151‑1 (MQ=255)
                                                                                                                                            cgcAGCATGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttcgttcg       >  1:1517660/1‑151 (MQ=255)
                                                                                                                                             gcagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGc      <  1:17468/151‑1 (MQ=255)
                                                                                                                                             gcagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGc      <  4:516535/151‑1 (MQ=255)
                                                                                                                                             gcagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGc      <  1:1273729/151‑1 (MQ=255)
                                                                                                                                              cagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCt     >  2:1067942/1‑151 (MQ=255)
                                                                                                                                              cagcaTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCt     >  1:358024/1‑151 (MQ=255)
                                                                                                                                                 caTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg  >  1:1172058/1‑151 (MQ=255)
                                                                                                                                                 caTGCACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg  <  1:1103525/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTGGGTGGTAAGGCACAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGG  >  NC_000913/4184950‑4185245

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