Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,722 T→A 9.0% intergenic (‑102/+3) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,7220TA 9.0% 703.2 / 3.2 242intergenic (‑102/+3)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base T (75/136);  new base A (29/2);  total (104/138)
Fisher's exact test for biased strand distribution p-value = 3.56e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

TGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAA  >  NC_000913/3491574‑3491869
                                                                                                                                                    |                                                                                                                                                    
  aaGCTCCCCGCGCGGGGACTTGGGCTTAGCGGTGAGGTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACACGTATAAGCGGCAATTCGGGGTTAGGCGGGGTGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                                   <  1:1147813/150‑1 (MQ=255)
  aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                                   >  2:1235283/1‑150 (MQ=255)
  aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                                   <  2:2368363/150‑1 (MQ=255)
        ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGAATGCCCCCGCGCGAAAACATTACGGc                                                                                                                                             >  1:2830053/1‑150 (MQ=255)
         cAGAGACGGGACTTTTGCTTATCTGTGAGTTTCTTCACTTGATGCCTCCCTCAACGTGTCAATTTCCCAAAAGTATAAGCGGAAAATCCGGGTTACGCCGGGAGAGCACAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCa                                                                                                                                            <  1:2552420/150‑1 (MQ=255)
         cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGATTGCGCCCCCGCGAAATCAATACAGCa                                                                                                                                            >  1:3455310/1‑150 (MQ=255)
         cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCa                                                                                                                                            <  2:2578464/150‑1 (MQ=255)
          agtgaggggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGTTTGCCCGCGCGAGAAATCATTACGGCAg                                                                                                                                           <  1:2609697/144‑1 (MQ=255)
          agagACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCACTTTGCCCGCGCGAGAAATCAGTaacgccg                                                                                                                                           >  1:4626845/1‑144 (MQ=255)
          agagACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCCGGCAACATTGCCCGCCCGAGAAAAAAATACGGCAg                                                                                                                                           >  1:1738899/1‑150 (MQ=255)
               cGGGACTTTTGCTTATCGCTGAGTTTCTTCACTTGATGCCTCCCAGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACtt                                                                                                                                      <  2:2787077/150‑1 (MQ=255)
                gggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGGGATATTGGCCGCGCGGGAAAACAATAAGGCAgggcgtt                                                                                                                                     >  1:999518/1‑145 (MQ=255)
                   aCTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                                                  <  2:2965424/150‑1 (MQ=255)
                     ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAATCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGGAACCATGACCCCGCGCGAAAAAAATAAGGCAGGACGTTAGCgg                                                                                                                                >  1:4454471/1‑150 (MQ=255)
                     ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGAGTGCCCCCGCGCGCAAACAAAACCGCACGACTTTAGCgg                                                                                                                                >  2:4139400/1‑150 (MQ=255)
                        tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTCAGCGGaaa                                                                                                                             <  2:693321/150‑1 (MQ=255)
                                ggTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCCGAGTGGGCGCGCGGGAAATAATTAAGGCCGGTCGTTATCGGAAAGGGTAAca                                                                                                                     >  1:1324928/1‑149 (MQ=255)
                                 gTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTa                                                                                                                   <  2:4738224/150‑1 (MQ=255)
                                 gTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTa                                                                                                                   <  1:2025253/150‑1 (MQ=255)
                                      tttCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTTGCCGCGCGCGAAATCATTAAGGCAGGACTTTAGCGGGAAGGATAACAAACGGt                                                                                                              >  1:3090939/1‑150 (MQ=255)
                                      tttCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGt                                                                                                              <  1:723663/150‑1 (MQ=255)
                                       ttcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGtt                                                                                                             <  1:117948/150‑1 (MQ=255)
                                        tcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGttt                                                                                                            <  2:3970378/150‑1 (MQ=255)
                                          ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCGCGCGAGAAATCAATACGGCAGGACGTTAGCGGAAAGGATAACTACCCGTTTTg                                                                                                          >  2:3764312/1‑150 (MQ=255)
                                          ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCCCCCGCGAAATCAATAACGCAGGACGTTAGCGGAAAGGATAACTACCCGTTTTg                                                                                                          >  1:4984436/1‑150 (MQ=255)
                                           tcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGCTTGCCCGCGCGCGGAATAATTTCGGCCGGGCGTTTGCGGGAACGGTAACGGCCGGgttttt                                                                                                         >  1:524766/1‑148 (MQ=255)
                                            cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                                        <  1:3894445/150‑1 (MQ=255)
                                             aCTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACg                                                                                                       <  2:4651059/150‑1 (MQ=255)
                                              cTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGGCCGATTTGCCCCGCGCGAAAAAAAAAACGCCGGTCTTTAGAGGTAAGGATAAAAACCCGGTTTGTCgg                                                                                                      >  1:938841/1‑150 (MQ=255)
                                              cTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACgg                                                                                                      <  1:1411706/150‑1 (MQ=255)
                                               ttGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGc                                                                                                     <  2:3889232/150‑1 (MQ=255)
                                                  aTGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGGCCGGTTTGCCGCGCGGGTAAACATTAAGGCAGGGCGTTTGCGGGAAGGATAACAAACGGGTTTTTCGGccca                                                                                                  >  1:2082119/1‑150 (MQ=255)
                                                    gCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGAGTGGCCGCCCGCGAAAACAATAAGGCCGGGCTTTAGCGGGAAGGATAAAGACCGGGGTTggcggccccag                                                                                                >  1:3331176/1‑145 (MQ=255)
                                                    gCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGCGAAAAAATTAAGGCAGGACTTTAGCGGAAAGGATAACAACCCGGTTTGTAGGcacaat                                                                                                >  1:3331150/1‑148 (MQ=255)
                                                     ccTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGCCCCCGCGCGAAATCAATAAGGCAGGCCGGTATCGGAAAGGAGAACAACCCGGTTTGTCGGCACAGTc                                                                                               >  2:908021/1‑150 (MQ=255)
                                                      cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACGTTAGCGGAAAGGATAACTACCGGGTTTGACGGCACATTCt                                                                                              >  2:3752884/1‑150 (MQ=255)
                                                       tCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTg                                                                                             <  2:4546429/150‑1 (MQ=255)
                                                       tCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTg                                                                                             <  2:3283860/150‑1 (MQ=255)
                                                            gAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTAcn                                                                                        <  1:2698517/150‑2 (MQ=255)
                                                               cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGTAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                                     <  1:3661760/150‑1 (MQ=255)
                                                                   tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAAc                                                                                 <  1:1674902/150‑1 (MQ=255)
                                                                   tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTACCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAAc                                                                                 >  1:579648/1‑150 (MQ=255)
                                                                   tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAAATTGCCCGCCCGAGAAAAAAATACGGCAGGACATTATCGGAAAGGAAAAAAAACCGTTTTTACGGCACAATCTGGTAAGGACCCAc                                                                                 >  2:4853427/1‑150 (MQ=255)
                                                                     cATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCACCCATTCCCGCGCGCGAAATCAATAAGGCCCGACATTTTCCGAAACGGTAACAAACGGGTTTGTCGGCCCAGTATGGGACCGACCAACgg                                                                               >  2:1830183/1‑149 (MQ=255)
                                                                       tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                                             <  2:5042394/150‑1 (MQ=255)
                                                                          gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                                          <  2:3830858/150‑1 (MQ=255)
                                                                           ccACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGGCCCCGCGAGAAATAATTAAGGCAGGACGTTAGCGGAAAGGCTAACAAACGGTTTTGTCGGCACAATATGGTACGGACCAACCTCAagag                                                                         >  2:1538452/1‑150 (MQ=255)
                                                                             acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTc                                                                       <  1:2912213/150‑1 (MQ=255)
                                                                             acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTc                                                                       <  1:1498155/150‑1 (MQ=255)
                                                                             acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTc                                                                       >  2:2888632/1‑150 (MQ=255)
                                                                               aaGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCgg                                                                     <  1:5066362/150‑1 (MQ=255)
                                                                                       gCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAAACAATAACGCAGGACGTTATCGGAAAGGAAAACTACCCGTTTTGACGGCACAATATGGTACGGACCAAAGACATGAGTAGGCACCTggg                                                             >  2:4149666/1‑150 (MQ=255)
                                                                                        cGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGa                                                            <  1:1015058/150‑1 (MQ=255)
                                                                                         ggCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaa                                                           <  1:4261005/150‑1 (MQ=255)
                                                                                         ggCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaa                                                           <  1:235362/150‑1 (MQ=255)
                                                                                         ggCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAACAATACCGCAGGACGTTAGCGGAAAGGATAACAACCCGTTTTGACGGCACAATCTGGTACCGACCAACCTCAAGAGTCCGCACCAGGGaa                                                           >  2:2797851/1‑150 (MQ=255)
                                                                                          gCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaaa                                                          <  1:4158006/150‑1 (MQ=255)
                                                                                           cAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAAt                                                         <  2:3439514/150‑1 (MQ=255)
                                                                                           cAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAAt                                                         <  1:2193290/150‑1 (MQ=255)
                                                                                            aaaTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCATGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATc                                                        <  1:908403/150‑1 (MQ=255)
                                                                                              aTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCtt                                                      <  2:2696931/150‑1 (MQ=255)
                                                                                                 cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                                   <  1:2844422/150‑1 (MQ=255)
                                                                                                 cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                                   <  1:3256059/150‑1 (MQ=255)
                                                                                                 cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                                   <  1:4666454/150‑1 (MQ=255)
                                                                                                 cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                                   <  2:1738899/150‑1 (MQ=255)
                                                                                                  gggTTACGCCGGGAGAGCAGAAGCGCGGGCCCGGTTGCCCCCCCCCGAAATCAATAACGCACGACGTTAGCGGGAAGGATAACAAACCGGTTTTTCGGCCCAATCAGGTACGGCCCCACCGCAAGAGGCGGCCCCCGGGCAATCatatag                                                  >  2:624528/1‑148 (MQ=255)
                                                                                                          ccGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTc                                          <  2:3021360/150‑1 (MQ=255)
                                                                                                             ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                                       <  2:1305196/150‑1 (MQ=255)
                                                                                                             cgagaGCACAACCGCGCCCAACATTGCCCCCGCGAGAAAACATTACGGCACGACAATACCGCAAAGGATAACTACCCGAAATGACCGCACATTCTGGTACGGACCAACGTCAAGAGTCGGCACCTGGGAAATCTAATCCCCAACCACCAt                                       <  2:2204306/149‑1 (MQ=255)
                                                                                                               agagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                                     <  2:3688784/150‑1 (MQ=255)
                                                                                                                   cAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt                                 <  1:3515313/150‑1 (MQ=255)
                                                                                                                    aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                                >  2:1615752/1‑150 (MQ=255)
                                                                                                                    aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                                <  1:2446530/150‑1 (MQ=255)
                                                                                                                     gAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca                               <  1:4015022/150‑1 (MQ=255)
                                                                                                                       aGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc                             <  2:2460966/150‑1 (MQ=255)
                                                                                                                       aGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc                             <  2:209058/150‑1 (MQ=255)
                                                                                                                        gcgcGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacca                            >  2:396253/1‑150 (MQ=255)
                                                                                                                        gcgcGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacca                            <  2:2452171/150‑1 (MQ=255)
                                                                                                                            gggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTAGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt                        <  2:2032934/150‑1 (MQ=255)
                                                                                                                             ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                       <  2:3124245/150‑1 (MQ=255)
                                                                                                                             ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                       <  2:759273/150‑1 (MQ=255)
                                                                                                                               cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc                     <  1:358664/150‑1 (MQ=255)
                                                                                                                                aaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCa                    >  2:1017516/1‑150 (MQ=255)
                                                                                                                                aaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCa                    >  1:3825668/1‑150 (MQ=255)
                                                                                                                                   aTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAt                 >  1:1466415/1‑150 (MQ=255)
                                                                                                                                     tGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATAc               >  1:241687/1‑150 (MQ=255)
                                                                                                                                      gCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGTCGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATAcc              <  1:109761/150‑1 (MQ=255)
                                                                                                                                        ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc            <  2:4724784/150‑1 (MQ=255)
                                                                                                                                        ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc            <  2:2072085/150‑1 (MQ=255)
                                                                                                                                         cgcgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgcg           >  2:565900/1‑150 (MQ=255)
                                                                                                                                          gcgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgcga          >  2:941404/1‑149 (MQ=255)
                                                                                                                                          gcgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGt          >  2:1909245/1‑150 (MQ=255)
                                                                                                                                           cgcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTa         >  2:4345758/1‑150 (MQ=255)
                                                                                                                                             cgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAAc       >  2:3557396/1‑150 (MQ=255)
                                                                                                                                             cgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAAc       >  2:1051008/1‑150 (MQ=255)
                                                                                                                                                gaAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGa    >  2:1162131/1‑150 (MQ=255)
                                                                                                                                                gaAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGa    <  1:1905619/150‑1 (MQ=255)
                                                                                                                                                 aaaTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGaa   <  2:4831355/150‑1 (MQ=255)
                                                                                                                                                 aaaTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCATAACCGaa   <  2:1402265/150‑1 (MQ=255)
                                                                                                                                                  aaTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGaaa  <  2:3175187/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAA  >  NC_000913/3491574‑3491869

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