Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 951,219 T→A 5.6% intergenic (‑139/+53) pflA ← / ← pflB pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase/formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913951,2190TA 5.6% 315.4 / 4.6 110intergenic (‑139/+53)pflA/pflBpyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase/formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
Reads supporting (aligned to +/- strand):  ref base T (53/50);  new base A (3/4);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTC  >  NC_000913/951075‑951368
                                                                                                                                                |                                                                                                                                                     
tGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTg                                                                                                                                                  <  2:894992/150‑1 (MQ=255)
  aCATTGCGGTGTTTCTCCTGATGGGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCGACGTCGAGTCTGTTTTGGCAGTCTCCTTAAAGTATAGATAGCTGACAAAATAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGa                                                                                                                                                <  1:1110097/150‑1 (MQ=255)
     ttGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGcc                                                                                                                                             <  1:1462158/150‑1 (MQ=255)
      ttctgtgtTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAATGCCCCACTTTCGTGGAGCCt                                                                                                                                            <  2:1535169/146‑1 (MQ=255)
       gCGGTGGTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCtt                                                                                                                                           <  2:468773/150‑1 (MQ=255)
          gtggTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGGGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTAt                                                                                                                                        <  2:1601861/150‑1 (MQ=255)
           tgtTTCTCCAGATGTGGCCCATCTGAGGCAGTGTTGGTGCGCAGCTCGAAGTCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTGGTGGAGCCTTTAtt                                                                                                                                       >  2:137487/1‑150 (MQ=255)
           ggtTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTAtt                                                                                                                                       <  2:176118/149‑1 (MQ=255)
            gtTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTAATg                                                                                                                                      <  2:802911/150‑1 (MQ=255)
            gtTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGGCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTg                                                                                                                                      <  1:366353/150‑1 (MQ=255)
             tttCTCCCGATGGGGCCCATCTGAGGCCGTGTTGGTGCGCAGCGCGAAGGCTACGGCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCTCCACGAACGGGGGGCCttttttag                                                                                                                                     <  2:831080/150‑3 (MQ=255)
             tttCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAACGCTCCACGAAAGTGGGGCCttttttag                                                                                                                                     <  1:263931/150‑3 (MQ=255)
              ttCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCTCCACGAAAGTGGGGCCttttttagc                                                                                                                                    <  1:193235/150‑4 (MQ=255)
              ggcctcAGATGGGCCACATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTa                                                                                                                                    >  1:1692002/6‑150 (MQ=255)
               tctcCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCTCCACGAAAGTGGGGCCttttttagcg                                                                                                                                   <  2:1465550/150‑5 (MQ=255)
               tctcCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTAc                                                                                                                                   <  2:1177373/150‑1 (MQ=255)
                ctcCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACg                                                                                                                                  >  2:1491669/1‑150 (MQ=255)
                ctcCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCGCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTCCg                                                                                                                                  >  1:216337/1‑150 (MQ=255)
                   cAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCtt                                                                                                                               >  2:811076/1‑150 (MQ=255)
                      aTGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCttttt                                                                                                                            <  2:35072/150‑1 (MQ=255)
                       tgtgGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTa                                                                                                                           >  1:1432514/1‑150 (MQ=255)
                          ggCCCATCTGAGGCCGTGTTGGTGCGCATCTCGAAGGATACGTCGAGTCAGTTTTTGCAGTCACCTTAAAATATAGATAGATGACAAAAAAGGATCTCGCGCAAAAAAAGGCCACAGTTACGGGCAGCCTATATTGTACGCATTATACtt                                                                                                                        >  2:1695861/1‑149 (MQ=255)
                           gCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGt                                                                                                                       >  1:995114/1‑150 (MQ=255)
                           gCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGt                                                                                                                       <  1:1137073/150‑1 (MQ=255)
                           gCCCATCTGAGGCCGTGTTGGTGCGCAGCGCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGt                                                                                                                       <  1:1131297/150‑1 (MQ=255)
                             ccATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTCTAGATAGCTGTCAAAAATGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTAc                                                                                                                     <  1:80220/150‑1 (MQ=255)
                                 cTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGAtt                                                                                                                 >  2:548979/1‑150 (MQ=255)
                                    cggCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCa                                                                                                              <  2:1235154/149‑1 (MQ=255)
                                    aGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCa                                                                                                              >  2:1407887/1‑150 (MQ=255)
                                    aGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCa                                                                                                              >  2:1476511/1‑150 (MQ=255)
                                     ggCCGTGTTGGTGCGCAGCGCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAACAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAg                                                                                                             <  1:337955/150‑1 (MQ=255)
                                      gCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAATGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGt                                                                                                            >  2:1137104/1‑150 (MQ=255)
                                         gtgtTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCaa                                                                                                         >  1:1148698/1‑150 (MQ=255)
                                         gtgtTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCaa                                                                                                         >  2:1527549/1‑150 (MQ=255)
                                             tGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACTCTTTTTACTGTACGATTTCAGTCAAATCt                                                                                                     <  1:709857/150‑1 (MQ=255)
                                             tGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCt                                                                                                     >  2:1251811/1‑150 (MQ=255)
                                             tGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCt                                                                                                     <  2:537915/150‑1 (MQ=255)
                                               gTGCGCAGCGCGAAGGCTACGTCGAGTCTGTTATGGCAGTCACCTTAAAGTATAGATAGCTGCCAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTaa                                                                                                   <  2:1043194/150‑1 (MQ=255)
                                                   gcAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTAc                                                                                               >  1:1500189/1‑150 (MQ=255)
                                                   gcAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTAc                                                                                               <  1:199928/150‑1 (MQ=255)
                                                   gcAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTAc                                                                                               >  2:1608514/1‑150 (MQ=255)
                                                   gcAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTAc                                                                                               >  1:1194348/1‑150 (MQ=255)
                                                   gcAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCACGCGCTAAAAAAGGCCCCACTATCGTGGAGCCTTTATTGTAAGCTTTTTACTGAACGAATACAGTCAAATCAAATTAc                                                                                               >  2:1338956/1‑150 (MQ=255)
                                                      gCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATa                                                                                            <  2:49307/150‑1 (MQ=255)
                                                       cTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAg                                                                                           >  2:990426/1‑150 (MQ=255)
                                                         cGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGAt                                                                                         >  2:1486000/1‑150 (MQ=255)
                                                           aaGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAATGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTg                                                                                       <  1:912343/150‑1 (MQ=255)
                                                           aaGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTg                                                                                       <  1:1256449/150‑1 (MQ=255)
                                                            aGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGa                                                                                      >  2:1261908/1‑150 (MQ=255)
                                                             ggCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAg                                                                                     <  2:989880/150‑1 (MQ=255)
                                                             ggCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAg                                                                                     >  2:758565/1‑150 (MQ=255)
                                                              gCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGt                                                                                    >  1:631772/1‑150 (MQ=255)
                                                              gCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGt                                                                                    <  1:687380/150‑1 (MQ=255)
                                                              gCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACATTAGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTAAGACAAATCTAATTACAAAGATTGAGt                                                                                    >  1:1563158/1‑150 (MQ=255)
                                                               cTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCTCTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTg                                                                                   <  2:805323/150‑1 (MQ=255)
                                                               cTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTg                                                                                   >  2:636951/1‑150 (MQ=255)
                                                                  cGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAg                                                                                >  1:1279907/1‑150 (MQ=255)
                                                                  cGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTAAAATCTAATTACATAGATTGAGTGAAg                                                                                <  2:552678/150‑1 (MQ=255)
                                                                     cGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTa                                                                             >  2:743483/1‑150 (MQ=255)
                                                                         tCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAg                                                                         <  1:1208935/150‑1 (MQ=255)
                                                                          cTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGt                                                                        >  1:396252/1‑150 (MQ=255)
                                                                             ttttGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGtaat                                                                     >  2:690457/1‑150 (MQ=255)
                                                                                    aGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTcc                                                              >  2:1658451/1‑150 (MQ=255)
                                                                                    aGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTcc                                                              <  1:664817/150‑1 (MQ=255)
                                                                                    aGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTcc                                                              <  2:193265/150‑1 (MQ=255)
                                                                                       cACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgc                                                           >  2:670683/1‑150 (MQ=255)
                                                                                        aCCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgct                                                          <  1:609305/150‑1 (MQ=255)
                                                                                         ccTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                         >  2:1723126/1‑150 (MQ=255)
                                                                                         ccTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                         >  2:599882/1‑150 (MQ=255)
                                                                                          cTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgc                                                        >  2:1033525/1‑150 (MQ=255)
                                                                                          cTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgc                                                        <  1:1036234/150‑1 (MQ=255)
                                                                                                gtatttatagCTGACAAAAAAGGCTCTCGCGCTAAATAAGGCCCCACTTTCTTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATGTAATTACATAGATTGAGTGAATGTACGAGTAATAACGTCCTGCTGCTGTTCt                                                  <  1:1509281/144‑1 (MQ=255)
                                                                                                gTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCt                                                  <  2:300857/150‑1 (MQ=255)
                                                                                                gTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCt                                                  <  2:1195936/150‑1 (MQ=255)
                                                                                                gTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTATCGTGGAGCCATTATTGTACGCTTATTACTGTACCATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACCTCCTGCTGCTGTTCt                                                  >  2:725349/1‑150 (MQ=255)
                                                                                                  atagatagCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCttt                                                <  2:447861/150‑1 (MQ=255)
                                                                                                    agatagCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATTCATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCttttg                                              >  2:372643/1‑148 (MQ=255)
                                                                                                    agatagCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAg                                              <  1:1428295/150‑1 (MQ=255)
                                                                                                    agatagCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGAATGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAg                                              <  2:701198/150‑1 (MQ=255)
                                                                                                    agatagCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGATGTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAg                                              >  1:523109/1‑150 (MQ=255)
                                                                                                           tGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGa                                       <  1:584474/150‑1 (MQ=255)
                                                                                                            gACAAAAAAGGCTCTCGCGCTTAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAg                                      >  1:90213/1‑150 (MQ=255)
                                                                                                             aCAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                     <  1:670683/150‑1 (MQ=255)
                                                                                                                aaaaaGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCGTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGa                                  >  2:1213703/1‑150 (MQ=255)
                                                                                                                 aaaaGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaa                                 >  2:1151566/1‑150 (MQ=255)
                                                                                                                      gCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                            <  1:803278/150‑1 (MQ=255)
                                                                                                                       ctctCGCGCTAAAAAAGGCCACACTTTCGTTGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTAc                           <  2:520586/150‑1 (MQ=255)
                                                                                                                         ctcGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTg                         <  2:1432514/150‑1 (MQ=255)
                                                                                                                             cgcTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTa                     >  2:222406/1‑150 (MQ=255)
                                                                                                                             cgcTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTa                     >  2:227458/1‑150 (MQ=255)
                                                                                                                             cgcTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAAAGTACTGCGTa                     >  2:181883/1‑150 (MQ=255)
                                                                                                                                  aaaaaGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAg                <  2:758962/150‑1 (MQ=255)
                                                                                                                                  aaaaaGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTATTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACCTACTGCGTATCCAg                >  1:879169/1‑150 (MQ=255)
                                                                                                                                     aaGGCCCCACGTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCATATCTAATTACATAGATTGAGTGAAGGTACGAGGAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCTTAGCCAGATa             >  1:389665/1‑150 (MQ=255)
                                                                                                                                         ccccACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg         >  2:1579235/1‑150 (MQ=255)
                                                                                                                                         ccccACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg         >  2:569139/1‑150 (MQ=255)
                                                                                                                                         ccccACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTAAGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg         >  2:129635/1‑150 (MQ=255)
                                                                                                                                             aCTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATg     <  2:90213/150‑1 (MQ=255)
                                                                                                                                               tttCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt   <  1:758565/150‑1 (MQ=255)
                                                                                                                                                ttCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGGGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTc  <  1:1527549/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGATAGCTGACAAAAAAGGCTCTCGCGCTAAAAAAGGCCCCACTTTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTC  >  NC_000913/951075‑951368

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: