Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,270,497 C→T noncoding (96/377 nt) rnpB ← RNase P catalytic RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,4970CT96.1% 144.0 / NA 51noncoding (96/377 nt)rnpBRNase P catalytic RNA component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/27);  minor base G (0/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAA  >  NC_000913/3270261‑3270738
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
catCTAGGCCACCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTGACCGCCCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCAACCTGccc                                                                                                                                                                                                                                        <  2:368559/251‑1 (MQ=255)
          atcAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGccc                                                                                                                                                                                                                              <  1:303400/249‑1 (MQ=255)
           gCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTGTCTCTTTTGCCCTGGTCGGGGG‑GTTCCCCCCCCGGCGTTATCTTGCAACCCGCCCTCTGGAGcccc                                                                                                                                                                                                                             >  1:49109/1‑250 (MQ=255)
                            cAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGTCGGGGGGTTTCCCCCCC‑CGCGGTTTCTGGCACCCCGCCCCATTGAGCCCCGACCTTtctcccctccc                                                                                                                                                                                                            >  1:323498/1‑250 (MQ=255)
                             aaGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACTGGGCGTGGGGTTTCCCCCCC‑GGCGTTATTTGGGCCCCTGCCCttt                                                                                                                                                                                                                                     >  1:456725/1‑223 (MQ=255)
                             aaGCAGCCTACCCGGGTGCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCAGCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTAt                                                                                                                                                                                                                                     <  2:456725/225‑1 (MQ=255)
                                        ccGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTCTCTTTTGCTCTGGTGGTGGGGTTTCCCCCCCCGCCGTT‑TCTGGCACCCTGCCCCCTGGGGCCCCGGCTTTCCTCCCCCCCCCCCCCCTccccc                                                                                                                                                                                                >  1:101081/1‑251 (MQ=255)
                                            gTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGATTGCTCTCTGTTTCACTGGGCGTGGG‑GTTCCCCCCCCGGCGGTATCTGGCACCCCGCCCTCTGGGGCCCGGGCTTTCCCCCCCCCCCCCCCTCTccccccc                                                                                                                                                                                              >  1:193141/1‑247 (MQ=255)
                                              tCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑aga                                                                                                                                                                                          <  1:213790/250‑1 (MQ=255)
                                                     gggCCGTACCTTATGACCCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCCCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGgacgacg                                                                                                                                                                                  <  2:323498/251‑1 (MQ=255)
                                                           tACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTTTTGTCCTGGGCGTGGGGTTTCCCCCCC‑GGCGTTATCTGTCACCCTGCCCTCTGGAGGCCGGACTTTTCTCCCCCCCTCCCCGCTCCCCCCC‑AGAGGAccacgccga                                                                                                                                                                              >  1:356786/1‑249 (MQ=255)
                                                             ccTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTTGTCGTGGGGTTTCCCCCCCAGG‑GGTATCTTGCACCCCGCCCCATGGCGCCCCGACTTTCCCCCCCTCCCCCCCCCTTCCCCCC‑AGAGGACGACGACCAAgcc                                                                                                                                                                          >  1:389419/1‑250 (MQ=255)
                                                                    aaCCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGGGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCGTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTAGCGTGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCTTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑ACGACTg                                                                                                                                                                  <  2:62550/251‑1 (MQ=255)
                                                                        ccTATTGGGCCGTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACTCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCGGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCTTCTCCCCCGA‑AGAGGACGACGACGAAGCT‑GCGACTGTCTg                                                                                                                                                              <  2:359420/251‑1 (MQ=255)
                                                                           aTTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTgg                                                                                                                                                             <  1:504327/249‑1 (MQ=255)
                                                                                     gcaccGGGTGGAGTTTACCGTGCCCCGGACTGTTACCCGCCTCGCGGTGCGCTCTTACCGCACCCTTTCACCCTGACCTGATCCCGCTTGCTCTGGCCATCGGCGGTTTGCTCTCTGTGGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCATCCTGCCCTATGGAGCCCGGACTTTTCGCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGGCGAAGCG‑GCGGCTGTCTTATCAGATTCGgcg                                                                                                                                                 <  2:569475/248‑1 (MQ=255)
                                                                                          gggTGGAGTTTACCGTGCCTCGGACTGTTACCAGCCTCGCGGTGCGCTCTTACCGCCCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAg                                                                                                                                            <  2:519785/251‑1 (MQ=255)
                                                                                                  tttACCGTGCCACGGACTGTTACCAGCCGCGAGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCCCTGGTCGGGGG‑TTTTCCCCCCCGGCGTTATCTTGCACCCTGCCCTCTTGAGCCCGGACTTTCCTCCCCCCCCCCCCCCTCCCCCCC‑CGAGGACGACGCCGAAACG‑GCGACTGTCTGGTCAGCCTCGGCGCGCACTATAGAgg                                                                                                                                    >  1:148330/1‑251 (MQ=255)
                                                                                                   ttACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTTCTCTCTTTTGCACTGGCCGTGGG‑GTTCCCCCCCCGGCCTTATCCGGCACCCTGCCCTATAGAGCCCGGACTTTCCTCCCCTCCCCCCGTCTCCCCCGC‑AGAGGACGACCACGCAACG‑GCGACTGTCTTGGCAACGTCGGCGGGCAGGATAGAggg                                                                                                                                   >  1:343966/1‑251 (MQ=255)
                                                                                                     aCCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCCCCGGTCGTGGG‑GTTCCCCCCCCGGCGGTTTCTGGCGCCCTGCCCCCTGGAGGCCGGGCCTTCCCCCCCCCCCCCCCGCTCCCCCCACAG‑GGGCGACGACGGAGCC‑GCGGCCAGCTGGTCGGCTTCCGGGCGCAGTCTAGAggggt                                                                                                                                 >  1:90806/1‑249 (MQ=255)
                                                                                                                         cAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCcg                                                                                                             <  1:369049/251‑1 (MQ=255)
                                                                                                                            ccGCGCGGGTCGGTCTTACCGCCCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCGCTCTGTTGCACTGGGCGGGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACGTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGt                                                                                                          <  2:194/251‑1 (MQ=255)
                                                                                                                             cgcgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTTAACTGGTCGTGGG‑TTTTACCCCCCGGCGGTTTCTGGCCCCCCGCCCCCTTGCGGCCGGGCTTTCCCCCCCCCCGCCCGTCTCCCCCCC‑AGAGGACGCAGAAGCA‑CGAGCGGCTGTCTGGTCAACCTCGGCGCGCCGGATAGAGGGTGTTCGCGcccctggacccccaccgg                                                                                                         >  1:399794/1‑236 (MQ=255)
                                                                                                                              gcgcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTg                                                                                                        <  1:122966/251‑1 (MQ=255)
                                                                                                                                   gTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCAt                                                                                                   <  2:101081/251‑1 (MQ=255)
                                                                                                                                               cGCACCCTTTCCCCCTTACCCGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAg                                                                                       <  2:8804/251‑1 (MQ=255)
                                                                                                                                                 cACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAgtg                                                                                     <  1:502280/251‑1 (MQ=255)
                                                                                                                                                  ccccTTTCACCCTTGCCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAgtgt                                                                                    <  2:90806/250‑1 (MQ=255)
                                                                                                                                                           cccTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGat                                                                           <  1:305784/251‑1 (MQ=255)
                                                                                                                                                                    gATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACGGGTCGTGGG‑TTTCCCCCCCAGGCGGTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACggga                                                                  <  2:468844/251‑4 (MQ=255)
                                                                                                                                                                      ccccGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAg                                                                <  1:408304/250‑1 (MQ=255)
                                                                                                                                                                                cctggCCCTCGGCGTTTTTCTCTCTGTTGCTCTGGTCGTGTG‑TTTCCCCCCCAGGCGTTATCTGGCCCCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCTGCGCTCAGTATAGAGTTTTTGCGCGCCCTTGTCTCCCCGCGTTGCGCTTTCCAATCGCCATTGTCGCGGCGATATTACGTGaatctcgtcagat                                                      <  2:551527/248‑12 (MQ=255)
                                                                                                                                                                                   ggCCATCGGCGGTTTGCCCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATAtt                                                  <  1:482887/251‑1 (MQ=255)
                                                                                                                                                                                          ggcggTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCAGTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATaa                                           <  2:49109/251‑1 (MQ=255)
                                                                                                                                                                                            cggTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGc                                         >  1:31920/1‑251 (MQ=255)
                                                                                                                                                                                             ggTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATACGcc                                        >  2:37554/1‑251 (MQ=255)
                                                                                                                                                                                             ggTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCACCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGcc                                        >  1:493430/1‑251 (MQ=255)
                                                                                                                                                                                                           tGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCtt                           <  1:72572/250‑1 (MQ=255)
                                                                                                                                                                                                                  gTCGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAAc                   <  2:493430/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAATCCCCGCGGAATGCTTCGTCCAACg                  <  1:118367/250‑1 (MQ=255)
                                                                                                                                                                                                                     gTGGG‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATACGCCCCGCGGAATGCTTCGTCCAACGt                 >  1:190737/1‑250 (MQ=255)
                                                                                                                                                                                                                       ggg‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTAcc              >  1:320544/1‑251 (MQ=255)
                                                                                                                                                                                                                         g‑TTTCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCt             >  1:342015/1‑250 (MQ=255)
                                                                                                                                                                                                                         g‑TTTCCCCCCCAGGCGTTAGCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTa            <  1:72907/251‑1 (MQ=255)
                                                                                                                                                                                                                            ttCCCCCCCAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATg          >  1:228862/1‑251 (MQ=255)
                                                                                                                                                                                                                             tCCCCCCCAGGAGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATCGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCAAGTGTCGCGGCGCTATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGc         >  2:223995/1‑251 (MQ=255)
                                                                                                                                                                                                                              cccccccAGGTGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑TCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCTATTCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCTCGCGGAATGCTTCGTCCAACGTACCTATGc         <  2:406072/250‑1 (MQ=255)
                                                                                                                                                                                                                              cccccccAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCt        <  2:559217/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccccAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGt     >  1:565297/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccAGGCGTTATCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCCGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTc    >  2:344824/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cAGGCGTTATCCGGCACGCTGCCCGATGGAGCCCAGAATTACCTCCCCTCCTCCCGACTCCCCCGA‑ACAGGACGACGACTTAGCA‑GCCACTGCCTCGTCAGCTTCAGAGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCATCGTTGCGCATTCCAATCGCCGGTGGCGCGACGCTATCACGAGACGC‑CGTGCCGCTATTGTCATAATCACCGCGAAATGCTACGTccgacataactctactgctcac  >  2:190852/1‑231 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
CATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCTCCCCTCCGCCCGTCTCCCCCGA‑AGAGGACGACGACGAAGCG‑GCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACG‑GCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCTATGCTGGTCAA  >  NC_000913/3270261‑3270738

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