Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,214 A→T 67.8% D622V (GAC→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,2140AT67.8% 158.9 / 95.7 135D622V (GAC→GTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (18/26);  new base T (44/47);  total (62/73)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

CTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTC  >  NC_000913/4187067‑4187363
                                                                                                                                                   |                                                                                                                                                     
cTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTcc                                                                                                                                                     <  1:541045/150‑1 (MQ=255)
  gTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAg                                                                                                                                                   >  1:76084/1‑150 (MQ=255)
    tGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGAt                                                                                                                                                 >  2:1699384/1‑150 (MQ=255)
      gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGTGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCa                                                                                                                                               <  1:1172971/150‑1 (MQ=255)
      gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCa                                                                                                                                               <  2:401867/150‑1 (MQ=255)
      gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTAATTTTTGCGGTCCAGATCa                                                                                                                                               <  2:81485/150‑1 (MQ=255)
       ccGTGCCATTCTGTGGATGATTGTAACGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCAt                                                                                                                                              <  2:398150/150‑1 (MQ=255)
        cGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATg                                                                                                                                             <  1:999071/150‑1 (MQ=255)
            ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTaca                                                                                                                                         <  1:215672/150‑1 (MQ=255)
             cATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTACATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTacac                                                                                                                                        <  2:948783/150‑1 (MQ=255)
               ttCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCg                                                                                                                                      >  1:268737/1‑150 (MQ=255)
                  tgtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCt                                                                                                                                   >  2:1142981/1‑150 (MQ=255)
                   gtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCtt                                                                                                                                  <  1:1672995/150‑1 (MQ=255)
                     ggATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGTGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCg                                                                                                                                <  1:578959/150‑1 (MQ=255)
                     ggATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCg                                                                                                                                >  2:1395436/1‑150 (MQ=255)
                      gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGc                                                                                                                               <  2:76084/150‑1 (MQ=255)
                       atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGcc                                                                                                                              <  1:1674345/150‑1 (MQ=255)
                       atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGcc                                                                                                                              <  1:232936/150‑1 (MQ=255)
                        tgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTAACGCATTCTCGGTCTGAAACCGACCGTTATTTTTGAGGTCCATATCATGTACACCGGCTTCTCCt                                                                                                                             >  1:1094551/1‑150 (MQ=255)
                         gatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTa                                                                                                                            >  2:687335/1‑150 (MQ=255)
                              tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAg                                                                                                                       >  1:399277/1‑150 (MQ=255)
                              tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAg                                                                                                                       <  1:725264/150‑1 (MQ=255)
                               aCCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                      >  2:1443647/1‑150 (MQ=255)
                                ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                     <  2:1637680/150‑1 (MQ=255)
                                ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                     <  1:482658/150‑1 (MQ=255)
                                ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                     <  2:1362295/150‑1 (MQ=255)
                                 cGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAgcgc                                                                                                                    >  2:133965/1‑150 (MQ=255)
                                        tCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTgg                                                                                                             <  1:1329983/150‑1 (MQ=255)
                                         cTGCCGTACTCCATCGTCAACCAGGCGCTGGGTAAAACAGCAATCTCCAAAATGCTGAACACCAGCTACCGCATTCTCGGTCTGACACCGACCGTTATTTTTGCGGTCCAGATCATGGACACCGGCTTCGCCTATGCAGCGCGTTCTGGt                                                                                                            <  1:74984/150‑1 (MQ=255)
                                          tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTg                                                                                                           <  1:1379686/150‑1 (MQ=255)
                                           gCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGc                                                                                                          >  1:1671190/1‑150 (MQ=255)
                                            ccTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAACCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCa                                                                                                         >  1:1055147/1‑150 (MQ=255)
                                              ttACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATc                                                                                                       <  2:1374987/150‑1 (MQ=255)
                                               tACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCt                                                                                                      >  1:1459271/1‑150 (MQ=255)
                                               tACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCt                                                                                                      >  2:75478/1‑150 (MQ=255)
                                                aCTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTg                                                                                                     >  2:520380/1‑150 (MQ=255)
                                                aCTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACAGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTg                                                                                                     >  1:217527/1‑150 (MQ=255)
                                                  tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                                   >  1:4165/1‑150 (MQ=255)
                                                   ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                                  <  2:134512/150‑1 (MQ=255)
                                                    cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTgg                                                                                                 <  2:1459271/150‑1 (MQ=255)
                                                    cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGGTAGTTTTGCGGACCAGCTCATGTACACCGGCGTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTgg                                                                                                 >  2:1227505/1‑150 (MQ=255)
                                                     aTCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGt                                                                                                >  2:1652088/1‑150 (MQ=255)
                                                      tCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGtg                                                                                               <  1:337071/150‑2 (MQ=255)
                                                        gTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATc                                                                                             >  1:647182/1‑150 (MQ=255)
                                                            aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                         <  1:1395436/150‑1 (MQ=255)
                                                            aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                         <  1:1191448/150‑1 (MQ=255)
                                                                      gggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCAt                                                                               >  2:1136277/1‑150 (MQ=255)
                                                                      gggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCAt                                                                               >  2:1026436/1‑150 (MQ=255)
                                                                         tAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATccc                                                                            >  2:1241696/1‑150 (MQ=255)
                                                                           aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                          <  1:473454/150‑1 (MQ=255)
                                                                             aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACCATGTCATcgatacc                                                                        <  1:1522829/150‑7 (MQ=255)
                                                                              aaGCAATCTCCAAAATGCTGCACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgaga                                                                       <  1:278134/150‑1 (MQ=255)
                                                                               aGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                      <  1:1395093/150‑1 (MQ=255)
                                                                               aGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                      >  1:698013/1‑150 (MQ=255)
                                                                                  aaTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGCTCATATACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAgaagaa                                                                   <  2:1081924/150‑1 (MQ=255)
                                                                                   aTCTCCATACTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGaaa                                                                  <  2:1720294/150‑1 (MQ=255)
                                                                                   aTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCTATGACATGGTCATCCCGGAGAAGaaa                                                                  >  2:900006/1‑150 (MQ=255)
                                                                                     ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                                >  2:219237/1‑150 (MQ=255)
                                                                                     ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                                <  1:985464/150‑1 (MQ=255)
                                                                                     ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGCTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                                <  1:1604852/150‑1 (MQ=255)
                                                                                      tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                               >  1:1213255/1‑150 (MQ=255)
                                                                                       ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACg                                                              >  1:1449812/1‑150 (MQ=255)
                                                                                        cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGa                                                             >  1:1664389/1‑150 (MQ=255)
                                                                                         aaaaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGaa                                                            <  1:1336704/150‑1 (MQ=255)
                                                                                             tGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCGGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatca                                                        <  1:1142981/150‑1 (MQ=255)
                                                                                             tGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTACATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatca                                                        <  2:1664096/150‑1 (MQ=255)
                                                                                               cTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcatc                                                      >  1:692654/1‑150 (MQ=255)
                                                                                                tGAACAGCTGCTACCGCATTCTCGGTCTGATACCGACCGTTATTTTTGCGGTCCAGATGATGTACACCGGCTTGGCCTATGGAGCGCGTTCTGGTGCATCTGTTGGAATCGATGACATGGTCAACCCGGAGAAGAAACACGAACTCAtct                                                     <  2:947907/150‑1 (MQ=255)
                                                                                                tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcatca                                                     >  2:1408879/1‑149 (MQ=255)
                                                                                                     acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAg                                                <  2:1204426/150‑1 (MQ=255)
                                                                                                     acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAg                                                <  1:1223881/150‑1 (MQ=255)
                                                                                                      ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                               <  1:1381558/150‑1 (MQ=255)
                                                                                                        tGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                             >  2:654327/1‑150 (MQ=255)
                                                                                                           tACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACTTGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaa                                          >  1:113555/1‑150 (MQ=255)
                                                                                                             ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagc                                        >  2:560895/1‑150 (MQ=255)
                                                                                                             ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagc                                        >  2:109833/1‑150 (MQ=255)
                                                                                                              cGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagca                                       >  1:1309332/1‑150 (MQ=255)
                                                                                                              cGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                             >  1:1516863/1‑144 (MQ=255)
                                                                                                              cGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                             <  2:1516863/144‑1 (MQ=255)
                                                                                                                         tCTGAAACCGACCGTTATTTTTGTGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGa                            <  2:1578429/150‑1 (MQ=255)
                                                                                                                         tCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGa                            <  2:268737/150‑1 (MQ=255)
                                                                                                                          cTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaa                           >  2:726394/1‑150 (MQ=255)
                                                                                                                           tGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                          >  1:1498297/1‑150 (MQ=255)
                                                                                                                           tGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                          >  1:707500/1‑150 (MQ=255)
                                                                                                                            gAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTAAGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCACGGAGAAGAAACACGAAATGATCTCGGAGGCTGAAGCAGAAGTTGCTGAAAt                         >  1:56668/1‑150 (MQ=255)
                                                                                                                               accgaccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCa                      >  2:255293/1‑150 (MQ=255)
                                                                                                                               accgaccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCa                      <  2:774027/150‑1 (MQ=255)
                                                                                                                               accgaccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCa                      <  1:1467281/150‑1 (MQ=255)
                                                                                                                                   accgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg                  <  1:1026436/150‑1 (MQ=255)
                                                                                                                                    ccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGc                 <  1:255293/150‑1 (MQ=255)
                                                                                                                                    ccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGc                 <  1:1439609/150‑1 (MQ=255)
                                                                                                                                    cagTTTTTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCAGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGc                 <  2:1094551/148‑1 (MQ=255)
                                                                                                                                         aTTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTc            <  2:4165/150‑1 (MQ=255)
                                                                                                                                          tttttGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTcc           <  2:1449812/150‑1 (MQ=255)
                                                                                                                                               gCGGTCCAGATCATGTACACCGGCTTCGCCTATGCCGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCTTCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGGTGCTGAAATTCAGGAGCAGTTCCAgtct      >  2:1152339/1‑150 (MQ=255)
                                                                                                                                               gCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtct      >  1:65527/1‑150 (MQ=255)
                                                                                                                                               gCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtct      <  2:902640/150‑1 (MQ=255)
                                                                                                                                               gCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtct      <  2:1388651/150‑1 (MQ=255)
                                                                                                                                                cGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctg     >  2:1369601/1‑150 (MQ=255)
                                                                                                                                                  gACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctggt   >  2:1476911/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTC  >  NC_000913/4187067‑4187363

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: