Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,965 A→T 7.4% W7R (TGG→AGG)  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,9650AT 7.4% 783.0 / 4.8 405W7R (TGG→AGG) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (161/178);  new base T (2/44);  total (164/241)
Fisher's exact test for biased strand distribution p-value = 1.36e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 423 total aligned reads displayed.

TGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCT  >  NC_000913/281826‑282104
                                                                                                                                           |                                                                                                                                           
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGTACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  2:1580249/140‑1 (MQ=255)
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  2:1320667/140‑1 (MQ=255)
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  2:1537976/140‑1 (MQ=255)
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  2:905898/140‑1 (MQ=255)
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  1:512633/140‑1 (MQ=255)
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  2:359407/140‑1 (MQ=255)
tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  1:777328/140‑1 (MQ=255)
agAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                             <  1:279789/139‑1 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  1:714745/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  2:629683/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  1:599399/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  2:520995/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  2:1168966/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  1:1611575/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  1:377629/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATa                                                                                                                                      >  1:1546634/1‑140 (MQ=255)
       aGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCCCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGGAAAGACATGGGATCGCTCGCTTCCCAGGGCATa                                                                                                                                      >  2:218565/1‑140 (MQ=255)
        gcgcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATaa                                                                                                                                     >  1:253474/1‑140 (MQ=255)
         cgcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGACCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCCTAAg                                                                                                                                    >  1:694617/1‑140 (MQ=255)
            tGGTGCCACTTGTAGCCGGTGGCTGGTGAAAGGCCGAAGCGACGGCAGAGGGAACGGCTGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACTTGGTTTCTCTCGCATCCCTGGGCTTAAGCGa                                                                                                                                 <  2:263192/140‑1 (MQ=255)
            tGGAGCCAATTGTAGCCGGTGGCTGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTACGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGa                                                                                                                                 <  2:728544/140‑1 (MQ=255)
              gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACt                                                                                                                               >  1:1318237/1‑140 (MQ=255)
              gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACt                                                                                                                               >  1:622838/1‑140 (MQ=255)
                   aCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTCTCTCTCGCATCCCAGGGCATAAGCGACTCCATa                                                                                                                          >  2:697944/1‑140 (MQ=255)
                         tgCCGGTGGCAGGTGAATTGCCGAAGCGACGTCAGAGGTAACTGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCATTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACggg                                                                                                                    <  2:799991/139‑1 (MQ=255)
                         cgCCGGTGGCTGGTTAAATGCCGAAGGGACGGCAGTGGGAAGGTATGTTCGCCCCGTCCTGCGAGGCGAACAAAAAAACCACAGTACGTAATGAATTGGTATCTCTCGAATCCCAGGGATTAAGCGACTCCATAAACggg                                                                                                                    <  1:1165971/139‑1 (MQ=255)
                         aGCCGGTGGCAGGTCAAAAGCCGAAGCGAGCGCAGAGGCAACGTATGTTCGCCCCGACCTGCCAGGCCAACAAAACAAACTAAGAACGAAAAGACATGGTATCACACCCAACCCAGGCCATAAGCGACTCCATAAACggg                                                                                                                    <  2:881689/140‑1 (MQ=255)
                           ccGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACggggt                                                                                                                  >  2:1396094/1‑138 (MQ=255)
                               tGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAAGGGTTGTTCGCCCCGTCCTGCGTGGCGAACAAAACAAACTCAGTACTAAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTa                                                                                                              <  2:412781/140‑1 (MQ=255)
                               tGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTa                                                                                                              >  1:740569/1‑140 (MQ=255)
                               tGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGGATCGCGCGCATCCCAGGGCCTAAGCGACGCCCTCAACgggggcgga                                                                                                              >  2:1078571/1‑134 (MQ=255)
                                   tggTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTACTTAATGACAGGTTATCTCTCGTATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGcc                                                                                                          <  1:1402704/139‑1 (MQ=255)
                                   cggTGAAATGCCGACGCGCCGGCAGCGGGAACGGCTGTTCGCCCCGCCCTGCGAGGCTAACAAACCAAACTCCGTACGTAAGGACTTGGTATCTCTCGCATCCCCGGGCAAAAGCGACTCCATAAACGGGTTCTTATGcc                                                                                                          <  1:919997/139‑1 (MQ=255)
                                   aGGTTAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGcc                                                                                                          <  2:1041633/140‑1 (MQ=255)
                                   aGGTGAAATGCCGACGCGCCGGCTGAGGGAACGGATGTACGCCCCGTCCTGCGAGGCTAAAAAAACAAACTCAGTACTTAATGACATGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGcc                                                                                                          <  1:239088/140‑1 (MQ=255)
                                     gttaattGCCGAAGCGACGGCAGTGGTAACGGATGTTCGCCCCGTCCTGCGTGGCTAAAAAAACAAACACAGTACTTAATGAAATGGTATCTCTCGCATCCCAGGTAATACGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                        <  1:303307/134‑1 (MQ=255)
                                     gttaattGCCGAAGCGACGGATGTGGGAACGGATTTTCGCCCCTTCCTGCGCGGCCAACAAACCAAACTCGTTACTTAATGACATGTTATCTCTCGCAACCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                        <  1:65470/134‑1 (MQ=255)
                                     gTGAAATGCCGAAGCGACGGCAGAGGGACCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCAGAAGCGACTCCAGAAACGGGGTCTTATGccgt                                                                                                        >  2:652761/1‑138 (MQ=255)
                                     gTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAAAAAAACAAACTCATTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                        <  2:508988/140‑1 (MQ=255)
                                     gTGAAAAGCACAAGCGACGGCAGAGGGAACGGATGATCGCCCCGACCTGCGAGGCGAACAAATCAAACTCTCTCCTTAATTACATGGTATCTCTCCCATCCCAGGGCAAAAGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                        <  2:1691322/140‑1 (MQ=255)
                                           ttccGAAGCGAGGGCTGAGGGAACGTATGTTCGCCCCGACCTGCGAGGCTAACAAAACAACATCAGTACGTAATGACATGTTAACACACGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTATTATGACTTAGTTGt                                                                                                  <  2:249986/138‑1 (MQ=255)
                                                aaGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgtg                                                                                             >  2:314048/1‑140 (MQ=255)
                                                  gCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGACCGTAATGACATGTTATCTCGCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTc                                                                                           <  2:1628876/140‑1 (MQ=255)
                                                      cGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAcc                                                                                       <  2:1476133/140‑1 (MQ=255)
                                                       ggCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCTAAAAAAACAAACACAGTACGTAATGACATGGTATCTCTCGCAGCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCa                                                                                      <  1:681064/140‑1 (MQ=255)
                                                          agagGTAACGGATGTACGCCCCGTCCTGCGAGGCGAACAAACAAAACTCAGTACGTAATGAAATGGTATCACTCGCATCCCAGGGCATAAGCGACTCCTTACCCGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGt                                                                                   <  2:1245625/140‑1 (MQ=255)
                                                            agGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATCAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTTCCATGTcc                                                                                 >  2:775721/1‑140 (MQ=255)
                                                                aaCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGa                                                                             >  1:1120115/1‑140 (MQ=255)
                                                                 aCGGATTTTCGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTAGGTACTGAAATGGTTACTCTCGCATCCCTGGGCAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGaa                                                                            <  1:1225628/140‑1 (MQ=255)
                                                                 aCGGATGTTCGCCCGGTCCTGCGAGGCTAACAAAAAAAACACTGTACGTATTGACATGGTATCTCTCGCAACCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGaa                                                                            <  1:450831/140‑1 (MQ=255)
                                                                 aCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATTACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGaa                                                                            <  2:1495654/140‑1 (MQ=255)
                                                                  cGGTTGTGCGCCCCGTCCTGCGTGGCTAAAAAAACAAACACATTACTTAATCACATGTTATCTCACCCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAAc                                                                           <  1:989788/140‑1 (MQ=255)
                                                                  cGGATGTTCGCCCCGACCTGCGAGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAAc                                                                           <  1:576286/140‑1 (MQ=255)
                                                                    gATTTTCGCCCCGTCCTGCGTGGCTAACAAAACAAACTCAGTACTTAATGAAATGGTATCTCTCGCATCCCAGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACaa                                                                         <  1:151642/140‑1 (MQ=255)
                                                                    gATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACaa                                                                         <  2:1546634/140‑1 (MQ=255)
                                                                         tCGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGtt                                                                    <  2:101996/140‑1 (MQ=255)
                                                                         tCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGtt                                                                    <  1:324840/140‑1 (MQ=255)
                                                                              ccttCCTGCGAGGCTAAAAAAACAAACTCATTACTTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTa                                                               <  2:403928/137‑1 (MQ=255)
                                                                              ccGTCCTGCGAGGCGAACACACCAAACTCAGAACTTAATGACATGTTATCTCTCGCACCCCAGGGCAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACtc                                                               <  1:955215/140‑2 (MQ=255)
                                                                              ccGTCCTGCGAGGCGAACAACCCCAACTCAGTACGTAATTACTTGGTTTCACTCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTa                                                               <  2:350786/140‑1 (MQ=255)
                                                                              ccGTCCTGCGAGGCGAACAAAACAACCTCAGTACGTAATGACATGGTATCTCTCGCAACCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTa                                                               <  1:1055385/140‑1 (MQ=255)
                                                                              ccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGTTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTa                                                               <  1:1454728/140‑1 (MQ=255)
                                                                                gTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATg                                                             >  2:1294556/1‑140 (MQ=255)
                                                                                 tCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGt                                                            >  2:827302/1‑140 (MQ=255)
                                                                                 tCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGt                                                            >  1:1688029/1‑140 (MQ=255)
                                                                                      cGAGGCTAACAAACCAAACACTTTACTTAATTACAGGTTATCTCTCGAATCCCCGGGCTTAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCg                                                       <  2:191164/140‑1 (MQ=255)
                                                                                        tggCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGa                                                     <  2:89108/139‑1 (MQ=255)
                                                                                        tggCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGa                                                     <  2:594267/139‑1 (MQ=255)
                                                                                        cggCGCACAAAACCACCTCAGTCCGCACTGCCCTGGTCTCTCTCGCAGCCCTGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGa                                                     <  1:1078571/139‑1 (MQ=255)
                                                                                           cGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCg                                                  >  2:754020/1‑140 (MQ=255)
                                                                                           cGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCg                                                  >  2:549768/1‑140 (MQ=255)
                                                                                           cGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCGTAGTTGTCAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCg                                                  >  2:1062390/1‑140 (MQ=255)
                                                                                             aacaaaccaaCCTCCGTACGTAATGACTGGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTa                                                <  1:533751/140‑1 (MQ=255)
                                                                                             aacaaaacaaCCTCATTACGTACGGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTa                                                <  1:1001681/140‑1 (MQ=255)
                                                                                             aacaaaaaaaACTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTa                                                <  2:225037/140‑1 (MQ=255)
                                                                                               caaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTaca                                              <  2:456155/140‑1 (MQ=255)
                                                                                                 aaacaaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACAcc                                            <  1:841534/140‑1 (MQ=255)
                                                                                                    caaacTCAGTACGTAATGACTTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCa                                         <  1:1585579/140‑1 (MQ=255)
                                                                                                    caaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCa                                         <  1:422172/140‑1 (MQ=255)
                                                                                                    caaacCCCGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCa                                         <  1:1434615/140‑1 (MQ=255)
                                                                                                      aacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCaaa                                       <  2:868242/140‑1 (MQ=255)
                                                                                                             tACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGgaga                                <  1:1345866/140‑1 (MQ=255)
                                                                                                             cccGTAATCCCATGGAAGCACCCGCATCCCCGGGCAAACGCGACTCCAAACACGCGATCTAATGCCTTAGTTGAAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGgaga                                <  2:1257860/138‑1 (MQ=255)
                                                                                                                  aaTGACTTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGa                           <  1:163075/140‑1 (MQ=255)
                                                                                                                  aaTGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGa                           >  1:196856/1‑140 (MQ=255)
                                                                                                                  aaTGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGa                           <  2:1134461/140‑1 (MQ=255)
                                                                                                                  aaTGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCGTATGCGTTAGTTGTAAGTGTCTACCATGTCCACGAACAAGTGGTCTCTTTGTCCCCGGACCGGACACCCCAAGGGGGCGTGGGGa                           >  2:1217551/1‑140 (MQ=255)
                                                                                                                       cagggTATCTCTCGAATCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCa                      <  1:1258126/137‑1 (MQ=255)
                                                                                                                       cATGGTCTCTCTGGCATCCCAGGGCATTAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCa                      <  2:302649/140‑1 (MQ=255)
                                                                                                                       cATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTACTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGAGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCa                      >  2:393412/1‑140 (MQ=255)
                                                                                                                         tGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCAcc                    >  1:1386530/1‑140 (MQ=255)
                                                                                                                          ggTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCg                   >  2:1087065/1‑140 (MQ=255)
                                                                                                                          ggTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCg                   >  2:1620299/1‑140 (MQ=255)
                                                                                                                          ggTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCg                   >  2:1366806/1‑140 (MQ=255)
                                                                                                                              tctctcGCAGCCCAGGGCATAAGCGACTCCATAACCGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGcc               <  1:1499061/140‑1 (MQ=255)
                                                                                                                                 ctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATc            >  1:1111232/1‑140 (MQ=255)
                                                                                                                                    gCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCttt         >  2:962160/1‑140 (MQ=255)
                                                                                                                                     cATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCtttt        <  2:582663/140‑1 (MQ=255)
                                                                                                                                     cATCCCAGGGCATAAGCGACTCCATAAACGGATTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTa                                                >  1:899649/1‑100 (MQ=255)
                                                                                                                                      aTCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTc       >  2:1402239/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCT  >  NC_000913/281826‑282104

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