Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,912,751 C→A 7.0% intergenic (‑175/+17) htpX ← / ← prc putative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,7510CA 7.0% 383.5 / 13.3 137intergenic (‑175/+17)htpX/prcputative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base C (71/55);  new base A (8/2);  total (79/58)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 138 total aligned reads displayed.

AAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTT  >  NC_000913/1912602‑1912896
                                                                                                                                                     |                                                                                                                                                 
aaaaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTc                                                                                                                                                   >  2:954260/1‑150 (MQ=255)
 aaaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCt                                                                                                                                                  >  2:103/1‑150 (MQ=255)
 aaaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCt                                                                                                                                                  >  2:1583641/1‑150 (MQ=255)
   aCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTtg                                                                                                                                                <  1:1697032/150‑1 (MQ=255)
     gCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTTAAAAATCAGGCACAATTTCTtgtg                                                                                                                                              >  1:1657204/1‑150 (MQ=255)
        atatGCGATACCCACATCGTATGTGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCt                                                                                                                                           <  2:1500218/150‑1 (MQ=255)
         tatGCGATACCCAAATCGTATGTGTTACTCGACTAATTTCAAGTCTGTATAGTGCGTAAGTAAAAAGAAAGACACATCTTTACATTTTGTAGCATCTGCTTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAATTCTTGTGCCtg                                                                                                                                          <  1:289184/150‑1 (MQ=255)
           tGCTATACCCACATCGTTTGTGTTACGAGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgat                                                                                                                                        <  2:410091/150‑1 (MQ=255)
           tGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgat                                                                                                                                        >  2:935238/1‑150 (MQ=255)
            gCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgatt                                                                                                                                       <  2:1449665/150‑1 (MQ=255)
            gCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgatt                                                                                                                                       >  1:1351362/1‑150 (MQ=255)
              gATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgattga                                                                                                                                     >  1:1719719/1‑150 (MQ=255)
               aTACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCtgattgat                                                                                                                                    >  2:335608/1‑150 (MQ=255)
                tACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAATTCACCAGAATGACACAAATTTACATTTTGTAGCTTCTGATGTACGGCATCATGTCGCTGTTAAAAAATCAGGCACTCTTTCTTGAGCCTGATTGata                                                                                                                                   <  1:32643/150‑1 (MQ=255)
                tACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGata                                                                                                                                   <  2:1571683/150‑1 (MQ=255)
                   ccACATCCTATGGGTTACGCGACTATTTTCAAGTCTGGAGAGTGCGGAAGTCACCAGAAAGACACAACGTTACATTTTGTAGCATCTGAATTACGGCATCTTGTCGCTGTAAAAAAAGCCGGCACAATTTCTTGTGCCTGATTGATATTa                                                                                                                                <  1:897066/150‑1 (MQ=255)
                    cacaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTaa                                                                                                                               >  1:1388396/1‑149 (MQ=255)
                           aaTGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAATTCACCAGAAATACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACgg                                                                                                                        <  1:1622576/149‑1 (MQ=255)
                               ggTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGc                                                                                                                    >  1:1663000/1‑150 (MQ=255)
                                 ttACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCgg                                                                                                                  <  2:1173066/150‑1 (MQ=255)
                                  tACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCggg                                                                                                                 <  1:1303592/150‑1 (MQ=255)
                                  tACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCggg                                                                                                                 <  1:1310563/150‑1 (MQ=255)
                                     gcgACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTATACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgttg                                                                                                              <  2:86491/150‑1 (MQ=255)
                                      cgACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgttgt                                                                                                             >  1:321743/1‑150 (MQ=255)
                                         cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTTCCTGATTGATATTACTTGACGGGAGCGGGTTGTTcc                                                                                                          >  1:785565/1‑150 (MQ=255)
                                          tATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCcg                                                                                                         >  1:869534/1‑150 (MQ=255)
                                             tttCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCgg                                                                                                      >  1:605619/1‑150 (MQ=255)
                                             tttCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCgg                                                                                                      >  1:1243684/1‑150 (MQ=255)
                                                cAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTc                                                                                                   <  2:1411784/150‑1 (MQ=255)
                                                 aaGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCt                                                                                                  >  2:833959/1‑150 (MQ=255)
                                                  aGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTg                                                                                                 >  1:1466171/1‑150 (MQ=255)
                                                     cTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCt                                                                                              >  2:298262/1‑150 (MQ=255)
                                                      tGGATAGTGCGTAAGGCACCAGAAAGTCACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCtt                                                                                             <  1:396712/150‑1 (MQ=255)
                                                      tGGATAGTGCGGAAGTCACCAGAAAGACACACCTGTACATTGTGGAGCATCGGAGTTACGGCATCGTGGCGCTGTTAAAAAATCAGGCACAATGTCTTGTGACGTATTGATATTACTTGACGGGAGCGGGTTGGTCCGCGGGTCAGGCtg                                                                                             >  2:1056445/1‑149 (MQ=255)
                                                         aTAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCttttt                                                                                          <  1:335608/150‑1 (MQ=255)
                                                          tAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTc                                                                                         >  2:961797/1‑150 (MQ=255)
                                                           aGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCa                                                                                        <  1:1176872/150‑1 (MQ=255)
                                                             tGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAACTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAg                                                                                      <  2:125458/150‑1 (MQ=255)
                                                             tGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAg                                                                                      >  2:233632/1‑150 (MQ=255)
                                                             tGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAg                                                                                      <  2:574164/150‑1 (MQ=255)
                                                                gTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCtt                                                                                   >  2:532340/1‑150 (MQ=255)
                                                                 tAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATAACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTc                                                                                  <  2:1226934/150‑1 (MQ=255)
                                                                  aaGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTTTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCg                                                                                 <  2:1659147/150‑1 (MQ=255)
                                                                  aaGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCg                                                                                 <  1:931489/150‑1 (MQ=255)
                                                                  aaGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCg                                                                                 <  2:978072/150‑1 (MQ=255)
                                                                      cacaAGCAATACACAACTTTACATTTTCTAGCATCTGATTTCCGGGATCTTGTCGCTGTTATAAAATCTGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAg                                                                             <  2:86286/150‑1 (MQ=255)
                                                                      cACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCGGGTAAAAAATCAGGCCCACTTTCTGGGGCCGGCTGGAGATTACTTGCCGGGAGCGGGTTGTTCCGCGGGGCTGGCTTTTTCAAGCTTCTCCAg                                                                             >  2:323205/1‑150 (MQ=255)
                                                                      cACCAGAAAGACACAACTTTAACTTTTGTATCATCTGATTTACGGCATCTTGGCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAg                                                                             <  1:1119980/150‑1 (MQ=255)
                                                                        ccAGAAAGACACAACTTTACATTTTGTAGCATCTGAATTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGAt                                                                           >  1:1651857/1‑150 (MQ=255)
                                                                         cAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGCCGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATc                                                                          <  1:1583641/150‑1 (MQ=255)
                                                                          aGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCg                                                                         >  1:1325255/1‑150 (MQ=255)
                                                                           gAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGa                                                                        >  2:901479/1‑150 (MQ=255)
                                                                                 aacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  <  1:671492/149‑1 (MQ=255)
                                                                                  acaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGat                                                                 <  1:132827/150‑1 (MQ=255)
                                                                                   caACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGata                                                                >  2:1405921/1‑150 (MQ=255)
                                                                                     aCTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACATGACGGGAGAAGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt                                                              >  1:222313/1‑150 (MQ=255)
                                                                                      cTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCAAGATCGAGTGCAatcggt                                                             >  2:1475353/1‑146 (MQ=255)
                                                                                       tttACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCa                                                            >  2:783743/1‑150 (MQ=255)
                                                                                         tACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTTCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                          <  2:785565/150‑1 (MQ=255)
                                                                                         tACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                          <  1:456225/150‑1 (MQ=255)
                                                                                         tACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                          <  1:1215598/150‑1 (MQ=255)
                                                                                         tACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                          <  1:935238/150‑1 (MQ=255)
                                                                                               ttGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCa                                                    <  1:743697/150‑1 (MQ=255)
                                                                                                   aGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACTGGAGCGGGTTGTTCCGAGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACATTCTCATCCa                                                >  1:353018/1‑150 (MQ=255)
                                                                                                    gCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAg                                               >  2:1114619/1‑150 (MQ=255)
                                                                                                        cTGCTTTACGGCATCTTGTCGCTGTGAAAAAATCAGGCACAATTGCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGGTGTTCCGCTGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                           <  1:83116/150‑1 (MQ=255)
                                                                                                         tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTTCGATATTCACCGTCTCATCCAGATAAg                                          >  2:120190/1‑150 (MQ=255)
                                                                                                         tGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAg                                          <  2:1637143/150‑1 (MQ=255)
                                                                                                          gATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAgg                                         >  1:400829/1‑150 (MQ=255)
                                                                                                          gATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAgg                                         >  1:122217/1‑150 (MQ=255)
                                                                                                          gATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAgg                                         >  1:1218187/1‑150 (MQ=255)
                                                                                                              tACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATcc                                     >  2:929398/1‑150 (MQ=255)
                                                                                                               aCGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCg                                    <  2:583046/150‑1 (MQ=255)
                                                                                                                cGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCgg                                   >  2:1569180/1‑150 (MQ=255)
                                                                                                                  gCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCt                                 >  1:27082/1‑150 (MQ=255)
                                                                                                                  gCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCt                                 >  1:993040/1‑150 (MQ=255)
                                                                                                                        tGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGt                           <  1:532340/150‑1 (MQ=255)
                                                                                                                        ggTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGt                           <  2:1346898/149‑1 (MQ=255)
                                                                                                                         gTCGCTGTTAAAAAATCAGGCACAATTTGTTGTGCCTGATTGATATTACTTGACGGTAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTTCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTa                          <  2:1020758/150‑1 (MQ=255)
                                                                                                                            gCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATc                       >  1:704462/1‑150 (MQ=255)
                                                                                                                               gTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCttt                    >  1:75056/1‑150 (MQ=255)
                                                                                                                                ttAAAAAATCAGGGACATTTTCTTGTGCGTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCATTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTc                   <  2:878354/150‑1 (MQ=255)
                                                                                                                                 tAAATAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGGTGTTCCGCGGGGCTGGCTTTTTGAAGCTTCGCCAGATCGAGTGCGATATTCTCCGTCTCATCCAGATATGGATCCGGCTCCTGGTAATCTTTCg                  <  1:142011/150‑1 (MQ=255)
                                                                                                                                 tAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCg                  >  1:1459088/1‑150 (MQ=255)
                                                                                                                                   aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGt                >  1:86846/1‑150 (MQ=255)
                                                                                                                                    aaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTTCTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTa               >  1:591857/1‑150 (MQ=255)
                                                                                                                                    aaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTa               >  2:866826/1‑150 (MQ=255)
                                                                                                                                      aaTCAGTCACAAGAAATTGTGCCAGATTGATATAACTTGACGGGAGCGGTGTGTTCCGAGGGTCTGGCTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCCCCGTCTCATCCATATAAGGATCCGGCTCCTAGTAATCTTTCGGTAGaa            >  2:1589128/1‑149 (MQ=255)
                                                                                                                                       aTCAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat            >  1:1010535/1‑150 (MQ=255)
                                                                                                                                       aTCAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat            >  2:1491572/1‑150 (MQ=255)
                                                                                                                                         cAGGCAGAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTGGCTCCGCGGGTCTGGCTTTTTCAAGCTGCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTTGTAATCTTTCGGTAGatca          <  2:114893/150‑1 (MQ=255)
                                                                                                                                         cAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          >  2:448881/1‑150 (MQ=255)
                                                                                                                                         cAGGCACAATTTCTTGTGACTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGGCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          <  1:1727981/150‑1 (MQ=255)
                                                                                                                                         cAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          >  2:42651/1‑150 (MQ=255)
                                                                                                                                          aGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatcat         >  2:1536905/1‑150 (MQ=255)
                                                                                                                                          aGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatcat         >  2:128138/1‑150 (MQ=255)
                                                                                                                                            gCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATcc       <  2:1663000/150‑1 (MQ=255)
                                                                                                                                             cacaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCa      >  1:156421/1‑150 (MQ=255)
                                                                                                                                             cacaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCa      >  1:1844/1‑150 (MQ=255)
                                                                                                                                                aaTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGtt   <  2:1444025/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTT  >  NC_000913/1912602‑1912896

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: