Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,881,358 T→A 5.6% D151V (GAC→GTC)  dmlR ← DNA‑binding transcriptional regulator DmlR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,3580TA 5.6% 278.9 / 5.3 89D151V (GAC→GTC) dmlRDNA‑binding transcriptional regulator DmlR
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base A (2/3);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

GATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAA  >  NC_000913/1881209‑1881503
                                                                                                                                                     |                                                                                                                                                 
gatgatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGt                                                                                                                                                   >  1:999130/1‑150 (MQ=255)
gatgatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGt                                                                                                                                                   >  1:645702/1‑150 (MQ=255)
gatgatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGt                                                                                                                                                   <  2:52001/150‑1 (MQ=255)
   gatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGt                                                                                                                                                   >  2:161198/1‑147 (MQ=255)
   gatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGt                                                                                                                                                   <  1:161198/147‑1 (MQ=255)
    atCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATc                                                                                                                                               >  1:800445/1‑150 (MQ=255)
     tCGCGCTCTTTAATGAGCAGGCACGGCAATGAAGAAAGATCGTTTAAATGTTTCGGCGCGAGATGTTGCTAAATAAATTCCGGAGACGCGCAGAGAATAAGGTAATTGGTCGCCAGTTTGAGGTCGATTAAATTAGGGGCGATGTCATCg                                                                                                                                              <  2:333354/150‑1 (MQ=255)
      cgcgCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGc                                                                                                                                             >  2:168879/1‑150 (MQ=255)
      cgcgCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGc                                                                                                                                             >  2:432697/1‑150 (MQ=255)
       gcgcTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGcc                                                                                                                                            <  2:645702/150‑1 (MQ=255)
          ctctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccct                                                                                                                                         <  1:1143086/150‑4 (MQ=255)
          ctctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccct                                                                                                                                         <  1:838550/150‑4 (MQ=255)
            ctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGc                                                                                                                                       >  2:1566957/1‑150 (MQ=255)
              ttAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGa                                                                                                                                     <  2:849237/150‑1 (MQ=255)
              ttAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGa                                                                                                                                     >  1:431427/1‑150 (MQ=255)
              ttAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccctaatt                                                                                                                                     <  1:1018132/150‑8 (MQ=255)
                aaTGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAat                                                                                                                                   <  1:146111/150‑1 (MQ=255)
                   gACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                                                                >  1:545974/1‑150 (MQ=255)
                    aCCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCg                                                                                                                               >  2:1290697/1‑150 (MQ=255)
                     ccAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                                                                              <  2:217/150‑1 (MQ=255)
                     ccAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                                                                              >  2:287203/1‑150 (MQ=255)
                        ggCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGAt                                                                                                                           >  1:659637/1‑150 (MQ=255)
                          cACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCg                                                                                                                         <  1:1267341/150‑1 (MQ=255)
                          cACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCg                                                                                                                         <  1:794421/150‑1 (MQ=255)
                          cACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCg                                                                                                                         <  2:1249034/150‑1 (MQ=255)
                              gCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGacac                                                                                                                     <  1:206004/150‑1 (MQ=255)
                               cAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACAcc                                                                                                                    <  1:766083/150‑1 (MQ=255)
                                  tGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTc                                                                                                                 >  2:1404144/1‑150 (MQ=255)
                                   gCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCa                                                                                                                <  1:1566957/150‑1 (MQ=255)
                                   gCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCa                                                                                                                >  2:212314/1‑150 (MQ=255)
                                    cAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAt                                                                                                               >  2:438358/1‑150 (MQ=255)
                                    cAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAt                                                                                                               <  2:1220291/150‑1 (MQ=255)
                                      caaaGATCGGTAAAATGTTTCGGCGCGCCATGTTGCGAAAAAAATTCCGGCGACGCGCAGATAATACGGTAATTGGGCGCCAGTTTGGGGGCTATTAAATTTGGGGCTTTGTCATCGCCAATGCTAATATCTATATCGACACCTTCATTg                                                                                                             <  2:424261/149‑1 (MQ=255)
                                        aaGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAc                                                                                                           >  1:74004/1‑150 (MQ=255)
                                         aGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAcc                                                                                                          <  1:357361/150‑1 (MQ=255)
                                         aGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAcc                                                                                                          >  2:1490746/1‑150 (MQ=255)
                                             cGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAAt                                                                                                      <  2:328157/150‑1 (MQ=255)
                                               gTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCg                                                                                                    <  2:1250833/150‑1 (MQ=255)
                                               gTTAAATGTTTCGGCGCACCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCg                                                                                                    <  2:1438105/150‑1 (MQ=255)
                                                   aaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCa                                                                                                >  1:467803/1‑150 (MQ=255)
                                                     tGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGt                                                                                              >  2:1444233/1‑150 (MQ=255)
                                                      gTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTc                                                                                             <  1:342234/150‑1 (MQ=255)
                                                       tttCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCg                                                                                            >  2:704731/1‑150 (MQ=255)
                                                         tCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGt                                                                                          >  1:987488/1‑150 (MQ=255)
                                                           ggCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCt                                                                                        <  1:22997/150‑1 (MQ=255)
                                                            gcgcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCtt                                                                                       >  2:601215/1‑150 (MQ=255)
                                                             cgcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTc                                                                                      >  2:1046501/1‑150 (MQ=255)
                                                               cgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCca                                                                                    >  2:1663239/1‑150 (MQ=255)
                                                                 ccATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCcaca                                                                                  >  2:463545/1‑150 (MQ=255)
                                                                        ttgtgttAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAAc                                                                           <  1:181611/144‑1 (MQ=255)
                                                                        gcgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAAc                                                                           <  1:432697/150‑1 (MQ=255)
                                                                        gcgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTGCATTGCCCAAAACGACCAGGCGGTCTACCCCATCAAAAc                                                                           >  1:772435/1‑150 (MQ=255)
                                                                         cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGac                                                                                                                       <  1:1159156/105‑1 (MQ=255)
                                                                         cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGac                                                                                                                       >  2:1159156/1‑105 (MQ=255)
                                                                         cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACg                                                                          >  1:1154015/1‑150 (MQ=255)
                                                                         cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACg                                                                          >  2:1506654/1‑150 (MQ=255)
                                                                          gcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGt                                                                         <  2:846384/150‑1 (MQ=255)
                                                                                aaTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGACTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTcc                                                                   <  1:456824/150‑1 (MQ=255)
                                                                                 aTTCCGGAGACGCGAAGAGAATACGGTAATTGGTCGACAGTTTGCGGGCGATTAAATTAGGGGCGATGTAATCGCAAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCa                                                                  >  2:583472/1‑150 (MQ=255)
                                                                                  ttCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAg                                                                 <  1:601215/150‑1 (MQ=255)
                                                                                    ccGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGCCACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGtt                                                               >  1:257021/1‑150 (MQ=255)
                                                                                      ggCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTAATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGc                                                             <  2:1540826/150‑1 (MQ=255)
                                                                                       gCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGAAAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCg                                                            <  2:1523310/150‑1 (MQ=255)
                                                                                        cGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCgg                                                           >  1:1551269/1‑150 (MQ=255)
                                                                                         gACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGa                                                          >  2:741415/1‑150 (MQ=255)
                                                                                           cgcgCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATa                                                        <  1:168879/150‑1 (MQ=255)
                                                                                             cgcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAAGTCCAGTTGAGGATAcc                                                      >  1:270180/1‑149 (MQ=255)
                                                                                              gcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGAGATGTAATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGc                                                     <  2:1556419/150‑1 (MQ=255)
                                                                                               cAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGcc                                                    <  2:74004/150‑1 (MQ=255)
                                                                                                 gagaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCtt                                                  >  1:545602/1‑150 (MQ=255)
                                                                                                   gaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCTATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCg                                                <  1:1089255/150‑1 (MQ=255)
                                                                                                     aTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGTGGATACGCCTTCGcc                                              <  1:1046501/150‑1 (MQ=255)
                                                                                                     aTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGcc                                              >  2:1653258/1‑150 (MQ=255)
                                                                                                                   cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGc                                >  2:902350/1‑150 (MQ=255)
                                                                                                                   cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGc                                >  1:219436/1‑150 (MQ=255)
                                                                                                                   cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGc                                <  1:617325/150‑1 (MQ=255)
                                                                                                                          ttGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGcc                         <  2:551035/150‑1 (MQ=255)
                                                                                                                           tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCgg                                                           >  1:31279/1‑115 (MQ=255)
                                                                                                                           tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCgg                                                           <  2:31279/115‑1 (MQ=255)
                                                                                                                           tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGcca                        <  1:1404144/150‑1 (MQ=255)
                                                                                                                             cGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGTGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccacc                      >  1:286522/1‑150 (MQ=255)
                                                                                                                                gCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCaccacc                   <  1:704731/150‑1 (MQ=255)
                                                                                                                                gCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCaccacc                   <  1:542378/150‑1 (MQ=255)
                                                                                                                                 cGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCt                  >  2:1247398/1‑150 (MQ=255)
                                                                                                                                  gATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccacc                      >  2:1625925/1‑145 (MQ=255)
                                                                                                                                  gATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccacc                      <  1:1625925/145‑1 (MQ=255)
                                                                                                                                   aTTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgt                <  1:1490746/150‑1 (MQ=255)
                                                                                                                                     tAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtcg              >  1:185634/1‑150 (MQ=255)
                                                                                                                                       aaTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtcgtc            >  1:355910/1‑150 (MQ=255)
                                                                                                                                         ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCa          <  1:212314/150‑1 (MQ=255)
                                                                                                                                         ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCa          >  1:338419/1‑150 (MQ=255)
                                                                                                                                         gcattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATACCGCCGGCGCCCCCACCTGTCGTCCa          >  2:846413/6‑150 (MQ=255)
                                                                                                                                          cattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTAGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACAACCTGTCGTCCaa         >  1:1432411/5‑150 (MQ=255)
                                                                                                                                                 gATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCaaa  <  1:741415/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GATGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAA  >  NC_000913/1881209‑1881503

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