Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 572,088 A→T 5.0% D212V (GAT→GTT)  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,0880AT 5.0% 1671.9 / 10.1 527D212V (GAT→GTT) ybcMDLP12 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (273/210);  new base T (2/41);  total (275/252)
Fisher's exact test for biased strand distribution p-value = 3.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 533 total aligned reads displayed.

ATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAA  >  NC_000913/571939‑572236
                                                                                                                                                     |                                                                                                                                                    
aTTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGa                                                                                                                                                      >  2:2023587/1‑150 (MQ=255)
       aGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAgg                                                                                                                                               <  2:1356995/150‑1 (MQ=255)
       aGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAgg                                                                                                                                               <  1:8705374/150‑1 (MQ=255)
         tCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGtt                                                                                                                                             >  2:5205685/1‑150 (MQ=255)
          cGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTa                                                                                                                                            <  1:6122333/150‑1 (MQ=255)
           gttgttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTaa                                                                                                                                           >  2:3624600/1‑150 (MQ=255)
               ttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataa                                                                                                                                       >  2:8935475/1‑150 (MQ=255)
                tataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataat                                                                                                                                      <  1:3194404/150‑1 (MQ=255)
                tataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataat                                                                                                                                      <  2:7198516/150‑1 (MQ=255)
                tataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataat                                                                                                                                      <  2:7018244/150‑1 (MQ=255)
                tataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataat                                                                                                                                      <  2:8482323/150‑1 (MQ=255)
                tataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataat                                                                                                                                      >  2:8636656/1‑150 (MQ=255)
                tataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataat                                                                                                                                      >  2:3451130/1‑150 (MQ=255)
                 ataACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATg                                                                                                                                     <  2:6031291/150‑1 (MQ=255)
                  taACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGc                                                                                                                                    >  1:1293723/1‑150 (MQ=255)
                    aCCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGt                                                                                                                                  <  1:5249194/150‑1 (MQ=255)
                         ctgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCt                                                                                                                             <  1:8219306/150‑1 (MQ=255)
                          tgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTg                                                                                                                            >  1:3936906/1‑150 (MQ=255)
                          tgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTg                                                                                                                            <  2:4267120/150‑1 (MQ=255)
                           gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                                                           <  1:7290098/150‑1 (MQ=255)
                            ctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCt                                                                                                                          >  1:8135358/1‑150 (MQ=255)
                             tATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCta                                                                                                                         <  1:2254700/150‑1 (MQ=255)
                              aTCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCtat                                                                                                                        <  1:2895049/150‑1 (MQ=255)
                               tCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCtata                                                                                                                       >  1:9463325/1‑150 (MQ=255)
                               tCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCtata                                                                                                                       <  1:6455518/150‑1 (MQ=255)
                               tCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCtata                                                                                                                       >  2:9383701/1‑150 (MQ=255)
                               tCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCtata                                                                                                                       <  1:8485740/150‑1 (MQ=255)
                                 cGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACa                                                                                                                     >  2:6319188/1‑150 (MQ=255)
                                 cGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACa                                                                                                                     >  2:8760176/1‑150 (MQ=255)
                                   aaCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACaaa                                                                                                                   <  2:5405244/150‑1 (MQ=255)
                                   aaCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACaaa                                                                                                                   <  1:4117460/150‑1 (MQ=255)
                                        ggCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAacgacg                                                                                                              <  2:2999838/150‑1 (MQ=255)
                                         gCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACgg                                                                                                             >  2:2460638/1‑150 (MQ=255)
                                             aaaaTGGACAGCAAACAAGGTTGCTCGATATCTCAACATTTCTGTTTCTACATTACATCGCCGTCTACCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTAAAAAAATGCGATGTCTGCTATACAAACGACGGTaaa                                                                                                         >  2:4064502/1‑150 (MQ=255)
                                               aaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAAc                                                                                                       <  2:1567121/150‑1 (MQ=255)
                                                    aCAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCtata                                                                                                  >  1:3373191/1‑150 (MQ=255)
                                                      aGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAg                                                                                                >  1:1996411/1‑150 (MQ=255)
                                                        cAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCg                                                                                              >  2:9052859/1‑150 (MQ=255)
                                                         aaaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCga                                                                                             >  2:7311690/1‑150 (MQ=255)
                                                         aaaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCga                                                                                             <  2:4907960/150‑1 (MQ=255)
                                                          aaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCgag                                                                                            >  1:168866/1‑150 (MQ=255)
                                                            cAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGAt                                                                                          >  1:4608906/1‑150 (MQ=255)
                                                             atgtttGCTCTATATATCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGTGGGGGTATGTTTCCAAAGTATACTGGACGATGTTAGTTTAAATAATGCGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGAtt                                                                                         <  2:9277928/146‑1 (MQ=255)
                                                             aaGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGAtt                                                                                         >  2:6439977/1‑150 (MQ=255)
                                                               ggTTGCTATATATCTCTAATTTTCTTTTTCTAAATTACATCGCCTTCTAGCACGGGTGGGGGTATGTTTCAAAAGTTTACTGGACGTTGTGTGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGc                                                                                       <  1:5604331/150‑1 (MQ=255)
                                                                 ttGCTCTATTTCTCTACTTTTCTGTTTATACATTACAACGCCGTCTAGCATGGGTGGGGGTATGTTCCAAAAGTATACTGGACGTTGTGTGTTTAAATATTGCGTTGTCTGCTATAAAAACGACGGTAAAACCTATAAGCGAGATTGCCa                                                                                     <  1:4706175/150‑1 (MQ=255)
                                                                   gCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCaga                                                                                   >  2:1009808/1‑150 (MQ=255)
                                                                       gATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGaaa                                                                               >  2:2808019/1‑150 (MQ=255)
                                                                        atatCTCTACATTTCTTTTACTAAATTACATCGCCGTCTATCAAGGGTGGGGGTATGTTTCAAAATTATCCTGGAGGTTTTGTGTTTAAATAATGGGTTTTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGaaaa                                                                              <  1:3570831/150‑1 (MQ=255)
                                                                          atCTCTACTTTTCTGTTTCTACTTTACATCGCCGTCTAGCAAGCGTGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGTGGTTAAATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATg                                                                            <  2:2547931/150‑1 (MQ=255)
                                                                          atCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATCCAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATg                                                                            <  2:9547790/150‑1 (MQ=255)
                                                                            ctctACATTTCTGTTTCTACATTACACCGCCGTCTAGCAAGCGTGGGGGTATTTTTCCAACGTATACTGGACGATGTGAGTTTAAATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGt                                                                          <  1:1462598/150‑1 (MQ=255)
                                                                             tctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGTTAAGTTTCAAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGtt                                                                         <  2:5872422/150‑1 (MQ=255)
                                                                               tACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTtat                                                                       >  1:9028887/1‑150 (MQ=255)
                                                                                   tttCTGTTTCTAAATTACATCGCCGTCTAGCAAGCGAGGGGGTAATTTTCCAAAGTATACTGGACGATGTGAGTTTAAATAATGCGTTGTATGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAgt                                                                   <  2:8388803/150‑1 (MQ=255)
                                                                                     tCTGTTTCTACATTACATCGCCGTCTAGCAATCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATTCGTTGTCTGCTATACAAACGACGGTAAAACCGATAAGAGAGATTGCCAGAGAAAATGGTTATAAgtgt                                                                 >  2:4774599/1‑150 (MQ=255)
                                                                                          ttCTACATTACATCGCCTTCTAGCAAGCGTGGGGTTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTtc                                                            <  1:2599141/150‑1 (MQ=255)
                                                                                           tCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTtct                                                           >  2:7648653/1‑150 (MQ=255)
                                                                                            cTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTtctc                                                          >  2:6133994/1‑150 (MQ=255)
                                                                                                attacatCGCCGTCTTGCACGCGTGGGGGTAAGTTTCAAATGTTTACTGGACTTTGTGTGTTTAAATAATGCTTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                                                      <  1:8955313/150‑1 (MQ=255)
                                                                                                attacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                                                      >  1:2674599/1‑150 (MQ=255)
                                                                                                 ttacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGtttt                                                     >  2:3305021/1‑150 (MQ=255)
                                                                                                 ttacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGtttt                                                     >  1:1166008/1‑150 (MQ=255)
                                                                                                 ttacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGtttt                                                     >  2:2943870/1‑150 (MQ=255)
                                                                                                  tacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTa                                                    >  2:4572517/1‑150 (MQ=255)
                                                                                                    catCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATATTGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACt                                                  <  2:270598/150‑1 (MQ=255)
                                                                                                        gCCGTCTCGCATGCGTGGGGGTACGTTTCCAACGTATCGGGGAGGTGGTGCGGTTAAATATTGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGaaa                                              <  2:6550117/150‑1 (MQ=255)
                                                                                                        gCCGTCTCGCACGCGTGGGGGTATGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAAGGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGaaa                                              <  1:2411865/150‑1 (MQ=255)
                                                                                                            tCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGTTGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGAtt                                          <  2:5675256/150‑1 (MQ=255)
                                                                                                             cTAGCACGCGAGGGGGTAAGTTTCCAACGTATACTGGACGATGTGCGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGAttt                                         <  1:9241009/150‑1 (MQ=255)
                                                                                                             cTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGAttt                                         <  1:4979644/150‑1 (MQ=255)
                                                                                                              tAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTc                                        >  2:7079002/1‑150 (MQ=255)
                                                                                                               aGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCa                                       >  2:1913413/1‑150 (MQ=255)
                                                                                                                  aaGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTTAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATaa                                    >  2:1287884/1‑150 (MQ=255)
                                                                                                                  aaGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGACTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATaa                                    <  1:2481475/150‑1 (MQ=255)
                                                                                                                   aGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAAt                                   >  1:8567970/1‑150 (MQ=255)
                                                                                                                    gcgcGGGGTTAAGTTTCCAAATTATACTGGACGATGTGTGGTTAAATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATc                                  <  1:9291037/150‑1 (MQ=255)
                                                                                                                        gggggTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGttt                              <  2:2204773/150‑1 (MQ=255)
                                                                                                                        gggggTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGttt                              <  1:3065628/150‑1 (MQ=255)
                                                                                                                         ggggTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGtttt                             >  1:4273874/1‑150 (MQ=255)
                                                                                                                           ggTAAGTTACCAACTTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTaa                           <  2:2599791/150‑1 (MQ=255)
                                                                                                                             tAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAata                         >  1:1071095/1‑150 (MQ=255)
                                                                                                                             tAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAata                         <  1:711486/150‑1 (MQ=255)
                                                                                                                              aaGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAatat                        >  1:275231/1‑150 (MQ=255)
                                                                                                                              aaGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAatat                        >  2:8066068/1‑150 (MQ=255)
                                                                                                                                gTTTCCAATGTATAGTGGACGTTGTGAGGTTAAATAATGCGTGGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATaa                      <  1:7469372/150‑1 (MQ=255)
                                                                                                                                gTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATaa                      <  1:6821126/150‑1 (MQ=255)
                                                                                                                                gTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATAAAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATaa                      <  2:6762470/150‑1 (MQ=255)
                                                                                                                                   tCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAacac                   <  2:1247156/150‑1 (MQ=255)
                                                                                                                                    ccAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACAcc                  >  1:7713759/1‑150 (MQ=255)
                                                                                                                                     cAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCa                 >  1:5486908/1‑150 (MQ=255)
                                                                                                                                     cAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCa                 >  1:1784642/1‑150 (MQ=255)
                                                                                                                                     cAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCa                 >  2:3845072/1‑150 (MQ=255)
                                                                                                                                          tataCTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAagaga            <  1:5741364/150‑1 (MQ=255)
                                                                                                                                           ataCTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAAAACCAagagag           >  1:8027479/1‑150 (MQ=255)
                                                                                                                                            taCTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATCTAACACCAAgagaga          >  1:8369300/1‑150 (MQ=255)
                                                                                                                                            taCTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAgagaga          <  1:1905406/150‑1 (MQ=255)
                                                                                                                                              cTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATa        >  2:1793149/1‑150 (MQ=255)
                                                                                                                                               tGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATaa       <  2:748766/150‑1 (MQ=255)
                                                                                                                                                   cGATGTTATTTTAAATAATGCGTTGTCTGCTATAAAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGaaa   <  2:8207800/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGTTTTAATATAACACCAAGAGAGATAAGAAAA  >  NC_000913/571939‑572236

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: