Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 287,437 T→A 5.2% S217T (TCC→ACC)  yagH → CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,4370TA 5.2% 249.3 / 3.4 77S217T (TCC→ACC) yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base A (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

CACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATGGCC  >  NC_000913/287289‑287574
                                                                                                                                                    |                                                                                                                                         
cACGCGCCCGCCCGACCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACCCCGAACGCGCGGATACGGATCGCCACGCGGGTTGGTACAACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGtt                                                                                                                                           <  5:280491/149‑1 (MQ=255)
     tctcGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCaa                                                                                                                                       >  3:75025/1‑148 (MQ=255)
     tctcGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCaa                                                                                                                                       >  2:81101/1‑148 (MQ=255)
             gAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCTCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGa                                                                                                                              <  8:112595/149‑1 (MQ=255)
               gcgcAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACg                                                                                                                            <  4:325357/149‑1 (MQ=255)
                   aaaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGcc                                                                                                                        >  8:371744/1‑149 (MQ=255)
                   aaaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTACCAAAAATATCGACGGGcc                                                                                                                        >  4:82855/1‑149 (MQ=255)
                   aaaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTACCAAAAATATCGACGGGcc                                                                                                                        >  4:78317/1‑149 (MQ=255)
                    aaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCg                                                                                                                       >  4:337356/1‑149 (MQ=255)
                    aaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTCCGTTCCAAAAATATCGACGGGCCg                                                                                                                       >  4:342699/1‑149 (MQ=255)
                     aaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGt                                                                                                                      <  7:263921/149‑1 (MQ=255)
                            tttACTGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGc                                                                                                                <  3:70688/148‑1 (MQ=255)
                            tttACCGGCACGCCGCTCTGCTACACCGAAGGTGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGc                                                                                                                <  2:31294/148‑1 (MQ=255)
                                 cGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCAc                                                                                                           <  8:314496/148‑1 (MQ=255)
                                 cGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCAc                                                                                                           >  7:314496/1‑148 (MQ=255)
                                  ggCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCAccc                                                                                                         <  2:227661/149‑1 (MQ=255)
                                      cgccgcTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGAc                                                                                                     <  6:151339/149‑1 (MQ=255)
                                        ccgcTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCAGCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGt                                                                                                   <  2:104410/149‑1 (MQ=255)
                                          gcTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTaa                                                                                                 <  2:364260/149‑1 (MQ=255)
                                          gcTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTaa                                                                                                 >  6:217353/1‑149 (MQ=255)
                                            tctGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAAc                                                                                                >  2:35317/1‑148 (MQ=255)
                                            tctGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAAc                                                                                                >  3:348930/1‑148 (MQ=255)
                                              tGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAg                                                                                                                 >  8:201099/1‑129 (MQ=255)
                                              tGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAg                                                                                                                 <  7:201099/129‑1 (MQ=255)
                                              tGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgat                                                                                             >  6:31115/1‑149 (MQ=255)
                                               gCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGGTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatg                                                                                            <  6:141509/149‑1 (MQ=255)
                                                   cacCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatgat                                                                                          <  4:217933/147‑1 (MQ=255)
                                                   cacCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGAtgatga                                                                                        >  2:175646/1‑149 (MQ=255)
                                                   cacCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCACCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatgat                                                                                          <  6:243604/147‑1 (MQ=255)
                                                    acCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACga                                                                                              >  6:216950/1‑142 (MQ=255)
                                                    acCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACga                                                                                              <  5:216950/142‑1 (MQ=255)
                                                    acCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGAc                                                                                       <  2:246845/149‑1 (MQ=255)
                                                     ccGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAACTATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGAcc                                                                                      <  6:9446/149‑1 (MQ=255)
                                                     ccGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGt                                                                                                   <  1:171759/136‑1 (MQ=255)
                                                     ccGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGt                                                                                                   >  2:171759/1‑136 (MQ=255)
                                                      cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCCGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGCCGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGAc                                                                                       >  4:287718/1‑147 (MQ=255)
                                                         aGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCACCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCt                                                                                  <  8:257245/149‑1 (MQ=255)
                                                          ggCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTg                                                                                 <  3:270347/149‑1 (MQ=255)
                                                                  ccTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGc                                                                         <  5:299447/149‑1 (MQ=255)
                                                                    tGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCTCGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatg                                                                                            <  5:110048/128‑1 (MQ=255)
                                                                    tGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatg                                                                                            >  6:110048/1‑128 (MQ=255)
                                                                      tATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTc                                                                                                                                          <  7:306568/80‑1 (MQ=255)
                                                                      tATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTc                                                                                                                                          >  8:306568/1‑80 (MQ=255)
                                                                      tATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTACCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGga                                                                     <  3:82855/149‑1 (MQ=255)
                                                                      tATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTACCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGga                                                                     <  3:78317/149‑1 (MQ=255)
                                                                        tCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGgag                                                                    <  1:89189/148‑1 (MQ=255)
                                                                          gCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGaa                                                                  <  1:249588/148‑1 (MQ=255)
                                                                               gcgGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCt                                                            >  1:63097/1‑149 (MQ=255)
                                                                               gcgGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCt                                                            <  5:31115/149‑1 (MQ=255)
                                                                                     tGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGACGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCag                                                        >  6:204010/1‑147 (MQ=255)
                                                                                      ggTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaag                                                     >  4:36634/1‑149 (MQ=255)
                                                                                      ggTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACAAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaa                                                      <  5:204185/148‑1 (MQ=255)
                                                                                        tactacCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAagag                                                   >  3:286007/1‑149 (MQ=255)
                                                                                           tacCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCgg                                                >  8:143997/1‑149 (MQ=255)
                                                                                            acCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGc                                               >  6:22732/1‑149 (MQ=255)
                                                                                               tGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCAc                                            >  6:27699/1‑149 (MQ=255)
                                                                                               tGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCAc                                            >  2:28467/1‑149 (MQ=255)
                                                                                                      cgccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCCCCCGGACGGAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgc                                     >  2:92709/1‑149 (MQ=255)
                                                                                                        ccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgctg                                   <  5:27699/149‑1 (MQ=255)
                                                                                                               cggcACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGa                            >  6:276869/1‑149 (MQ=255)
                                                                                                                ggcACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGAc                           <  1:74264/149‑1 (MQ=255)
                                                                                                                ggcACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCCCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGAc                           >  1:10240/1‑149 (MQ=255)
                                                                                                                     cAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATa                      >  1:90818/1‑149 (MQ=255)
                                                                                                                     cAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATa                      >  1:341558/1‑149 (MQ=255)
                                                                                                                     cAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATa                      <  7:143997/149‑1 (MQ=255)
                                                                                                                       gCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACg                    >  2:166945/1‑149 (MQ=255)
                                                                                                                       gCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTC‑‑‑GCTGCAGACGCATACGGGt                 <  1:28735/149‑1 (MQ=255)
                                                                                                                           cGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACggg                  <  4:75025/147‑1 (MQ=255)
                                                                                                                           cGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGt                 <  1:166945/148‑1 (MQ=255)
                                                                                                                             aGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGCCGGGCCGTACGAGCGGCACCCGGACCTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTAGCAGCTGCAGACGCATACGGGTGaa              >  1:88420/1‑149 (MQ=255)
                                                                                                                               cACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATg            <  7:321457/149‑1 (MQ=255)
                                                                                                                                  gCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTa         <  3:287718/149‑1 (MQ=255)
                                                                                                                                    cgtcgtGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACa       <  5:204010/149‑1 (MQ=255)
                                                                                                                                      tcgtGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATg     <  5:264703/149‑1 (MQ=255)
                                                                                                                                      tcgtGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGAAGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATg     <  5:276869/149‑1 (MQ=255)
                                                                                                                                      gcgtGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATg     <  6:55105/148‑1 (MQ=255)
                                                                                                                                        gtggtgCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATgg    <  4:348930/148‑1 (MQ=255)
                                                                                                                                         tggtgCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATGGcc  <  2:90818/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
CACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGCTGCAGACGCATACGGGTGAATGGTACATGGCC  >  NC_000913/287289‑287574

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