Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 572,115 C→A 6.2% A221D (GCT→GAT)  ybcM → DLP12 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913572,1150CA 6.2% 1349.4 / 9.8 434A221D (GCT→GAT) ybcMDLP12 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (209/183);  new base A (2/39);  total (211/223)
Fisher's exact test for biased strand distribution p-value = 2.24e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 446 total aligned reads displayed.

GCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGTTTTAAGA  >  NC_000913/571966‑572264
                                                                                                                                                     |                                                                                                                                                     
gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                                                                                       >  1:3969931/1‑150 (MQ=255)
gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                                                                                       <  1:8166665/150‑1 (MQ=255)
gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                                                                                       <  1:1990556/150‑1 (MQ=255)
gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                                                                                       <  2:2866427/150‑1 (MQ=255)
 ctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGTGTTGTCTGCt                                                                                                                                                      >  2:7178017/1‑150 (MQ=255)
 ctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCt                                                                                                                                                      <  2:713131/150‑1 (MQ=255)
 ctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCt                                                                                                                                                      <  2:4600994/150‑1 (MQ=255)
  tATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCta                                                                                                                                                     <  1:7240682/150‑1 (MQ=255)
   aTCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCtat                                                                                                                                                    <  1:7070988/150‑1 (MQ=255)
       gAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACaa                                                                                                                                                >  2:7382889/1‑150 (MQ=255)
        aaCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACaaa                                                                                                                                               <  1:6520079/150‑1 (MQ=255)
         aCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAac                                                                                                                                              <  1:607851/150‑1 (MQ=255)
         aCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAac                                                                                                                                              <  2:1722873/150‑1 (MQ=255)
           ccGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAacga                                                                                                                                            >  2:1943526/1‑150 (MQ=255)
             ggCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAacgacg                                                                                                                                          >  1:6407840/1‑150 (MQ=255)
              gCACAAAGTGGACCGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACTTTACATCGCCGTCTTGCACGCGAGGGGGTAAGTTTCCAAAGTTTACTGGACGATGTGAGTTTAAATATTGCGTTGTCTGCTATACAAACGACgg                                                                                                                                         <  2:5920169/150‑1 (MQ=255)
                acaAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTa                                                                                                                                       >  2:3243274/1‑150 (MQ=255)
                acaAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACAACGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTa                                                                                                                                       >  1:1003312/1‑150 (MQ=255)
                     aTGGACAGCAAACAAGGTTGCTCTATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCATGCGTGGGGTTATGTTTCCAAAGTATACTGGACGTTGTGTGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAAcc                                                                                                                                  <  2:897242/150‑1 (MQ=255)
                      tGGACAGCAAACAAGGTTGCTCGATATCTCTACTTTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGTGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCt                                                                                                                                 <  2:7550340/150‑1 (MQ=255)
                        gACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCtat                                                                                                                               >  1:465786/1‑150 (MQ=255)
                        gACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCtat                                                                                                                               >  2:2719925/1‑150 (MQ=255)
                         aCAGCAAACAAGGTTGCTCGATATCTCTACATTTGTGTTTCTACTTTACATCGCCGTCTAGCAAGCGTGGGGGTAAGTTTCCAAAGTATACTGGAGGATGTGAGGTTAAATAATGCGTGGTCTGCTATACAAACGACGGTAAAACCtata                                                                                                                              <  2:1345553/150‑1 (MQ=255)
                          cAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATaa                                                                                                                             >  1:8448162/1‑150 (MQ=255)
                             cAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGTTTAAATAATGCTTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCg                                                                                                                          <  2:2192122/150‑1 (MQ=255)
                              aaaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCga                                                                                                                         >  2:6914897/1‑150 (MQ=255)
                               accttgtttGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCgag                                                                                                                        >  1:4102734/8‑150 (MQ=255)
                               aaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCgag                                                                                                                        >  1:9347195/1‑150 (MQ=255)
                                aCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCgaga                                                                                                                       >  1:373814/1‑150 (MQ=255)
                                  aaGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGAtt                                                                                                                     >  1:6006509/1‑150 (MQ=255)
                                  aaGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGAtt                                                                                                                     >  1:7993765/1‑150 (MQ=255)
                                  aaGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGAtt                                                                                                                     >  2:121369/1‑150 (MQ=255)
                                     gTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGcc                                                                                                                  >  2:7524713/1‑150 (MQ=255)
                                      ttGCTCTATATCTCTACATTTCTGTTTATACATTACATCGCCTTCTAGCACGCGAGGGGTTAAGTTTCCAAATTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCa                                                                                                                 <  2:4032434/150‑1 (MQ=255)
                                      ttGCTCTATAACTCTACATTTCTGTTTCTACATTACCTCGCCGTCTAGCACGGGTGGGGTTAAGTTTCCAAATTATAGTGTACGTTGTGTGTTTAATTATTGCTTTTTCTTATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCa                                                                                                                 <  1:1524425/150‑1 (MQ=255)
                                        gCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCaga                                                                                                               <  2:6266101/150‑1 (MQ=255)
                                          tCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAgaga                                                                                                             >  2:3646758/1‑150 (MQ=255)
                                               atCTCTAATTTTCTGTTTATACATTACATCGCCGTTTAGCAAGCGAGGGGTTAATTTACAAATGTATACTGGACGATTTGTGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATg                                                                                                        <  2:8642070/150‑1 (MQ=255)
                                                  tctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGtt                                                                                                     >  1:2335430/1‑150 (MQ=255)
                                                  tctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGtt                                                                                                     >  1:5391102/1‑150 (MQ=255)
                                                  tctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGtt                                                                                                     >  1:2662225/1‑150 (MQ=255)
                                                    tACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTtat                                                                                                   >  1:3018372/1‑150 (MQ=255)
                                                    tACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTtat                                                                                                   >  1:2484437/1‑150 (MQ=255)
                                                     aaaTTTCTGTTTCTACATTACATCGCCGTCTAGCATGGGTGGGTGTATTTTTCCAAAGTATACTGGACTTTGTGTGGTTAAATAAGGCGTTGTCTGATATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTtata                                                                                                  <  1:2203423/148‑1 (MQ=255)
                                                      cATTTGTTTTTATACTTTACATCGCCTTCTAGCAAGGGTGGGGTTAATTTTCAAATGTATACTGGACGTTGTGTGTTTAAATAATGGGTGGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATaa                                                                                                 <  1:1789447/150‑1 (MQ=255)
                                                         ttCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAgtg                                                                                              >  1:7754928/1‑150 (MQ=255)
                                                             gTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCt                                                                                          >  1:6484921/1‑150 (MQ=255)
                                                                tCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTtct                                                                                       >  1:5613952/1‑150 (MQ=255)
                                                                  tacattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAATGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCg                                                                                     <  2:6599924/150‑1 (MQ=255)
                                                                  tacattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACTTTGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCg                                                                                     <  2:483169/150‑1 (MQ=255)
                                                                  tacattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCg                                                                                     >  2:8009380/1‑150 (MQ=255)
                                                                   acattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATTTGAGTTTAAATATTGCGTTGTCTGATATAAAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGt                                                                                    <  1:8342569/150‑1 (MQ=255)
                                                                   acattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGt                                                                                    >  2:4719786/1‑150 (MQ=255)
                                                                    cattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGtt                                                                                   >  2:5080511/1‑150 (MQ=255)
                                                                     attacatCGCCGTCTAGCAAGCGTGGGGGTAATTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                                                                                  <  1:8215842/150‑1 (MQ=255)
                                                                     attacatCGCCGTCTAGCAAGCGAGGGGTTAAGTTTCCAAATTATACTGGACGTTGTGAGGTTAAATATTGGGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                                                                                  <  1:8750956/150‑1 (MQ=255)
                                                                       tacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTa                                                                                >  1:8909239/1‑150 (MQ=255)
                                                                          atCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTg                                                                             >  2:5324336/1‑150 (MQ=255)
                                                                           tCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGa                                                                            >  1:6639151/1‑150 (MQ=255)
                                                                           tCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGa                                                                            >  1:5810860/1‑150 (MQ=255)
                                                                            cGCCGTCTTGCATGCGAGGGTGTATGTTTCCAAAGTATACTGGACGATGTTAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGaa                                                                           <  2:8137376/150‑1 (MQ=255)
                                                                              ccGTCTAGCAAGCGAGGGGGTAAGTTTCCAATGTATACTGGACGATGTTAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAg                                                                         <  1:1314162/150‑1 (MQ=255)
                                                                                 tCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGAtt                                                                      >  1:4444655/1‑150 (MQ=255)
                                                                                   tAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTc                                                                    >  2:7751912/1‑150 (MQ=255)
                                                                                    aGCAAGCGAGGGGTTAAGTTTCAAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCa                                                                   <  2:7651527/150‑1 (MQ=255)
                                                                                    aGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCa                                                                   <  2:6204882/150‑1 (MQ=255)
                                                                                    aGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCa                                                                   >  1:359491/1‑150 (MQ=255)
                                                                                              ggggTAATTTTCAAAAGTATACTGGACGATGTGAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGtttt                                                         <  1:7892019/150‑1 (MQ=255)
                                                                                               gggTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTa                                                        >  1:1641248/1‑150 (MQ=255)
                                                                                                ggtatttttCCAAAGTATACGGGACGATTTTTGGTTAAATATTGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTaa                                                       <  2:5512359/144‑1 (MQ=255)
                                                                                                ggTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTaa                                                       <  1:2921649/150‑1 (MQ=255)
                                                                                                ggTAAGTTTCCAAAGTATACTGGACGAGGTGAGGTTAAATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTaa                                                       <  1:8527756/150‑1 (MQ=255)
                                                                                                 gTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAat                                                      >  2:4155514/1‑150 (MQ=255)
                                                                                                 gTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAat                                                      <  1:7418094/150‑1 (MQ=255)
                                                                                                   aaGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAatat                                                    <  1:9184079/150‑1 (MQ=255)
                                                                                                   aaGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAatat                                                    >  2:4269967/1‑150 (MQ=255)
                                                                                                       ttCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAaca                                                <  1:6997601/150‑1 (MQ=255)
                                                                                                          cAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCa                                             >  1:1139725/1‑150 (MQ=255)
                                                                                                            aaGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAag                                           >  2:7894984/1‑150 (MQ=255)
                                                                                                               tataCTGGACGTTGTGAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAagaga                                        <  1:4223613/150‑1 (MQ=255)
                                                                                                               tataCTGGACGATGTTAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAagaga                                        <  2:9189949/150‑1 (MQ=255)
                                                                                                                    tgtactttgtgAGTTTAAATATTGCGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATaa                                   <  1:8125753/143‑1 (MQ=255)
                                                                                                                      gACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGa                                 >  2:8805671/1‑150 (MQ=255)
                                                                                                                            gtgAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCt                           <  2:703379/150‑1 (MQ=255)
                                                                                                                              taGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTc                         <  2:6072353/149‑1 (MQ=255)
                                                                                                                                ggTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCa                       <  2:2335430/150‑1 (MQ=255)
                                                                                                                                 gTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTTCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCag                      >  2:5210649/1‑150 (MQ=255)
                                                                                                                                  ttAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCaga                     >  1:6150434/1‑150 (MQ=255)
                                                                                                                                  ttAAATAATGCGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCaga                     <  2:1672218/150‑1 (MQ=255)
                                                                                                                                   tAAATATTGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCagag                    <  2:2484437/150‑1 (MQ=255)
                                                                                                                                   tAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCagag                    <  2:3018372/150‑1 (MQ=255)
                                                                                                                                    aaaTAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCagaga                   <  2:1132304/150‑1 (MQ=255)
                                                                                                                                    aaaTAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCagaga                   <  2:8435821/150‑1 (MQ=255)
                                                                                                                                     aataatGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCagagag                  <  2:647938/150‑1 (MQ=255)
                                                                                                                                             gTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAgtg          >  1:4922106/1‑150 (MQ=255)
                                                                                                                                                gTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGttt       <  2:3694853/150‑1 (MQ=255)
                                                                                                                                                   tGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGTTTTaa    >  1:3232681/1‑150 (MQ=255)
                                                                                                                                                    gCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGTTTTaag   >  2:1459115/1‑150 (MQ=255)
                                                                                                                                                    gCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGTTTTaag   <  2:5761236/150‑1 (MQ=255)
                                                                                                                                                    gCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGTTTTaag   >  1:3369064/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAAGCTTCCAGAGAGTAAAAGTGTTTTAAGA  >  NC_000913/571966‑572264

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