Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,214 A→G 92.7% D622G (GAC→GGC)  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,2140AG92.7% 263.0 / 6.2 83D622G (GAC→GGC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (34/42);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTC  >  NC_000913/4187069‑4187358
                                                                                                                                                 |                                                                                                                                                
gTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAg                                                                                                                                              >  1:157383/1‑150 (MQ=255)
gTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAg                                                                                                                                              >  2:1251912/1‑150 (MQ=255)
gTTGGCCGGGGCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAACGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTATGAAACCGACCGTTATTTTTGCGGGCCAg                                                                                                                                              <  2:550934/150‑1 (MQ=255)
    gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAACTGCTGAACACCTGCTACCGCATTCTCGGGCTGAAACCGACCGTTATTTTTGCGGACCAGATCa                                                                                                                                          <  2:183640/150‑1 (MQ=255)
     ccGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCAt                                                                                                                                         >  1:715306/1‑150 (MQ=255)
        ggCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTa                                                                                                                                      <  1:1157929/149‑1 (MQ=255)
         gCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTac                                                                                                                                     <  1:534809/150‑1 (MQ=255)
           cATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTacac                                                                                                                                   >  2:25123/1‑150 (MQ=255)
           cATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTacac                                                                                                                                   >  2:640086/1‑150 (MQ=255)
            aTTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACAcc                                                                                                                                  <  2:938998/150‑1 (MQ=255)
            aTTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGGCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACAcc                                                                                                                                  <  2:463888/150‑1 (MQ=255)
                     atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGAACACCGGCTTCGcc                                                                                                                         <  2:325800/150‑1 (MQ=255)
                       gatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTa                                                                                                                       <  1:618601/150‑1 (MQ=255)
                          tGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGc                                                                                                                    <  1:133256/150‑1 (MQ=255)
                          tGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGc                                                                                                                    <  2:120976/150‑1 (MQ=255)
                           gTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGACACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCa                                                                                                                   >  1:598708/1‑150 (MQ=255)
                           gTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCa                                                                                                                   >  2:919165/1‑150 (MQ=255)
                             aCCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGACCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                 >  1:697325/1‑150 (MQ=255)
                                  aaGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGtt                                                                                                            <  1:235542/150‑1 (MQ=255)
                                      tCTGCCTAACTCCATCCTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCATATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTgg                                                                                                        <  2:11507/150‑1 (MQ=255)
                                        tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTg                                                                                                      <  2:269736/150‑1 (MQ=255)
                                         gCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGc                                                                                                     >  1:894581/1‑150 (MQ=255)
                                         gCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCCCGATCTGGTGc                                                                                                     >  1:1230595/1‑150 (MQ=255)
                                         gCCTTACTCCATCGTCAACCAGGCGCTGCGTAAAAAAGCAATCTCCAAAATGCTGACCACCTGCTCCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGc                                                                                                     >  1:1006811/1‑150 (MQ=255)
                                           cTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCAt                                                                                                   <  1:677612/150‑1 (MQ=255)
                                            ttACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGTGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATc                                                                                                  <  2:1314544/150‑1 (MQ=255)
                                              aCTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTg                                                                                                >  1:606092/1‑150 (MQ=255)
                                                tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                              <  2:1032375/150‑1 (MQ=255)
                                                 ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                             >  2:253184/1‑150 (MQ=255)
                                                 ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                             <  2:263435/150‑1 (MQ=255)
                                                 ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                             <  2:551338/150‑1 (MQ=255)
                                                        cAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGCCCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTGTGGTGCATCTGTTGGGATCGa                                                                                      >  2:814554/1‑150 (MQ=255)
                                                          aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                    <  2:345738/150‑1 (MQ=255)
                                                            cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGCGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGaa                                                                                  <  2:188981/150‑2 (MQ=255)
                                                                cgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGGCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATgg                                                                              <  2:1049420/150‑1 (MQ=255)
                                                                cgcTGGGTAAAAAAGCAAACTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCACCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATgg                                                                              <  2:502502/150‑1 (MQ=255)
                                                                 gcTGTGTAAAAAACCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGt                                                                             <  2:187021/150‑1 (MQ=255)
                                                                    gggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCAt                                                                          >  2:318993/1‑150 (MQ=255)
                                                                      gTAAATAAGCAATCTCCAAAATGCTCAACACCTGCGACCGCATGCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATcc                                                                        <  2:184552/150‑1 (MQ=255)
                                                                       tAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATccc                                                                       >  1:355455/1‑150 (MQ=255)
                                                                         aaaaaGCAATCTCCACAATGCTGTTCACGTGCTACCTCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCATATCTTTTAAACCGGCTTGGCCTATGCTGCTCGTTCTGGTGCATCTGTTGGTATCGATGACATTGTCATCCCgg                                                                     <  2:719803/150‑1 (MQ=255)
                                                                         aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                     >  2:827407/1‑150 (MQ=255)
                                                                         aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                     <  2:506/150‑1 (MQ=255)
                                                                          aaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGACACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGGTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGga                                                                    <  1:1310561/150‑1 (MQ=255)
                                                                           aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTAGACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                   <  2:371525/150‑1 (MQ=255)
                                                                           aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                   <  2:1216150/150‑1 (MQ=255)
                                                                              gCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaag                                                                >  2:517071/1‑150 (MQ=255)
                                                                                  tctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCGGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAac                                                            <  2:1175888/150‑1 (MQ=255)
                                                                                   ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                           <  1:1027860/150‑1 (MQ=255)
                                                                                   ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                           >  2:788725/1‑150 (MQ=255)
                                                                                         aaTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAat                                                     <  2:752532/150‑1 (MQ=255)
                                                                                           tGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGTGCGTTCTGGTGCATCTGTTGGTATCGATGACATCGTCATCCCGGAGAAGAAACACGAAatca                                                   >  2:856889/1‑150 (MQ=255)
                                                                                             cTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcatc                                                 >  1:1306030/1‑150 (MQ=255)
                                                                                             cTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTAGTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGCAGCAACACGAGatcatc                                                 >  1:85083/1‑150 (MQ=255)
                                                                                              tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtct                                                <  2:1266974/150‑1 (MQ=255)
                                                                                               gAAGACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                               <  2:1136845/150‑1 (MQ=255)
                                                                                               gAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                               <  1:1264705/150‑1 (MQ=255)
                                                                                                aaCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTcc                                              >  1:373001/1‑150 (MQ=255)
                                                                                                aaCACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTcc                                              <  1:935828/150‑1 (MQ=255)
                                                                                                  cacCTGCTACCGCTTTCTCGGTCTGAACCCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCTGCGCGTTCTGGTGCTTCTGTTGGTATCGATGACATGGTCATCGCTGAGAAGAAACACGATATCATCTCCGa                                            <  1:984819/150‑1 (MQ=255)
                                                                                                   acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAg                                           <  2:1200850/150‑1 (MQ=255)
                                                                                                    ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          >  1:315983/1‑150 (MQ=255)
                                                                                                    ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          >  1:902808/1‑150 (MQ=255)
                                                                                                    ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                          >  1:745290/1‑150 (MQ=255)
                                                                                                    ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTACGAgg                                          >  2:338963/1‑150 (MQ=255)
                                                                                                     cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgc                                         >  2:1179855/1‑150 (MQ=255)
                                                                                                     cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgc                                         <  1:640086/150‑1 (MQ=255)
                                                                                                          aCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGTGAAGAAACACGAAATCATCTCCGAGgcagaag                                    >  1:105766/1‑150 (MQ=255)
                                                                                                                ttCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAACTCCGAGGcagaagcagaag                              >  1:339947/1‑150 (MQ=255)
                                                                                                                   tcGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTg                           <  2:804675/150‑1 (MQ=255)
                                                                                                                    cGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGc                          <  2:295958/150‑1 (MQ=255)
                                                                                                                      gTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGGCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTg                        <  2:718016/150‑1 (MQ=255)
                                                                                                                       tCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGa                       <  2:553546/150‑1 (MQ=255)
                                                                                                                        cTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaa                      >  2:568920/1‑150 (MQ=255)
                                                                                                                         tGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                     >  1:787507/1‑150 (MQ=255)
                                                                                                                          gAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                    >  2:16258/1‑150 (MQ=255)
                                                                                                                          gAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                    >  2:412920/1‑150 (MQ=255)
                                                                                                                          gAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                    >  2:142849/1‑150 (MQ=255)
                                                                                                                          gAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaatt                    >  2:1247764/1‑148 (MQ=255)
                                                                                                                            aaCCGACCGGTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTc                  <  2:333551/150‑1 (MQ=255)
                                                                                                                               cgaccgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAgg               <  2:721166/150‑1 (MQ=255)
                                                                                                                                  ccgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAGGACATGGTCATGCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGc            <  2:450263/150‑1 (MQ=255)
                                                                                                                                            tGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtc  <  2:606092/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTC  >  NC_000913/4187069‑4187358

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