Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,232 G→T 8.5% G628V (GGC→GTC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,2320GT 8.5% 372.4 / 16.7 117G628V (GGC→GTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (52/55);  new base T (5/5);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

Only 100 of 117 total aligned reads displayed.

GGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGC  >  NC_000913/4187088‑4187377
                                                                                                                                                |                                                                                                                                                 
ggATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTc                                                                                                                                               >  1:366786/1‑149 (MQ=255)
 gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCTCCCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCg                                                                                                                                              >  2:240494/1‑149 (MQ=255)
 gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTGTTGCGGACCAGATCATGTACACCGGCTTCg                                                                                                                                              >  6:264722/1‑149 (MQ=255)
        gTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATg                                                                                                                                        >  3:41871/1‑148 (MQ=255)
             gAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcgc                                                                                                                                  >  6:296306/1‑149 (MQ=255)
             gAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcgc                                                                                                                                  >  7:49322/1‑149 (MQ=255)
              aaaGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGcgcg                                                                                                                                 >  2:248940/1‑149 (MQ=255)
                       ccTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTgg                                                                                                                           <  8:280819/146‑1 (MQ=255)
                       ccTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTg                                                                                                                         <  4:41871/148‑1 (MQ=255)
                        cTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCa                                                                                                                       >  7:134116/1‑149 (MQ=255)
                           aCTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCAt                                                                                                                      <  4:321357/147‑1 (MQ=255)
                             tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTAGTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGt                                                                                                                  >  5:54190/1‑149 (MQ=255)
                               cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGACCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                                                >  8:282406/1‑149 (MQ=255)
                               cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                                                >  1:370376/1‑149 (MQ=255)
                               cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                                                >  7:199419/1‑149 (MQ=255)
                               cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGAGCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                                                 <  1:55544/148‑1 (MQ=255)
                                 tCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGGCTGAACCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTgg                                                                                                               >  5:384107/1‑148 (MQ=255)
                                    tCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATc                                                                                                           >  8:387463/1‑149 (MQ=255)
                                       aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTTCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                                        <  3:113010/149‑1 (MQ=255)
                                        ccAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAt                                                                                                        >  4:65365/1‑148 (MQ=255)
                                        ccAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                                       <  1:271534/149‑1 (MQ=255)
                                            gcgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACAt                                                                                                   <  8:175067/149‑1 (MQ=255)
                                             cgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATg                                                                                                  <  4:301234/149‑1 (MQ=255)
                                               cTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGt                                                                                                <  8:199419/149‑1 (MQ=255)
                                                tGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGt                                                                                                >  5:237095/1‑148 (MQ=255)
                                                   gttaaaaaTCAATCTCCAAAATGGTGTACACCTGCTACCGGATTCTTGGGCTGCAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGCGCATCTGTTGGTATCGATGACATGGTCATc                                                                                            <  3:14091/146‑1 (MQ=255)
                                                      aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCg                                                                                         <  7:387463/149‑1 (MQ=255)
                                                           gCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                                    <  5:348849/149‑1 (MQ=255)
                                                               tctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGaaa                                                                                >  2:48145/1‑149 (MQ=255)
                                                                 tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                                              >  8:146064/1‑149 (MQ=255)
                                                                   cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCg                                                                                                                                              <  1:17553/83‑1 (MQ=255)
                                                                   cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                                              <  2:251289/147‑1 (MQ=255)
                                                                   cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAac                                                                               <  8:49322/146‑1 (MQ=255)
                                                                   cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGCGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                                             <  3:17538/148‑1 (MQ=255)
                                                                   cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTAATTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                                             >  1:106700/1‑148 (MQ=255)
                                                                    aaaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                                            >  7:224981/1‑148 (MQ=255)
                                                                      aaTGCTGAACACCTGCTAGCGCATTCTCTGTCTGAAAGCGTCCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGATACACGaaa                                                                         <  2:231269/149‑1 (MQ=255)
                                                                      aaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaaa                                                                         <  8:134116/149‑1 (MQ=255)
                                                                            gAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtct                                                                   >  6:323837/1‑149 (MQ=255)
                                                                             aaCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                                                  >  6:189391/1‑149 (MQ=255)
                                                                             aaCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                                                  >  3:156074/1‑149 (MQ=255)
                                                                             aaCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                                                  <  2:77456/149‑1 (MQ=255)
                                                                                acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTcc                                                                 >  5:117041/1‑147 (MQ=255)
                                                                                 ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAg                                                              <  2:386565/149‑1 (MQ=255)
                                                                                  cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                                             <  5:347250/149‑1 (MQ=255)
                                                                                   tGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                                             >  5:316237/1‑148 (MQ=255)
                                                                                     cTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcag                                                          <  2:44441/149‑1 (MQ=255)
                                                                                     cTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                                           <  1:240738/148‑1 (MQ=255)
                                                                                      tACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACa                                                                                                    <  3:24052/106‑1 (MQ=255)
                                                                                       aCCGCATTCTCGGTCTGAACCCGACCGTTATTTTTGCGGACCAGATCATGTACACAGGCTTCGCCTATGCAGCGCGTTCTGGTGCATATGTTGGTATCGATGACATGGTCATCACGGAGAAGAAACACGAAATCATCTCCGAGgcag                                                          >  7:13441/1‑147 (MQ=255)
                                                                                        ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                                       >  8:298521/1‑149 (MQ=255)
                                                                                        ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                                       >  1:162802/1‑149 (MQ=255)
                                                                                        ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                                       >  6:3293/1‑149 (MQ=255)
                                                                                        ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaa                                                        <  1:248940/148‑1 (MQ=255)
                                                                                        ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACAGGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                                       <  6:145968/149‑1 (MQ=255)
                                                                                           cATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcag                                                    <  3:65365/149‑1 (MQ=255)
                                                                                           cATTCTCCGCCTGAAACCGACCGTTATTTTTGCGGCCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGGGCATCTGTTGGTATCGATGACATGGTCATCCCCGAGAAGAAACACGAAACCATCTCCGAGgcagacgcag                                                    <  6:117041/149‑1 (MQ=255)
                                                                                            aTTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcaga                                                   <  7:317725/149‑1 (MQ=255)
                                                                                             ttCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcagaa                                                  <  2:370376/149‑1 (MQ=255)
                                                                                              tctcGGTCTGAAGCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGcagaagcagaag                                                 >  7:264155/1‑149 (MQ=255)
                                                                                                  ggTCTGAAACCGACCGTTATTTTTGCGGACCTGATCATGTACACCGGCTTCTAATATGCAGCGCGTTCTTGTGCATGTGTTTGTATCGATGACATGGTCTTCCCTGAGAAGAAAGACGAAATCATCTCCGATGCAGAAGCAGAAGt                                                <  6:378239/146‑1 (MQ=255)
                                                                                                  ggTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                                    >  1:12165/1‑62 (MQ=255)
                                                                                                  ggTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                                    <  2:12165/62‑1 (MQ=255)
                                                                                                  ggTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGc                                             >  4:197308/1‑149 (MQ=255)
                                                                                                  ggTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGc                                             <  3:55958/149‑1 (MQ=255)
                                                                                                    tCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCt                                            >  4:152711/1‑148 (MQ=255)
                                                                                                    tCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTg                                           <  5:323837/149‑1 (MQ=255)
                                                                                                      tGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGa                                          >  5:369835/1‑148 (MQ=255)
                                                                                                        caaCCGCCCGTGATTTTTTAGGACCAGCTAATGTACACAGGCGTCGCCTATTAATAGAGTGATGGTGGATGTGTGGGTATCTATGAGATGGTCTTCCCGGAGAAGACACACGAACTCATCTCCGAGGCAGAAGCAGAAGTGGCTGACAt                                       >  8:32628/2‑149 (MQ=255)
                                                                                                        aaaCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                                        <  5:189391/148‑1 (MQ=255)
                                                                                                        aaaCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                                        <  2:366786/148‑1 (MQ=255)
                                                                                                        aaaCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTAGCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                                       <  2:245494/149‑1 (MQ=255)
                                                                                                            cgaccgTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                                      <  3:181465/98‑1 (MQ=255)
                                                                                                            cgaccgTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                                      >  4:181465/1‑98 (MQ=255)
                                                                                                            cgaccgTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   <  2:162802/149‑1 (MQ=255)
                                                                                                            cgaccgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   <  6:316237/149‑1 (MQ=255)
                                                                                                            cgaccgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   >  2:223320/1‑149 (MQ=255)
                                                                                                              accgTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   >  5:338780/1‑147 (MQ=255)
                                                                                                              accgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGACATTCAg                                   >  1:87112/1‑147 (MQ=255)
                                                                                                              accgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   >  7:164813/1‑147 (MQ=255)
                                                                                                                cgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg                                >  1:367744/1‑148 (MQ=255)
                                                                                                                 gTTATTTTTGCGGACCAGATCATGTACACCGTCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCa                              <  5:203802/149‑1 (MQ=255)
                                                                                                                  ttATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   >  6:378807/1‑143 (MQ=255)
                                                                                                                  ttATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                                   <  5:378807/143‑1 (MQ=255)
                                                                                                                  ttATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAg                             >  6:85658/1‑149 (MQ=255)
                                                                                                                     tttttGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTc                          >  5:7655/1‑149 (MQ=255)
                                                                                                                          gCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgt                      <  2:106700/148‑1 (MQ=255)
                                                                                                                            ggACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctg                   <  1:240494/149‑1 (MQ=255)
                                                                                                                               ccAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctggtc                <  8:224981/149‑1 (MQ=255)
                                                                                                                                 aGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctggtctg              <  6:384107/149‑1 (MQ=255)
                                                                                                                                   atgaTGTACACCGGCTTCGAATATGGAGCGAGTTCTGGTGCATCTGTTGGTATCTATGACATTGTCATCCCGGAGATGaaa                                                                                <  8:360963/78‑1 (MQ=255)
                                                                                                                                   aTCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGtctggtctgg             >  1:69908/1‑148 (MQ=255)
                                                                                                                                   aTCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGtctggtctgg             >  8:378857/1‑148 (MQ=255)
                                                                                                                                      aTGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCTATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTaa          <  1:391350/148‑1 (MQ=255)
                                                                                                                                          acacCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCg     >  6:274676/1‑149 (MQ=255)
                                                                                                                                           cacCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCgg    <  8:13441/149‑1 (MQ=255)
                                                                                                                                            acCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCggg   >  4:293451/1‑149 (MQ=255)
                                                                                                                                             ccGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGc  <  5:85658/149‑1 (MQ=255)
                                                                                                                                              cGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg                                >  8:15936/1‑118 (MQ=255)
                                                                                                                                              cGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg                                <  7:15936/118‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGC  >  NC_000913/4187088‑4187377

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