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

TTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGC  >  NC_000913/4362943‑4363427
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
tttgtGCTGATCGGCCAGGTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCATGCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACAGCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAgcg                                                                                                                                                                                                                                               <  2:543935/245‑1 (MQ=255)
         aTCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAAAGGATTGCCTCTTTGTCTTGCTAGAAGAGGCGT‑AGCTCGt                                                                                                                                                                                                                                   >  1:340448/1‑251 (MQ=255)
                          aaGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCcaca                                                                                                                                                                                                                  <  1:508475/250‑1 (MQ=255)
                            gggTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACg                                                                                                                                                                                                                >  1:117872/1‑251 (MQ=255)
                            gggTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACg                                                                                                                                                                                                                <  1:637843/251‑1 (MQ=255)
                                    aGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCa                                                                                                                                                                                                         >  1:560073/1‑250 (MQ=255)
                                      gATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCGTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAAcc                                                                                                                                                                                                      <  1:27455/251‑1 (MQ=255)
                                         cGGGATAAAACGTACAGGCAGCGATAAATCGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCATTGGCTAATCCCCACGTTCCGCAAccctc                                                                                                                                                                                                   >  1:220798/1‑248 (MQ=255)
                                            gATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCCTCCCGCAGACGTCGACGACGCCAGCCATCCATCTGGTGG‑CTAAGATAACGCAGCGCAGTGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGACGG‑AGACCATCCAGTGGAAAATCCACACGTTCataaacacactac                                                                                                                                                                                                >  2:323881/1‑239 (MQ=255)
                                                      aCAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCCCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGTTAATCCACAAGTTCAGCAACCATCTCCAGCGGGgtgt                                                                                                                                                                                      >  1:151694/1‑251 (MQ=255)
                                                       cAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGAGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCCTCAGCTGACGACGACAGCCCTCTATCCGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTACCTCTTTGTCTGGCTTAAAGCGGCGT‑AGCTCGTCCAGGGGATAATCCACACGTTCCACAACCAGCTcccgccgggtgtt                                                                                                                                                                                     >  2:110131/1‑240 (MQ=255)
                                                          gCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  <  2:316400/251‑1 (MQ=255)
                                                           cAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  >  1:321650/1‑250 (MQ=255)
                                                                aTAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTTTCTGTTGG‑CTAAGATAACGCAGCGCCACGTAGAGGATTTCCTCTTTGTTTGGCTAGAAAAGGCGG‑AGCTCGACCCGGGGGAAATTCACAAGTTTCGCCACCATCTCCAGCggggggtttaatagga                                                                                                                                                                             >  2:576026/1‑235 (MQ=255)
                                                                aTAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGtt                                                                                                                                                                                     <  1:576026/242‑1 (MQ=255)
                                                                          aGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCCGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGACGAGGATTGCCTCTTTGTCTGGCTTGAAACGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCt                                                                                                                                                                  >  1:673764/1‑251 (MQ=255)
                                                                            cccGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTTCGTTTCGTCCAGCAACAGCTGACGACGCCCAACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGAATGCCTCTTTTTCTGGCTAGAAGCGGCGTGAG‑TCGGCCAGTGGGTAATACCCACGGTTAGCAACCATCTCCAGCGTGGGGTTGGCAATCCCTTGTAATTCtt                                                                                                                                                                 >  2:659017/1‑249 (MQ=255)
                                                                             ccGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTGAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGGTAATCCACACGTTCAGCAACCATCTCCAGCGTGGGGTTGGCAAACCCTTGTAATTCtaat                                                                                                                                                               >  1:501213/1‑251 (MQ=255)
                                                                              cGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAACGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTAAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCtaat                                                                                                                                                               >  1:567780/1‑250 (MQ=255)
                                                                                gATAGCGGGTGTTTTTAACGACCTCCGATAACGCCTGATAACGTGCCAGCAGCTTGTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCGGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTATAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTaataat                                                                                                                                                            <  2:451049/251‑1 (MQ=255)
                                                                                     cgggtgttTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAg                                                                                                                                                       <  1:170822/247‑1 (MQ=255)
                                                                                     cGGTTGTTTTGAACGCACTCCTATAACGCCTGATTACGTGCCTGCAGTTTTTGGTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAg                                                                                                                                                       <  2:594919/251‑1 (MQ=255)
                                                                                       gttgttTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  >  1:381776/1‑222 (MQ=255)
                                                                                       gttgttTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGc                                                                                                                                                                                  <  2:381776/222‑1 (MQ=255)
                                                                                                aaCGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTcccc                                                                                                                                            >  1:569465/1‑250 (MQ=255)
                                                                                                     tcTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCCGACATCTATCTGTTGG‑CTAAGATTACGCAACGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAAGGGCGT‑AGCTCGTCCCGTGGATAATCCACACGTTCAGAAACCATCTCCAGCGTGGGGTTGGCAATCCTTTGTAATTCTAATAATTTCAGGGCTTCTCCCAggac                                                                                                                                        >  2:428717/2‑247 (MQ=255)
                                                                                                          gATAACGCCTGATAACGTGCCATCAGCTTTTGTGCGGCGGTTTGCGTTTCGTCCAGCATCAGCCGACGACGCCAGACAGCTATCTGTTGG‑CTGAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATCTCAGGGCTTCTCCCAGTACATCtgc                                                                                                                                  <  2:501213/251‑3 (MQ=255)
                                                                                                             aaCGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACg                                                                                                                               >  1:313492/1‑251 (MQ=255)
                                                                                                               cGCCTGATAACGTGCCAGCAGCTTTGTTTCGGCGGTTCGCGTGTCGGCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGT‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑TGCTCGTCCAGTGGATAATCCTCACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGtt                                                                                                                             <  2:55724/251‑1 (MQ=255)
                                                                                                                gCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTg                                                                                                                            >  1:320070/1‑251 (MQ=255)
                                                                                                                    gATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCAAACGTTCAGCCACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACg                                                                                                                        >  1:348782/1‑251 (MQ=255)
                                                                                                                    gATAACGTGCCAGCAGCTTTGGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCAc                                                                                                                         <  2:320070/250‑1 (MQ=255)
                                                                                                                        aCGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGGTTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCCGCAACCATCTCCAGCGTGGGGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTTCACGCTAt                                                                                                                    >  1:272165/1‑251 (MQ=255)
                                                                                                                          gTGCCAGCAGCTTT‑GTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTAtttt                                                                                                                 <  1:45502/251‑1 (MQ=255)
                                                                                                                                    cTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCgtggt                                                                                                                                                                                        >  1:92804/1‑171 (MQ=255)
                                                                                                                                    cTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGG‑AGCGCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCgtggt                                                                                                                                                                                        <  2:92804/171‑1 (MQ=255)
                                                                                                                                        tGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTc                                                                                                     <  2:313492/250‑1 (MQ=255)
                                                                                                                                        tGTTCGGCGGTTTGCGTTTCGGCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTc                                                                                                     <  2:569465/250‑1 (MQ=255)
                                                                                                                                          gtcggcgggttggtcgtttgccAGCAACTGCTGACGACGCCAGACAGCTATC‑GTTGGGCTAAGATATCTCATCGCAGCGTAGAGGATTGCTGCTTTGTTTGGCTAGAAGCGGCGT‑AGCTCGTTCAGTGGAGAAGCCACACGTTCAGCAATCATCTCCAGCGTGGTGTAGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTTCTGAATCACGGGGCACGCTATTTTTCGCCGTCTCTCCCa                                                                                                  <  2:256322/231‑1 (MQ=255)
                                                                                                                                                ggTTTGCGTTTCGTCCAGCATCAGCTGACGACGCAAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACCTCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCa                                                                                             >  2:470756/1‑250 (MQ=255)
                                                                                                                                                  tttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTCGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCAAGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAgt                                                                                          >  2:227306/1‑251 (MQ=255)
                                                                                                                                                   ttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTAGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTTCATCTTCACGTTGCACGCTATTTTCCTCAGTCTTTccccatgcacgtg                                                                                         >  2:285962/1‑241 (MQ=255)
                                                                                                                                                    tGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAgtgt                                                                                        >  1:541319/1‑251 (MQ=255)
                                                                                                                                                     gCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTc                                                                                       >  2:181887/1‑251 (MQ=255)
                                                                                                                                                      cGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCAACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCCCTGCAAGTGTCg                                                                                      >  2:622303/1‑251 (MQ=255)
                                                                                                                                                       gTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGAAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGA‑AGCTCGTCCAGTGGATAATCCACACGTTCCGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTCCCAGTGTCGt                                                                                     >  2:196909/1‑251 (MQ=255)
                                                                                                                                                            gTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACg                                                                                >  2:308215/1‑251 (MQ=255)
                                                                                                                                                              ccAGCATCCGCTGACGACGCCAGACATCTATCTGTTGG‑CTAACATACCGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAACGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATACCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCACACTGCAAGTGTCGTTCACGGt                                                                              >  2:55370/1‑251 (MQ=255)
                                                                                                                                                                  cATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCCCTGCAAGTGTCGTTCACGGTTGGc                                                                          >  2:214582/1‑251 (MQ=255)
                                                                                                                                                                   agcAGCTGATGACTCCAGACAGCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCCCATTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGTTAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCg                                                                         <  2:151694/249‑1 (MQ=255)
                                                                                                                                                                   aTCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCg                                                                         <  1:526809/251‑1 (MQ=255)
                                                                                                                                                                          gacgacGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATACCTTGTAATTCTAATAATTTCAGGGCTTCTCACAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCgcgc                                                                  >  1:651635/1‑251 (MQ=255)
                                                                                                                                                                          gacgacGCCAGACATCTATCCGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTGGTAATTCTAATAAGTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACTCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATcgcg                                                                   <  1:78806/250‑1 (MQ=255)
                                                                                                                                                                                gCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAAtgt                                                            <  1:470105/251‑1 (MQ=255)
                                                                                                                                                                                 ccAGACATCTATCTGTTGG‑CTAAGATAACGCAGGGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAAt                                                                                                                                                                               >  1:572560/1‑135 (MQ=255)
                                                                                                                                                                                 ccAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCTTAGAGGATTGCCTCTTTTTCTGGCTAGAAGCGGCGT‑AGCGCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAAt                                                                                                                                                                               <  2:572560/135‑1 (MQ=255)
                                                                                                                                                                                       atctatctGTTGG‑CTAAGATAACGCAGCGCATCGTAGCGAATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCTGAAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTTTTTTCCTCCGTCTTTCCCACTGCAAGTTTCGTTCACGGTTGGCGATCGCGACAATGTGCGCtt                                                      >  2:543270/1‑249 (MQ=255)
                                                                                                                                                                                         ctatctGTTGG‑CTAAGATATCGCAGCGCATCGTAGAGGATTGCCTCTTTTTCTGGCCAGAAGCGTCGT‑AGTTAGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCATCGTTGTGTTGGCaac                                                                                                                                                                               <  2:517017/127‑2 (MQ=255)
                                                                                                                                                                                         ctatctGTTGG‑CAAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGAGTTGGCAAt                                                                                                                                                                               >  1:517017/1‑127 (MQ=255)
                                                                                                                                                                                          tatctGTTGG‑CTAAGATAACGCAGCGCAACGTATATGGTTGTCTCTTTTTCTGGCTAGAAGCGGCGT‑AGTGCGTCCATGGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAAGCCCTTGTAATTCTAATAAGTTCAGGGCTTCTCCCAGTACATCTGCACGTTGCACGCTAGTTTCCTCCTTCGAGCCCACTACAAGTGGCGTTCACGGGTGGTGGTCGCTCAAATGTGTGCTGAAg                                                   <  2:673764/250‑1 (MQ=255)
                                                                                                                                                                                           atctGTTTG‑CTAAGAGAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAgg                                                  <  2:7083/250‑1 (MQ=255)
                                                                                                                                                                                                 tGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAg                                            <  1:622014/250‑1 (MQ=255)
                                                                                                                                                                                                 tGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTAGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCACTCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCCCGTTGCAAGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGc                                           >  2:402069/1‑251 (MQ=255)
                                                                                                                                                                                                       aaGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATc                                        <  1:181887/249‑1 (MQ=255)
                                                                                                                                                                                                              cGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCAAGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTAACGGTTGGCGTTTGCGCAAATGTGCGCTGTAGGTTTCAGCATCCATAAAGcc                               >  2:605286/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGTAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGAcgcgg                     >  2:118315/1‑250 (MQ=255)
                                                                                                                                                                                                                          tAGAGGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGGTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACCTCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCCAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCAGACGGTTTCAGCATCCATAAAGCCCGTGAcgcggg                    >  2:194000/1‑248 (MQ=255)
                                                                                                                                                                                                                            gagGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCtt                 >  1:114480/1‑251 (MQ=255)
                                                                                                                                                                                                                             agGATTGCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTtg                >  1:17692/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCCTCTTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCCCGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATCAAGCCCGGGACGCGTGCTTGTggctt          >  2:196901/1‑248 (MQ=255)
                                                                                                                                                                                                                                    ccgctTTGTCTGGCTAGCAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAAACATCTACAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGTAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGaggc         <  2:560073/248‑4 (MQ=255)
                                                                                                                                                                                                                                     ctctTTGTCTGGCTAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCt        >  1:194990/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gTCTGGCTAGAAGCGGCGT‑AGGTCGTCCAGTGGATATTCCACTCGTTCAGCAAACATCTCCAGCGTGGTGTTTGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTGCCGCCGTCTTTCCCACTGCAAGTGGCGTTCACGGGTGGCGATCGCGCAAATGTGCGCTGAAGGTTGCAGCAGTCATAAAGCCCGTGACGCGTGCTTGTGTATGTTCCTGGc  <  2:225755/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
TTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGATAACGCCTGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACATCTATCTGTTGG‑CTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGT‑AGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGC  >  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: ‑