Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,214 A→G 23.6% 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,2140AG23.6% 152.9 / 49.9 89D622G (GAC→GGC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base G (11/10);  total (48/41)
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.95e-01

CTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGT  >  NC_000913/4187067‑4187362
                                                                                                                                                   |                                                                                                                                                    
cTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTTCTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCTTTCTCGGTCTGAAACAGACCGTTATTTTTGCGGAcc                                                                                                                                                    >  1:657633/1‑150 (MQ=255)
cTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAcc                                                                                                                                                    <  1:995523/150‑1 (MQ=255)
 tGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGCCCGTTATTTTTGCGGACCa                                                                                                                                                   >  2:914415/1‑150 (MQ=255)
    tGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGAt                                                                                                                                                <  2:999931/150‑1 (MQ=255)
        ggTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAATAGCAATCTCCAAAATGCTGAGCACCTGCTACCGCATTCTCGGTCTGAAACCGGCCGTTTTTTTTGCGGGCCAGCTCATg                                                                                                                                            >  1:402565/2‑150 (MQ=255)
          tGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTCCTCCATCGGCAACCAGGCGCTGGGTAAAAAAGCAACCTCAAAAATACTGAACACCTGCTACCGCATACTCGGTCTGAACCCGACCGTAATTTTTGCGGACCGGAGCATGTa                                                                                                                                          >  1:788406/1‑150 (MQ=255)
          tGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCATAATGCTGAACACCTGCTACCGCTTTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTa                                                                                                                                          <  1:1264420/150‑1 (MQ=255)
          tGCCATTCAGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAAACTCCAAAATGCTGAACACCTCCTACCGCATTCTCGGACTGAACCCGACCGTGAGTTTTGCGGACCAGATCAAGTa                                                                                                                                          >  1:792243/1‑150 (MQ=255)
           gCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTac                                                                                                                                         <  1:794403/150‑1 (MQ=255)
            ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAATCCGACCGGTATTTTTGCGGACCAGATCATGTaca                                                                                                                                        >  1:832818/1‑150 (MQ=255)
            ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCAt                                                                                                                                             >  2:1301178/1‑145 (MQ=255)
            ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCAt                                                                                                                                             <  1:1301178/145‑1 (MQ=255)
            ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTaca                                                                                                                                        >  1:617593/1‑150 (MQ=255)
               ttCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTTAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCg                                                                                                                                     <  1:713733/150‑1 (MQ=255)
               ttCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCg                                                                                                                                     <  2:1207402/150‑1 (MQ=255)
                tCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCgg                                                                                                                                    <  1:545438/150‑1 (MQ=255)
                 cTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGc                                                                                                                                   >  1:1247069/1‑150 (MQ=255)
                 cTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGc                                                                                                                                   >  1:439371/1‑150 (MQ=255)
                        tgatTGTACCGAAATGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGGTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCt                                                                                                                            <  1:668596/150‑1 (MQ=255)
                        tgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCt                                                                                                                            <  2:313131/150‑1 (MQ=255)
                          tttGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTTGTTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTAt                                                                                                                          <  1:229600/149‑1 (MQ=255)
                            tGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGc                                                                                                                        <  1:1087257/150‑1 (MQ=255)
                             gTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCa                                                                                                                       >  1:21852/1‑150 (MQ=255)
                              tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAg                                                                                                                      >  1:243457/1‑150 (MQ=255)
                              tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAg                                                                                                                      <  1:1250557/150‑1 (MQ=255)
                                ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcg                                                                                                                    <  1:1221208/150‑1 (MQ=255)
                                 cGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAgcgc                                                                                                                   <  1:1131479/150‑1 (MQ=255)
                                      ggTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCt                                                                                                              >  2:933021/1‑150 (MQ=255)
                                       gTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTAATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTg                                                                                                             >  1:972328/1‑150 (MQ=255)
                                          tGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTg                                                                                                          <  1:1087671/150‑1 (MQ=255)
                                            ccTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCa                                                                                                        >  2:884286/1‑150 (MQ=255)
                                             cTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAACCCGCCGGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCAt                                                                                                       >  1:630658/1‑150 (MQ=255)
                                                  tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGCCCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                                  >  1:168922/1‑150 (MQ=255)
                                                  tCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATTATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                                  <  2:885225/150‑1 (MQ=255)
                                                     aTCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGt                                                                                               <  1:885021/150‑1 (MQ=255)
                                                      tCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTa                                                                                              >  2:463098/1‑150 (MQ=255)
                                                       cGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTAt                                                                                             >  2:1133290/1‑150 (MQ=255)
                                                       cGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAAAACCTTCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATAATGTACACCGGCTTCGCCTAAGCAGCGCGTTCTGGTTCATCTGTTGTTaa                                                                                             >  2:221029/1‑149 (MQ=255)
                                                        gTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGTTCATGTACACCGGCTTCGCCTAGGCAGCGCGTTCTGGTGGATCTGTTGGTATc                                                                                            >  2:822958/1‑150 (MQ=255)
                                                        gTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATc                                                                                            >  1:1142976/1‑150 (MQ=255)
                                                            aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                        <  2:640011/150‑1 (MQ=255)
                                                            aCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                        >  2:779072/1‑150 (MQ=255)
                                                                  cgcTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATgg                                                                                  <  1:775261/150‑1 (MQ=255)
                                                                    cTGGGTAAAAACGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGCCAGGGTc                                                                                >  1:547074/1‑150 (MQ=255)
                                                                        gTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCCGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATcc                                                                            >  1:1229819/1‑150 (MQ=255)
                                                                        gTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATcc                                                                            >  1:111356/1‑150 (MQ=255)
                                                                         tAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATccc                                                                           <  1:886229/150‑1 (MQ=255)
                                                                         tAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATccc                                                                           <  1:374150/150‑1 (MQ=255)
                                                                          aaaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACCATGTCATcgat                                                                          <  1:1301140/150‑4 (MQ=255)
                                                                           aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                         <  2:1178113/150‑1 (MQ=255)
                                                                           aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGGGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                         >  1:788935/1‑150 (MQ=255)
                                                                             aaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgag                                                                       <  2:796920/150‑1 (MQ=255)
                                                                              aaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGgaga                                                                      <  1:557398/150‑1 (MQ=255)
                                                                                gCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCTGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaag                                                                    <  2:524657/150‑1 (MQ=255)
                                                                                ccAATCTCCAAAAAGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaag                                                                    <  1:313364/149‑1 (MQ=255)
                                                                                     ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                               <  1:418831/150‑1 (MQ=255)
                                                                                     ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                               >  1:491144/1‑150 (MQ=255)
                                                                                       ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGACACACg                                                             >  2:1129209/1‑150 (MQ=255)
                                                                                        cAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGa                                                            <  1:1114928/150‑1 (MQ=255)
                                                                                           aaTGCTGAACACCTGCTACCTCTTTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAat                                                         <  1:1072046/150‑1 (MQ=255)
                                                                                            aTGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCAGCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatc                                                        >  1:329874/1‑150 (MQ=255)
                                                                                             tGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatca                                                       <  2:548281/150‑1 (MQ=255)
                                                                                              gCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcat                                                      >  1:888424/1‑150 (MQ=255)
                                                                                                tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtct                                                    >  1:379592/1‑150 (MQ=255)
                                                                                                tGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtct                                                    <  2:965066/150‑1 (MQ=255)
                                                                                                 gAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTGTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCAtctc                                                   >  1:1153985/1‑150 (MQ=255)
                                                                                                      ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGTAGAAACACGAAATCATCTCCGAgg                                              >  1:256428/1‑150 (MQ=255)
                                                                                                        tGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                            <  2:290459/150‑1 (MQ=255)
                                                                                                         gCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcag                                           >  2:918537/1‑150 (MQ=255)
                                                                                                             ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagc                                       >  2:768569/1‑150 (MQ=255)
                                                                                                               gCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAAAACGAAATCATCTCCGAGgcagacgcag                                     >  1:790349/1‑150 (MQ=255)
                                                                                                               gCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCACCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcag                                     >  1:1278898/1‑150 (MQ=255)
                                                                                                                      cGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGc                              >  2:1078453/1‑150 (MQ=255)
                                                                                                                        gTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTg                            <  1:421669/150‑1 (MQ=255)
                                                                                                                            gAAACCGACCGTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                        >  2:1200184/1‑150 (MQ=255)
                                                                                                                            gAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                        >  2:273444/1‑150 (MQ=255)
                                                                                                                            gAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                        >  1:219744/1‑150 (MQ=255)
                                                                                                                              aaCCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTc                      <  2:1159367/150‑1 (MQ=255)
                                                                                                                                ccgaccgTTATTTTTGCGGACCAGATCATGTACCCCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACTCGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg                    <  1:172423/150‑1 (MQ=255)
                                                                                                                                  gaccgTTATTTTTGCGGGCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGa                  <  2:972328/150‑1 (MQ=255)
                                                                                                                                   accgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAg                 <  1:1302218/150‑1 (MQ=255)
                                                                                                                                    ccgTTATTTTTGCGGACCAGATCATGTACACCTGGTTCGGCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAGGACATGGTCATCCCGGAGAAGAAACACGAATTCATCTCCGAGGCAGAAGCAGAAGTTGCTGAACTTCAGGAGc                <  1:216175/150‑1 (MQ=255)
                                                                                                                                    ccgGTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGc                <  1:112589/150‑1 (MQ=255)
                                                                                                                                     cgTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCa               >  1:563914/1‑150 (MQ=255)
                                                                                                                                       ttATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGACATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGt             >  1:824706/1‑150 (MQ=255)
                                                                                                                                        tATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGtt            >  2:1080046/1‑150 (MQ=255)
                                                                                                                                            tttGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAg        >  2:1085046/1‑150 (MQ=255)
                                                                                                                                                 ggACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCAAgtctgg   >  2:522548/1‑150 (MQ=255)
                                                                                                                                                  gACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACCCGAAATCATCGCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctggt  >  1:141492/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGT  >  NC_000913/4187067‑4187362

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: