Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,540 A→C 50.7% N1099T (AAC→ACC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5400AC50.7% 9.0 / 138.1 99N1099T (AAC→ACC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base C (27/23);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

Only 100 of 101 total aligned reads displayed.

AAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCG  >  NC_000913/4184395‑4184688
                                                                                                                                                 |                                                                                                                                                    
aaGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTAGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACAGTACGCCGGTAGACATCGTACTGAcccc                                                                                                                                                   >  4:142853/1‑149 (MQ=255)
   gTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCt                                                                                                                                                >  1:55391/1‑149 (MQ=255)
   gTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCt                                                                                                                                                >  7:352455/1‑149 (MQ=255)
    tataTCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGTTACGCAGTTAGACATCGTACTGAACCCGCTg                                                                                                                                               >  6:107848/1‑149 (MQ=255)
    tataTCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTg                                                                                                                                               <  8:225037/149‑1 (MQ=255)
      tatCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTgg                                                                                                                                              <  1:160032/148‑1 (MQ=255)
         cTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCg                                                                                                                                           >  3:173432/1‑148 (MQ=255)
          tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTa                                                                                                                                         >  7:155729/1‑149 (MQ=255)
            gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTATGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTgg                                                                                                                                              >  4:411800/1‑142 (MQ=255)
            gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTATGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTgg                                                                                                                                              <  3:411800/142‑1 (MQ=255)
            gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTAcc                                                                                                                                       >  2:47342/1‑149 (MQ=255)
                ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGtct                                                                                                                                   <  2:172056/149‑1 (MQ=255)
                    aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTa                                                                                                                               >  8:358453/1‑149 (MQ=255)
                        cGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGaa                                                                                                                           <  4:414697/149‑1 (MQ=255)
                        cGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGaa                                                                                                                           <  1:97778/149‑1 (MQ=255)
                          tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                         >  7:418942/1‑149 (MQ=255)
                              cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATg                                                                                                                             <  1:101970/141‑1 (MQ=255)
                              cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATg                                                                                                                             >  2:101970/1‑141 (MQ=255)
                              cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                                      <  4:173432/148‑1 (MQ=255)
                              cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                                      <  5:107848/148‑1 (MQ=255)
                               aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                                    >  7:280755/1‑149 (MQ=255)
                                       gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                            >  3:219057/1‑149 (MQ=255)
                                        aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                           <  7:358453/149‑1 (MQ=255)
                                        aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                           <  1:26506/149‑1 (MQ=255)
                                          aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                                          >  1:263771/1‑148 (MQ=255)
                                          aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                                         <  8:22744/149‑1 (MQ=255)
                                          aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                                          >  3:122631/1‑148 (MQ=255)
                                            gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTTGaaa                                                                                                       <  2:55391/149‑1 (MQ=255)
                                             aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                                      >  1:229178/1‑149 (MQ=255)
                                             aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                                      >  3:230919/1‑149 (MQ=255)
                                               ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                                     >  5:52769/1‑148 (MQ=255)
                                                gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                                   >  2:265886/1‑149 (MQ=255)
                                                  aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGCTCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAGCCTCGAAAcccccc                                                                                                 >  4:52814/1‑149 (MQ=255)
                                                   ggTCGTCACGGTAACAAGGGTGTAATTTCTAAGATAAACCCGATCGAAGATATGCCTTAAGATGAAAACGGTACTCAGTTAGACATCGTACTGAACCCTCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAACCTCGAAACCCACCt                                                                                                >  1:108259/1‑149 (MQ=255)
                                                       gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACCGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                             >  3:64920/1‑148 (MQ=255)
                                                        tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAAGGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTa                                                                                           <  3:142853/149‑1 (MQ=255)
                                                        tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                              <  4:238380/146‑1 (MQ=255)
                                                        tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                              >  3:238380/1‑146 (MQ=255)
                                                         cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTAAGATTAAAACGGTACGCCGGTAGACATCGTACTGATCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                             <  7:414062/146‑1 (MQ=255)
                                                          aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                         >  6:149226/1‑149 (MQ=255)
                                                           cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                                        <  6:52769/149‑1 (MQ=255)
                                                              tAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                                                                     >  5:188132/1‑149 (MQ=255)
                                                                  aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                                  >  2:459712/1‑148 (MQ=255)
                                                                  aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                                  >  8:146553/1‑148 (MQ=255)
                                                                   aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                                <  1:47342/149‑1 (MQ=255)
                                                                    gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATTATAACGGTACGCCGGTAGACTTCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                               <  5:180569/149‑1 (MQ=255)
                                                                    gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                               <  8:180426/149‑1 (MQ=255)
                                                                    gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                               <  7:17235/149‑1 (MQ=255)
                                                                       tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTa                                                                            >  7:142201/1‑149 (MQ=255)
                                                                           aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                               >  5:243080/1‑142 (MQ=255)
                                                                           aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                               <  6:243080/142‑1 (MQ=255)
                                                                             ttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                                                                       >  3:195652/1‑148 (MQ=255)
                                                                               cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                                     <  2:263771/148‑1 (MQ=255)
                                                                               cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                                     <  8:155729/148‑1 (MQ=255)
                                                                                tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTa                                                                                                                                         >  6:301606/1‑79 (MQ=255)
                                                                                tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTa                                                                                                                                         <  5:301606/79‑1 (MQ=255)
                                                                                tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                                >  1:425545/1‑136 (MQ=255)
                                                                                tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                                <  2:425545/136‑1 (MQ=255)
                                                                                  aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                                    >  2:261109/1‑146 (MQ=255)
                                                                                  aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                                 <  4:219057/149‑1 (MQ=255)
                                                                                    aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                               <  3:324682/149‑1 (MQ=255)
                                                                                        aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                           >  1:212861/1‑149 (MQ=255)
                                                                                        aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACCGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                           <  4:64920/149‑1 (MQ=255)
                                                                                         cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACg                                                          <  8:418942/149‑1 (MQ=255)
                                                                                         cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACg                                                          <  5:343055/149‑1 (MQ=255)
                                                                                             aTCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAt                                                      <  8:142201/149‑1 (MQ=255)
                                                                                              tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACACACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                                     <  1:197863/149‑1 (MQ=255)
                                                                                                gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt                                                   >  2:416036/1‑149 (MQ=255)
                                                                                                 aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAACCACCTGGGTATGGGTGCGAAAGGTATCGGCGACAAGATCAACGCCAAGcg                                                   >  6:281979/1‑147 (MQ=255)
                                                                                                 aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCGCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTg                                                  <  4:230919/149‑1 (MQ=255)
                                                                                                       tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAca                                             <  2:229178/148‑1 (MQ=255)
                                                                                                         ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                          <  8:280755/149‑1 (MQ=255)
                                                                                                         ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                           <  1:459712/148‑1 (MQ=255)
                                                                                                         ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                            >  2:77505/1‑147 (MQ=255)
                                                                                                           ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCGACGCCATGCTGAAAcagcag                                         >  4:206870/1‑148 (MQ=255)
                                                                                                           ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                                         >  5:111261/1‑148 (MQ=255)
                                                                                                           ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                                <  3:134793/109‑1 (MQ=255)
                                                                                                           ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                                >  4:134793/1‑109 (MQ=255)
                                                                                                           ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                        <  4:122631/149‑1 (MQ=255)
                                                                                                              cGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                                      <  6:304145/148‑1 (MQ=255)
                                                                                                               gATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                                     >  4:16069/1‑148 (MQ=255)
                                                                                                               gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                                        >  3:249546/1‑129 (MQ=255)
                                                                                                               gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                                        <  4:249546/129‑1 (MQ=255)
                                                                                                                aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                                   <  4:118230/149‑1 (MQ=255)
                                                                                                                     aaCGGTACGCCGGTAGACATCGTACTGCCCCCGCTGGGCGTCCCGTCTCGTATGCACATCGGTCAGATCCTCGCAACCCACCTGGGTATGGCGGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                              >  7:376373/1‑149 (MQ=255)
                                                                                                                     aaCGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                              >  3:228071/1‑149 (MQ=255)
                                                                                                                             gCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgc                      >  5:128405/1‑149 (MQ=255)
                                                                                                                                ggTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc                    >  8:235838/1‑148 (MQ=255)
                                                                                                                                 gTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCAGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                        <  2:258993/127‑1 (MQ=255)
                                                                                                                                 gTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCAGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                        >  1:258993/1‑127 (MQ=255)
                                                                                                                                      cATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa             <  4:228071/149‑1 (MQ=255)
                                                                                                                                       aTCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt            <  5:281979/149‑1 (MQ=255)
                                                                                                                                        tCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGATCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc           <  6:128405/149‑1 (MQ=255)
                                                                                                                                        tCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGCTCACATCCTCGAAACCCACCAGGCTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACACCAGCAAGAAGTCGCGAAACTGCGCGAATTCaaa           >  6:189148/1‑147 (MQ=255)
                                                                                                                                          gTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAATTCGCGCAGTTTCGCGACTTcttgc          <  7:470159/148‑5 (MQ=255)
                                                                                                                                          gTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCa         <  1:77505/149‑1 (MQ=255)
                                                                                                                                             cTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcg      >  7:42580/1‑149 (MQ=255)
                                                                                                                                             cTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcg      >  4:157731/1‑149 (MQ=255)
                                                                                                                                              tGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgt     >  1:42146/1‑149 (MQ=255)
                                                                                                                                                 aCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgtgcg  >  2:113542/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCG  >  NC_000913/4184395‑4184688

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