Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,363,184 C→T W42* (TGG→TAG)  yjdC ← putative DNA‑binding transcriptional regulator YjdC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,363,1840CT97.3% 207.5 / ‑3.9 73W42* (TGG→TAG) yjdCputative DNA‑binding transcriptional regulator YjdC
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (42/29);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGC  >  NC_000913/4362943‑4363427
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
tttgtGCTGATCGGCCAGGTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCATGCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACAGCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAgcg                                                                                                                                                                                                                                               <  2:543935/245‑1 (MQ=255)
         aTCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAAAGGATTGCCTCTTTGTCTTGCTAGAAGAGGCGTAGCTCGt                                                                                                                                                                                                                                   >  1:340448/1‑251 (MQ=255)
                          aaGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCcaca                                                                                                                                                                                                                  <  1:508475/250‑1 (MQ=255)
                            gggTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACg                                                                                                                                                                                                                >  1:117872/1‑251 (MQ=255)
                            gggTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACg                                                                                                                                                                                                                <  1:637843/251‑1 (MQ=255)
                                    aGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCa                                                                                                                                                                                                         >  1:560073/1‑250 (MQ=255)
                                      gATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCGTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAAcc                                                                                                                                                                                                      <  1:27455/251‑1 (MQ=255)
                                         cGGGATAAAACGTACAGGCAGCGATAAATCGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCATTGGCTAATCCCCACGTTCCGCAAccctc                                                                                                                                                                                                   >  1:220798/1‑248 (MQ=255)
                                            gATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCCTCCCGCAGACGTCGACGACGCCAGCCATCCATCTGGTGGCTAAGATAACGCAGCGCAGTGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGACGGAGACCATCCAGTGGAAAATCCACACGTTCataaacacactac                                                                                                                                                                                                >  2:323881/1‑239 (MQ=255)
                                                      aCAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCCCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGTTAATCCACAAGTTCAGCAACCATCTCCAGCGGGgtgt                                                                                                                                                                                      >  1:151694/1‑251 (MQ=255)
                                                       cAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGAGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCCTCAGCTGACGACGACAGCCCTCTATCCGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTACCTCTTTGTCTGGCTTAAAGCGGCGTAGCTCGTCCAGGGGATAATCCACACGTTCCACAACCAGCTcccgccgggtgtt                                                                                                                                                                                     >  2:110131/1‑240 (MQ=255)
                                                          gCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  <  2:316400/251‑1 (MQ=255)
                                                           cAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  >  1:321650/1‑250 (MQ=255)
                                                                aTAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTTTCTGTTGGCTAAGATAACGCAGCGCCACGTAGAGGATTTCCTCTTTGTTTGGCTAGAAAAGGCGGAGCTCGACCCGGGGGAAATTCACAAGTTTCGCCACCATCTCCAGCggggggtttaatagga                                                                                                                                                                             >  2:576026/1‑235 (MQ=255)
                                                                aTAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGtt                                                                                                                                                                                     <  1:576026/242‑1 (MQ=255)
                                                                          aGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCCGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGACGAGGATTGCCTCTTTGTCTGGCTTGAAACGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCt                                                                                                                                                                  >  1:673764/1‑251 (MQ=255)
                                                                            cccGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTTCGTTTCGTCCAGCAACAGCTGACGACGCCCAACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGAATGCCTCTTTTTCTGGCTAGAAGCGGCGTGAGTCGGCCAGTGGGTAATACCCACGGTTAGCAACCATCTCCAGCGTGGGGTTGGCAATCCCTTGTAATTCtt                                                                                                                                                                 >  2:659017/1‑249 (MQ=255)
                                                                             ccGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTGAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGGTAATCCACACGTTCAGCAACCATCTCCAGCGTGGGGTTGGCAAACCCTTGTAATTCtaat                                                                                                                                                               >  1:501213/1‑251 (MQ=255)
                                                                              cGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAACGGCGTAGCTCGTCCAGTGGATAATCCACACGTTAAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCtaat                                                                                                                                                               >  1:567780/1‑250 (MQ=255)
                                                                                gATAGCGGGTGTTTTTAACGACCTCCGATAACGCCTGATAACGTGCCAGCAGCTTGTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCGGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTATAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTaataat                                                                                                                                                            <  2:451049/251‑1 (MQ=255)
                                                                                     cgggtgttTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAg                                                                                                                                                       <  1:170822/247‑1 (MQ=255)
                                                                                     cGGTTGTTTTGAACGCACTCCTATAACGCCTGATTACGTGCCTGCAGTTTTTGGTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAg                                                                                                                                                       <  2:594919/251‑1 (MQ=255)
                                                                                       gttgttTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  >  1:381776/1‑222 (MQ=255)
                                                                                       gttgttTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  <  2:381776/222‑1 (MQ=255)
                                                                                                aaCGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTcccc                                                                                                                                            >  1:569465/1‑250 (MQ=255)
                                                                                                     tcTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCCGACATCTATCTGTTGGCTAAGATTACGCAACGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAAGGGCGTAGCTCGTCCCGTGGATAATCCACACGTTCAGAAACCATCTCCAGCGTGGGGTTGGCAATCCTTTGTAATTCTAATAATTTCAGGGCTTCTCCCAggac                                                                                                                                        >  2:428717/2‑247 (MQ=255)
                                                                                                          gATAACGCCTGATAACGTGCCATCAGCTTTTGTGCGGCGGTTTGCGTTTCGTCCAGCATCAGCCGACGACGCCAGACAGCTATCTGTTGGCTGAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATCTCAGGGCTTCTCCCAGTACATCtgc                                                                                                                                  <  2:501213/251‑3 (MQ=255)
                                                                                                             aaCGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACg                                                                                                                               >  1:313492/1‑251 (MQ=255)
                                                                                                               cGCCTGATAACGTGCCAGCAGCTTTGTTTCGGCGGTTCGCGTGTCGGCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGTCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTTGCTCGTCCAGTGGATAATCCTCACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGtt                                                                                                                             <  2:55724/251‑1 (MQ=255)
                                                                                                                gCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTg                                                                                                                            >  1:320070/1‑251 (MQ=255)
                                                                                                                    gATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCAAACGTTCAGCCACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACg                                                                                                                        >  1:348782/1‑251 (MQ=255)
                                                                                                                    gATAACGTGCCAGCAGCTTTGGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCAc                                                                                                                         <  2:320070/250‑1 (MQ=255)
                                                                                                                        aCGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGGTTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCCGCAACCATCTCCAGCGTGGGGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTTCACGCTAt                                                                                                                    >  1:272165/1‑251 (MQ=255)
                                                                                                                          gTGCCAGCAGCTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTAtttt                                                                                                                 <  1:45502/251‑1 (MQ=255)
                                                                                                                                    cTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCgtggt                                                                                                                                                                                        >  1:92804/1‑171 (MQ=255)
                                                                                                                                    cTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGGAGCGCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCgtggt                                                                                                                                                                                        <  2:92804/171‑1 (MQ=255)
                                                                                                                                        tGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTc                                                                                                     <  2:313492/250‑1 (MQ=255)
                                                                                                                                        tGTTCGGCGGTTTGCGTTTCGGCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTc                                                                                                     <  2:569465/250‑1 (MQ=255)
                                                                                                                                          gtcggcgggttggtcgtttgccAGCAACTGCTGACGACGCCAGACAGCTATCGTTGGGCTAAGATATCTCATCGCAGCGTAGAGGATTGCTGCTTTGTTTGGCTAGAAGCGGCGTAGCTCGTTCAGTGGAGAAGCCACACGTTCAGCAATCATCTCCAGCGTGGTGTAGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTTCTGAATCACGGGGCACGCTATTTTTCGCCGTCTCTCCCa                                                                                                  <  2:256322/231‑1 (MQ=255)
                                                                                                                                                ggTTTGCGTTTCGTCCAGCATCAGCTGACGACGCAAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACCTCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCa                                                                                             >  2:470756/1‑250 (MQ=255)
                                                                                                                                                  tttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTCGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCAAGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAgt                                                                                          >  2:227306/1‑251 (MQ=255)
                                                                                                                                                   ttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTAGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTTCATCTTCACGTTGCACGCTATTTTCCTCAGTCTTTccccatgcacgtg                                                                                         >  2:285962/1‑241 (MQ=255)
                                                                                                                                                    tGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAgtgt                                                                                        >  1:541319/1‑251 (MQ=255)
                                                                                                                                                     gCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTc                                                                                       >  2:181887/1‑251 (MQ=255)
                                                                                                                                                      cGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCAACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCCCTGCAAGTGTCg                                                                                      >  2:622303/1‑251 (MQ=255)
                                                                                                                                                       gTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGAAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGAAGCTCGTCCAGTGGATAATCCACACGTTCCGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTCCCAGTGTCGt                                                                                     >  2:196909/1‑251 (MQ=255)
                                                                                                                                                            gTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACg                                                                                >  2:308215/1‑251 (MQ=255)
                                                                                                                                                              ccAGCATCCGCTGACGACGCCAGACATCTATCTGTTGGCTAACATACCGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAACGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATACCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCACACTGCAAGTGTCGTTCACGGt                                                                              >  2:55370/1‑251 (MQ=255)
                                                                                                                                                                  cATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCCCTGCAAGTGTCGTTCACGGTTGGc                                                                          >  2:214582/1‑251 (MQ=255)
                                                                                                                                                                   agcAGCTGATGACTCCAGACAGCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCCCATTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGTTAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCg                                                                         <  2:151694/249‑1 (MQ=255)
                                                                                                                                                                   aTCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCg                                                                         <  1:526809/251‑1 (MQ=255)
                                                                                                                                                                          gacgacGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATACCTTGTAATTCTAATAATTTCAGGGCTTCTCACAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCgcgc                                                                  >  1:651635/1‑251 (MQ=255)
                                                                                                                                                                          gacgacGCCAGACATCTATCCGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTGGTAATTCTAATAAGTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACTCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATcgcg                                                                   <  1:78806/250‑1 (MQ=255)
                                                                                                                                                                                gCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAAtgt                                                            <  1:470105/251‑1 (MQ=255)
                                                                                                                                                                                 ccAGACATCTATCTGTTGGCTAAGATAACGCAGGGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAAt                                                                                                                                                                               >  1:572560/1‑135 (MQ=255)
                                                                                                                                                                                 ccAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCTTAGAGGATTGCCTCTTTTTCTGGCTAGAAGCGGCGTAGCGCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAAt                                                                                                                                                                               <  2:572560/135‑1 (MQ=255)
                                                                                                                                                                                       atctatctGTTGGCTAAGATAACGCAGCGCATCGTAGCGAATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCTGAAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTTTTTTCCTCCGTCTTTCCCACTGCAAGTTTCGTTCACGGTTGGCGATCGCGACAATGTGCGCtt                                                      >  2:543270/1‑249 (MQ=255)
                                                                                                                                                                                         ctatctGTTGGCTAAGATATCGCAGCGCATCGTAGAGGATTGCCTCTTTTTCTGGCCAGAAGCGTCGTAGTTAGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCATCGTTGTGTTGGCaac                                                                                                                                                                               <  2:517017/127‑2 (MQ=255)
                                                                                                                                                                                         ctatctGTTGGCAAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGAGTTGGCAAt                                                                                                                                                                               >  1:517017/1‑127 (MQ=255)
                                                                                                                                                                                          tatctGTTGGCTAAGATAACGCAGCGCAACGTATATGGTTGTCTCTTTTTCTGGCTAGAAGCGGCGTAGTGCGTCCATGGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAAGCCCTTGTAATTCTAATAAGTTCAGGGCTTCTCCCAGTACATCTGCACGTTGCACGCTAGTTTCCTCCTTCGAGCCCACTACAAGTGGCGTTCACGGGTGGTGGTCGCTCAAATGTGTGCTGAAg                                                   <  2:673764/250‑1 (MQ=255)
                                                                                                                                                                                           atctGTTTGCTAAGAGAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAgg                                                  <  2:7083/250‑1 (MQ=255)
                                                                                                                                                                                                 tGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAg                                            <  1:622014/250‑1 (MQ=255)
                                                                                                                                                                                                 tGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTAGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCACTCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCCCGTTGCAAGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGc                                           >  2:402069/1‑251 (MQ=255)
                                                                                                                                                                                                       aaGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATc                                        <  1:181887/249‑1 (MQ=255)
                                                                                                                                                                                                              cGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCAAGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTAACGGTTGGCGTTTGCGCAAATGTGCGCTGTAGGTTTCAGCATCCATAAAGcc                               >  2:605286/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGAcgcgg                     >  2:118315/1‑250 (MQ=255)
                                                                                                                                                                                                                          tAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGGTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACCTCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCCAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCAGACGGTTTCAGCATCCATAAAGCCCGTGAcgcggg                    >  2:194000/1‑248 (MQ=255)
                                                                                                                                                                                                                            gagGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCtt                 >  1:114480/1‑251 (MQ=255)
                                                                                                                                                                                                                             agGATTGCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTtg                >  1:17692/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCCTCTTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCCCGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATCAAGCCCGGGACGCGTGCTTGTggctt          >  2:196901/1‑248 (MQ=255)
                                                                                                                                                                                                                                    ccgctTTGTCTGGCTAGCAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAAACATCTACAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGTAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGaggc         <  2:560073/248‑4 (MQ=255)
                                                                                                                                                                                                                                     ctctTTGTCTGGCTAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCt        >  1:194990/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gTCTGGCTAGAAGCGGCGTAGGTCGTCCAGTGGATATTCCACTCGTTCAGCAAACATCTCCAGCGTGGTGTTTGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTGCCGCCGTCTTTCCCACTGCAAGTGGCGTTCACGGGTGGCGATCGCGCAAATGTGCGCTGAAGGTTGCAGCAGTCATAAAGCCCGTGACGCGTGCTTGTGTATGTTCCTGGc  <  2:225755/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
TTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGC  >  NC_000913/4362943‑4363427

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