Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,969 C→A 5.6% M5I (ATG→ATT yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,9690CA 5.6% 2591.4 / 12.4 802M5I (ATG→ATTyagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (452/273);  new base A (5/72);  total (457/345)
Fisher's exact test for biased strand distribution p-value = 7.53e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 814 total aligned reads displayed.

CCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGA  >  NC_000913/281820‑282118
                                                                                                                                                     |                                                                                                                                                      
ccTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGc                                                                                                                                                        >  2:4219567/1‑150 (MQ=255)
 cTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCa                                                                                                                                                       >  1:8674456/1‑150 (MQ=255)
  ttCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCAt                                                                                                                                                      >  2:12940401/1‑150 (MQ=255)
    ccTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATaa                                                                                                                                                    <  1:13787516/150‑1 (MQ=255)
     cTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAg                                                                                                                                                   >  2:11194863/1‑150 (MQ=255)
          cccAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACt                                                                                                                                              >  1:11064169/1‑150 (MQ=255)
          cccAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACt                                                                                                                                              <  2:9263834/150‑1 (MQ=255)
               cgcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATa                                                                                                                                         >  1:2875734/1‑150 (MQ=255)
               cgcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATa                                                                                                                                         >  2:5747522/1‑150 (MQ=255)
                gcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATaa                                                                                                                                        >  1:10142508/1‑150 (MQ=255)
                    gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGTATTTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACTTGTTATCTCTCGAATCCCTGGGCATAAGCGACTCCATAAACgg                                                                                                                                    <  2:9754910/150‑1 (MQ=255)
                    gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACgg                                                                                                                                    >  2:3719573/1‑150 (MQ=255)
                     aGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACggg                                                                                                                                   >  1:7245212/1‑150 (MQ=255)
                      gCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGt                                                                                                                                  >  2:2590276/1‑150 (MQ=255)
                       ccACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGtt                                                                                                                                 >  2:2758102/1‑150 (MQ=255)
                       cacgttGTAGCGGGTGGCTGGTGAACTGCCGACGGTACGGCAGAGGGACCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAAATAAGTCAGTAGTTATTTACTTGGTTTGTCTCGCCTCCCAGGGCATAAGCGACTCCATAAACGGGtt                                                                                                                                 <  2:14950576/146‑1 (MQ=255)
                         aCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCAATAAACGGGTTCt                                                                                                                               >  1:4514656/1‑150 (MQ=255)
                           ttGTAGCGGGTGGCAGGTTAATTGCCGACGCGAGGGCAGTGGGAACGGATGTTCGCCCCGTCCTGCGTGGCGAACAAAAAAAACTCCGTACGTAATGACCTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTa                                                                                                                             <  2:6974381/150‑1 (MQ=255)
                             gTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATg                                                                                                                           >  1:14135392/1‑150 (MQ=255)
                               aGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGcc                                                                                                                         >  2:7370924/1‑150 (MQ=255)
                                    gTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGt                                                                                                                    >  1:9217735/1‑150 (MQ=255)
                                       gctggTTAAATGCCGACGCGACGGCAGTGGGAAGGGATGTTCGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGt                                                                                                                 <  2:4840504/147‑1 (MQ=255)
                                         aGGTGAAATGCCTAAGCGAGGGCAGAGGGAAGGTTTGTTCGCCCCGTCCTGCTAGGCGAAAAAAACAAACTCTGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTaa                                                                                                               <  2:4467219/150‑1 (MQ=255)
                                          ggTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAg                                                                                                              >  1:1190660/1‑150 (MQ=255)
                                          ggTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAg                                                                                                              >  1:5100737/1‑150 (MQ=255)
                                          ggTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAg                                                                                                              >  1:8198586/1‑150 (MQ=255)
                                           gTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                             >  2:8865782/1‑150 (MQ=255)
                                           gTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                             >  1:10283409/1‑150 (MQ=255)
                                             gAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgtgt                                                                                                           >  1:1618176/1‑150 (MQ=255)
                                                 tGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAc                                                                                                       >  2:5555878/1‑150 (MQ=255)
                                                 tGCCGAAGCGACGGATGTGGGACCGGTTGTTCGCCCCGTCCTGCGTGGCTAAAAAAACAACATCAGTCCTTAATGACTTGTTATCCCTCGCACCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAc                                                                                                       <  1:3206979/150‑1 (MQ=255)
                                                  gCCGAAGCGACGGCAGAGGGAAGGGATGTTCGCCCCGTCCTGTGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAcc                                                                                                      <  2:1199453/150‑1 (MQ=255)
                                                  gCCGAAGCGACGGCAGAGGGAACGGTTGTGCGCCCCGTCCTGCGAGGCTAACAAAACAAACGCTGTACTTATTGACTTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAcc                                                                                                      <  2:13433959/150‑1 (MQ=255)
                                                   ccGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCa                                                                                                     >  2:9856336/1‑150 (MQ=255)
                                                   ccGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCa                                                                                                     >  1:2241501/1‑150 (MQ=255)
                                                        gCGACGGCAGCGGGAACGGATGTTCGCCCCGGCCTGCGAGGCGAACAAAACAAACTCAGTACTTAATGACATGTTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTATGTGTCTACCATGTcc                                                                                                <  2:15202374/150‑1 (MQ=255)
                                                            cGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGa                                                                                            >  2:6381121/1‑150 (MQ=255)
                                                             ggCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGaa                                                                                           >  1:5913505/1‑150 (MQ=255)
                                                               cAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACa                                                                                         >  1:3948001/1‑150 (MQ=255)
                                                                agagGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACaa                                                                                        >  2:6018649/1‑150 (MQ=255)
                                                                agagGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACaa                                                                                        >  1:14958870/1‑150 (MQ=255)
                                                                  tgggAAGGGTTGTTCGCCCCGTCCTGCGTGGCGAACAAAAAAAACCCCGTACTTAATGACATGGTATCTCTCGCATCCCTGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgt                                                                                      <  1:4881099/149‑1 (MQ=255)
                                                                    ggAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgtgt                                                                                    >  2:3116682/1‑150 (MQ=255)
                                                                      accGGTTGTGCGCCCCGTCCGGCGGGGCTAACAAAAAAAACTCTGTAGGTATTGACTGGGTATCTCGCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTc                                                                                  <  1:6414540/148‑1 (MQ=255)
                                                                      aaCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTc                                                                                  <  2:9111713/150‑1 (MQ=255)
                                                                      aaCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCGTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTc                                                                                  >  1:6847559/1‑150 (MQ=255)
                                                                       aCGTATGTTCGCCCCGACCTGCGAGGCTAAAAAAACAAACTCAGTACGTAATGACATGTTATCTCTCGCACCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCa                                                                                 <  1:13000871/150‑1 (MQ=255)
                                                                        cGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCAc                                                                                >  2:8008323/1‑150 (MQ=255)
                                                                         ggATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAAGGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCGTGTCCCCGAACTAGTGTTCACt                                                                               >  2:11875086/1‑150 (MQ=255)
                                                                               tCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTcc                                                                         >  2:671401/1‑150 (MQ=255)
                                                                               tCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTcc                                                                         >  1:10554033/1‑150 (MQ=255)
                                                                               tCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTcc                                                                         >  2:12491226/1‑150 (MQ=255)
                                                                                cGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTccc                                                                        <  2:5054685/150‑1 (MQ=255)
                                                                                 gCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTcccc                                                                       <  2:11520869/150‑1 (MQ=255)
                                                                                   cccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCgg                                                                     >  1:11469352/1‑150 (MQ=255)
                                                                                     cGTCCTGCGAGGCGAACAAAAAAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGAc                                                                   <  2:1944581/150‑1 (MQ=255)
                                                                                      gTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGAcc                                                                  >  1:1910660/1‑150 (MQ=255)
                                                                                       tCCGGCGAGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCACCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCg                                                                 <  2:9442172/150‑1 (MQ=255)
                                                                                          tgctatccGAAAAAAAAAAACTCAGTACGTAATGAAATGGAATATCTCGAAACCAAGGCAAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTac                                                              <  1:15167860/143‑1 (MQ=255)
                                                                                              aGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCC                                                          <  2:14301030‑M1/150‑1 (MQ=255)
                                                                                                   aacaaaacaaACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAG                                                    >  2:774260‑M1/1‑150 (MQ=255)
                                                                                                      aaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGG                                                 >  2:2921601‑M1/1‑150 (MQ=255)
                                                                                                       aaacaaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGga                                                >  2:10534053/1‑150 (MQ=255)
                                                                                                        aacaaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGgag                                               >  2:3937734/1‑150 (MQ=255)
                                                                                                          caaacCCAGTACGTAATGACATGGTATCACTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGagag                                             <  1:11680327/150‑1 (MQ=255)
                                                                                                           aaacTCAGTACGTATTGACATGGTATCTCTCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAgg                                            <  1:9791661/150‑1 (MQ=255)
                                                                                                               tCAGTACGTAATGACGTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGAc                                        >  1:12527384/1‑150 (MQ=255)
                                                                                                               tCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGGAGAGGGGa                                         >  1:11292697/1‑150 (MQ=255)
                                                                                                               tCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGAc                                        >  1:8981449/1‑150 (MQ=255)
                                                                                                               cccgTACGTAATGACCTGTTATCTCTCGCATCCCCGGGCAAACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGAc                                        <  2:1454767/147‑1 (MQ=255)
                                                                                                                 aGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTg                                      >  1:751890/1‑150 (MQ=255)
                                                                                                                       tAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAggg                                           >  1:362704/1‑139 (MQ=255)
                                                                                                                          tGACAGGGTATCTCCCGCAGCCCAGGGCTTAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGcc                             <  2:4143886/150‑1 (MQ=255)
                                                                                                                            aCATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCAt                           <  2:2326189/150‑1 (MQ=255)
                                                                                                                            aCATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCAt                           >  2:14003385/1‑150 (MQ=255)
                                                                                                                              aTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCt                         <  1:3719573/150‑1 (MQ=255)
                                                                                                                                ggTATCTCTCGCATCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCttt                       <  1:15291460/150‑1 (MQ=255)
                                                                                                                                ggTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAACCGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCttt                       <  2:15436005/150‑1 (MQ=255)
                                                                                                                                ggTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGtt                                                                                   >  1:7909374/1‑91 (MQ=255)
                                                                                                                                 gTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCtttt                      >  2:10790449/1‑150 (MQ=255)
                                                                                                                                 gTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCtttt                      >  1:13506992/1‑150 (MQ=255)
                                                                                                                                  tATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTc                     >  2:13646448/1‑150 (MQ=255)
                                                                                                                                      tctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccccc                 >  1:10722786/1‑150 (MQ=255)
                                                                                                                                      tctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccccc                 >  2:528061/1‑150 (MQ=255)
                                                                                                                                       ctcGCCTCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCt                <  1:5763631/150‑1 (MQ=255)
                                                                                                                                       ctcGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCt                <  2:11064169/150‑1 (MQ=255)
                                                                                                                                        tcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc               >  1:4652815/1‑150 (MQ=255)
                                                                                                                                        tcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc               >  1:2538558/1‑150 (MQ=255)
                                                                                                                                        tcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTAGGCCGTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc               >  1:9482784/1‑150 (MQ=255)
                                                                                                                                            aTCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccc           >  2:2073322‑M2/1‑147 (MQ=255)
                                                                                                                                               ccAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctt        <  1:12174147‑M2/150‑7 (MQ=255)
                                                                                                                                                cAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt       <  1:7018788‑M2/150‑8 (MQ=255)
                                                                                                                                                cAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt       >  2:859866‑M2/1‑143 (MQ=255)
                                                                                                                                                 aGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttg      >  1:12699790‑M2/1‑142 (MQ=255)
                                                                                                                                                 aGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttg      <  1:14588684‑M2/150‑9 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCCAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  2:8185438‑M2/1‑141 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  2:1623305‑M2/1‑141 (MQ=255)
                                                                                                                                                   ggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggg    >  2:6532316‑M2/1‑140 (MQ=255)
                                                                                                                                                    gCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgggg   <  1:547484‑M2/150‑12 (MQ=255)
                                                                                                                                                     cATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGga  <  1:9050545/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGA  >  NC_000913/281820‑282118

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: