Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,326,039 A→T 19.7% M432L (ATG→TTG)  yciQ → enhancer of membrane protein expression, putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,326,0390AT19.7% 129.9 / 2.4 51M432L (ATG→TTG) yciQenhancer of membrane protein expression, putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (5/35);  new base T (9/2);  total (14/37)
Fisher's exact test for biased strand distribution p-value = 2.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGTTTGCCTGACCGGTTATCTCACTGGCA  >  NC_000913/1325803‑1326267
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
tGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTAAGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                                                                                     <  2:743070/239‑1 (MQ=255)
   aaGCGTTATTTCGTAAAGGCGTTAAGGTCCTTCTTACGGGGAGCCGCAACAGGGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGAt                                                                                                                                                                                                                                   <  1:611679/238‑1 (MQ=255)
   aaGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTCGGTAACTTTGGTCCCGTGGGCTGGGGA‑GGAATGGGGCTGGTGGGCGGTATGTTTTTTTTGAtt                                                                                                                                                                                                                                  >  2:481130/1‑238 (MQ=255)
                ttagggCTATAAGGTCGTTTTTACGGGGAGCCGCACCAGATTTTTACGCAGGGTATTTTTACGGATGCGGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTTCCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTc                                                                                                                                                                                                                    <  1:774096/235‑1 (MQ=255)
                 aaaGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTc                                                                                                                                                                                                                    <  2:140047/238‑1 (MQ=255)
                      cGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCtttt                                                                                                                                                                                                               <  2:793162/238‑1 (MQ=255)
                               cGTTCTTACGGGGAGACCCAACAGAGTCTTAAGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGtt                                                                                                                                                                                                     <  1:118128/239‑1 (MQ=255)
                                gTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGGCATTCTTTATGGTAAATTGAGTCCTGTGGTCTGGGCG‑GGGGTGAGTGTTGTGGGGGGTATATTTTTTTTTTTCTTCTG‑GGTTTTTCCtttttttttttttt                                                                                                                                                                                                     >  1:14169/1‑238 (MQ=255)
                                 ttCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGTCCTGTAGGCTGGGCC‑GGAATTGGTTTGGGTGGGGATATGTTTTTTTTGGTTTTGTG‑GGTTTTTCCTTTTTTTTTTTTTTc                                                                                                                                                                                                    >  2:805379/1‑238 (MQ=255)
                                        cggggGACGCAACTGATTCTTTCGCAGGGCTTTTTTACGGGTGCAGAAATTTTAGCAGCCGCGTAAAAAGTTTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGc                                                                                                                                                                                            <  1:333111/238‑1 (MQ=255)
                                            agaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTTTTTACGGTAAATTTAGTCCTGTAGGGTTGGGC‑GGGATGAGGCTGGTGGGGGGGATGTTTTTTTTGTATTGCGT‑GGTTATTTCTTTTTTTTTTTTTTTCCtttgggtttt                                                                                                                                                                                        >  1:300537/1‑231 (MQ=255)
                                              aCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGCCATTCTTTTCGGTAACTTGAGTCCTGTTGGTTGGGCA‑GGGAATAGGCTTGTGGGGGGTATGTTTTTTTTGAACTGCCG‑GATTTTTCCTTTTTTTTTGTGTTTCCtttggacttttt                                                                                                                                                                                      >  2:154505/1‑229 (MQ=255)
                                                   accgattttttccccgggTATTTTTACGGAGGCAGAATTTTTATCTGCCGCGTAAAAAGTTTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGcc                                                                                                                                                                                 <  1:481130/224‑1 (MQ=255)
                                                         tCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGCCATTCTTTACGGGAACTTTGGTCCTGTAGGCTGGGCA‑GGGATTGGTCTGGGGGGCGGTATGTTTTTTTTGGTCTGCTG‑GATTTTTTCTTTTTTTTTTTTTTCCCTTTAGGTTTTGTTTGTCCCCgat                                                                                                                                                                           >  2:612765/1‑239 (MQ=255)
                                                            tacgaggggCATTTTTACGGATGCAGAAATTTTAATTGCCGGGTAAAAAGTCTTCGGTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCgatgat                                                                                                                                                                        <  2:705112/233‑1 (MQ=255)
                                                            tACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTTTTTTTGGTAACTTTGGTTCCGTAGGGCGGGCC‑GGGATGAGTCTGGTGGGCGGTCTGTTTATTTTGTTTGGCTG‑GGTTTTTCCTTTTTTTATTTTTTCCCTTTTAGTTTTGTTTGtcccggggggg                                                                                                                                                                        >  1:423365/1‑231 (MQ=255)
                                                                   ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCCTTTTTTTTGGTAACTTGGGTCCTGTTGGCTGGGCG‑GGGATGGGGCTGGGGGGCGGATTGTTTTTTTTTGTTTGGTG‑GGTTTTTCCTTTTTTTTTTTGTTTCCTTTTGCttttttttttcctcgcgtgtgtaaaa                                                                                                                                                                  >  1:550196/1‑220 (MQ=255)
                                                                   ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGGAAATTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTTCCTTTAGCTTTTGTTTGCCCCCGATGATGGCaaa                                                                                                                                                                  >  2:15690/1‑237 (MQ=255)
                                                                   ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGATTTACCGGTCTGATACTTTTTTGCAATGGGGGGGACTGGCAATAATGTCGGCTATTCTTTTCTGGTACTTGGGTCACGTTGGCTTGGGC‑GGGCTGGGTCTGGGTGGGGGGATGTTTTTTTTTATATCCTGTGTTTTTTCTCTTTTTTTTTTTTCCCCATGGGCTTTTGTTTTccccgcgttttgtcca                                                                                                                                                                   >  1:738162/1‑225 (MQ=255)
                                                                    gttttttttCGGATGGAGAATTTTTATCTGCGGTGTTAAAAGTTTTCGTTTTACCGGTTTGATTTTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGc                                                                                                                                                                 <  1:424254/235‑1 (MQ=255)
                                                                    gCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTGTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGc                                                                                                                                                                 <  2:364897/238‑1 (MQ=255)
                                                                      ttttttAGGGATGCAGAAGTTTTATCTGCCGCGTAAAAAGTTTTGGTTTTACAGGTCTGATTCTTTTTTGCAATGGGGTGGACGGGCAATTGTGGCGGTCATTCTTTACGGTAACTTGAGTACTGTATGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTGTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCt                                                                                                                                                               <  1:119589/237‑1 (MQ=255)
                                                                            tcgtattaaGAAATTTTATTGGTCGGGTAAAAAGTTTTCGTTTTACCGGTCTGATGCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTa                                                                                                                                                        <  1:173307/231‑1 (MQ=255)
                                                                                  gacgAAATTTTTTCTGGCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCgt                                                                                                                                                   <  1:8777/235‑1 (MQ=255)
                                                                                             agcTGCCGGGTAAAAAGGTTTCGTTTTACCGGTGTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACt                                                                                                                                       <  2:748131/237‑1 (MQ=255)
                                                                                             aTCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACt                                                                                                                                       <  2:413676/239‑1 (MQ=255)
                                                                                                              tCTTCGTTTTACCGGTTTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGAt                                                                                                                      <  2:529134/239‑1 (MQ=255)
                                                                                                               cTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGAt                                                                                                                      <  2:274124/238‑1 (MQ=255)
                                                                                                                           ggTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTg                                                                                                         <  2:130515/239‑1 (MQ=255)
                                                                                                                              cttgttcTTTTTGGCAATGGGGTTGACTGGCAATTTTGGCGGTCCTTCTTTACGGTAACTTGAGTACTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTAGTATGATCTGCTT‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGc                                                                                                       <  2:632583/232‑1 (MQ=255)
                                                                                                                                       ttttGCAATGGGGTGGACTGGTAATATTGGGGGTCATTCTTTGCGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATtct                                                                                             <  1:540469/239‑1 (MQ=255)
                                                                                                                                         ttGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTCACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTTGTGGGCGATCTGTTTATTATGATCTGGTG‑GGTTATTCCTTTTTTTTTTTTTTTCCCTTAGCTTTTTTTTGTCCCCGGTGATTACAAGCCTTTGGTTTATCCTGTTTTTATTCCTTTTTTTTTTCCACAGTTTTTTTTAGGCGGGGGCttttttttcc                                                                                            >  1:868295/1‑230 (MQ=255)
                                                                                                                                                ggggTGGACTGGAAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATc                                                                                    <  2:14169/239‑1 (MQ=255)
                                                                                                                                                        cTGGAAATTTTGGGGGTTATTCTTTTCGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTTTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCg                                                                             <  1:35320/238‑1 (MQ=255)
                                                                                                                                                         tGGGAATTTTGGCGGTCAATCTTTACGGTAACTTGAGTCCTGTAGGGTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCgg                                                                            <  2:423365/238‑1 (MQ=255)
                                                                                                                                                             aaTTTTGGGGTTTATTTTTTTCGGTAATTTGAGGCCTTTAGGGTGGGCA‑GGAATGGGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATg                                                                       <  1:780355/239‑1 (MQ=255)
                                                                                                                                                              atatTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATg                                                                       <  2:868295/238‑1 (MQ=255)
                                                                                                                                                                ttttgCGGTCTTTTTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAt                                                                    <  1:372410/235‑1 (MQ=255)
                                                                                                                                                                atTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAt                                                                    <  2:738162/239‑1 (MQ=255)
                                                                                                                                                                       gttaTTTTTTTCGTTAACTTGTTTCCTGTAGGCTGGGCA‑GGAATGATTCTGGTGGGCGATATGTTTATTATGATCTTCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCtttt                                                              <  1:805379/235‑1 (MQ=255)
                                                                                                                                                                        tCATTGTTTACGGTAACGTTAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTa                                                             <  2:216897/238‑1 (MQ=255)
                                                                                                                                                                         taTTGTTTGCGGTAACTTGAGTCCTGTAGGCTGGGCA‑TGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTAc                                                            <  2:468653/237‑1 (MQ=255)
                                                                                                                                                                         cATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTAc                                                            <  2:300537/238‑1 (MQ=255)
                                                                                                                                                                                         tttggTCCTGTTGGCTTGGCAGGGATTGTGT‑TGGGGGGGGTGTTGTGTATTTTTATCTGATG‑GATTAGTCCTTTTTTATGTTTTTCCCTTGAGCTTGTGTTTGCCCGCGATGATGACAAGCCTTGAGTTAAGCGTTTAATCATCACTTTGTTTCTACCACTGATTTGTTCAGGCGTGGCCGTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCggg                                           <  1:795836/235‑1 (MQ=255)
                                                                                                                                                                                         ttGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCggg                                           <  2:550196/239‑1 (MQ=255)
                                                                                                                                                                                             gTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTAt                                       >  2:289630/1‑239 (MQ=255)
                                                                                                                                                                                                          ggCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGtt                           >  2:157475/1‑238 (MQ=255)
                                                                                                                                                                                                                 aaTGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGTTTGCCTg                     >  2:813862/1‑238 (MQ=255)
                                                                                                                                                                                                                         tggtggGCGATATGTTTATTATGATCCGCTG‑GATTATTCCTTTTTTATTTTGTTCCCATGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTCCCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGTTTGCCTGACCGGTTAt            <  2:60435/239‑1 (MQ=255)
                                                                                                                                                                                                                             gggCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGTTTGCCTGACCGGTTATCTCa        <  2:69825/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tttGTTTTTTTTGTTCTGTGT‑GATTATTCCTTTTTTATTTTTTTCCCTTGAGCTTTGGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGTTTGCCTGACCGGTTATCTCACTGGCa  <  1:612765/237‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
TGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCA‑GGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTG‑GATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTATTCTTTTACTGGTATATGCCAGCGGGTTATTTTACCGCTGTTTGCCTGACCGGTTATCTCACTGGCA  >  NC_000913/1325803‑1326267

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