Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,923,235 C→T S149N (AGC→AAC)  rsmG ← 16S rRNA m(7)G527 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,2350CT96.2% 162.9 / ‑1.5 54S149N (AGC→AAC) rsmG16S rRNA m(7)G527 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (31/20);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAAC  >  NC_000913/3923002‑3923477
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
cTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTc                                                                                                                                                                                                                                      >  1:585169/1‑248 (MQ=255)
  ggTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCag                                                                                                                                                                                                                                    >  1:343639/1‑248 (MQ=255)
  ggTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACAACTGATTCGACCTGATATTCTTCGGGCAACAAAGAGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACAAGTTCACCATATCGTTCag                                                                                                                                                                                                                                    >  1:476942/1‑248 (MQ=255)
           tAACAGTCTAACCGGTCACTTTTTTATGATTTTTTTGATAAACATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTATCCCTTTCAGCGCGTAGAAACGTCTTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATAGCGttatgaaagtcaagagc                                                                                                                                                                                                                        >  2:344582/1‑236 (MQ=255)
                     aacGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATAACCACCAGATGACGTTAGCAATCGAGGGCTGGAATTTGAAGTTTAACCACTGATTAGACCTGATATTCTTCGGGCAACAAAGCGATTTAATCTTCCGGCATTTGCCCTTTAAGCGCGTAGAAACGGCCTTGCTCACCAGTAAGATGGTGTCACCAGTTCACCATAGCGTTCAGAGAGGCAAAAGCGCGGCTAAtt                                                                                                                                                                                                              >  2:247724/3‑251 (MQ=255)
                             aTTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTATTCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGTCACCAGTTCACCATATCGTTCAGAGCGGCAAAAGCGCGGCTACTTACGCCAt                                                                                                                                                                                                       >  2:169165/1‑250 (MQ=255)
                                             ttGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGa                                                                                                                                                                                       <  1:242898/250‑1 (MQ=255)
                                                 tAAAAATTAAATTTTATGTGCTTTAATCAACACCAGATGACGTTCGCCATCCTGGTCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGAGCAACATAGCGATTTCATCTTCCGCCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCCGGAAGATGGTGTCACCAGTTCACCATCTCGTTCTGAGAGGCAAAAGCGCGGCTAATTACGCCATCAACTGTCAGCTCTGAAGGa                                                                                                                                                                                   >  2:150362/1‑250 (MQ=255)
                                                  aaaaaTTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGaa                                                                                                                                                                                  <  2:343639/250‑1 (MQ=255)
                                                         aaattTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACtcttct                                                                                                                                                                          <  1:222654/251‑1 (MQ=255)
                                                             tttATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTAccc                                                                                                                                                                      >  2:3962/1‑251 (MQ=255)
                                                              ttATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTAccc                                                                                                                                                                      <  1:3962/250‑1 (MQ=255)
                                                               tATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCt                                                                                                                                                                     <  1:571041/250‑1 (MQ=255)
                                                               tATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTg                                                                                                                                                                    >  1:314402/1‑251 (MQ=255)
                                                                    gCTTTAAACACCACCAGATGCCGTTCGCCATCCAGGGCTGGAACCTGAAGGTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTGCGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGctct                                                                                                                                                                <  1:345572/250‑1 (MQ=255)
                                                                       ttAATCACCACCAGATGACGTTCGCCATCCCGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCCGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTGTACCCTGCTCTgaac                                                                                                                                                            >  2:637873/1‑248 (MQ=255)
                                                                          atccccaccaGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTgg                                                                                                                                                         <  1:602211/247‑1 (MQ=255)
                                                                          aTCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTgg                                                                                                                                                         >  2:344984/1‑251 (MQ=255)
                                                                            caccacCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCGAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGtt                                                                                                                                                       >  1:260153/1‑251 (MQ=255)
                                                                            caccacCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGtt                                                                                                                                                       >  1:60270/1‑251 (MQ=255)
                                                                            caccacCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGAATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGACGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGtt                                                                                                                                                       >  1:602941/1‑251 (MQ=255)
                                                                                    gATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATtc                                                                                                                                               >  1:634407/1‑251 (MQ=255)
                                                                                      tGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATtctc                                                                                                                                             >  1:609057/1‑251 (MQ=255)
                                                                                                cATCCAGGGCTGGAACCTGAAGTTTACCCACTGCTTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGc                                                                                                                                   >  1:235021/1‑251 (MQ=255)
                                                                                                   ccAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCTGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCa                                                                                                                                <  2:404174/251‑1 (MQ=255)
                                                                                                        gCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCCCCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTg                                                                                                                           >  1:219422/1‑251 (MQ=255)
                                                                                                             aaCCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCt                                                                                                                      <  1:344984/251‑1 (MQ=255)
                                                                                                                      tttAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAgg                                                                                                              >  2:236498/1‑250 (MQ=255)
                                                                                                                        tAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATGTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTGGCACCTGACGAAGGaa                                                                                                            <  2:643328/250‑1 (MQ=255)
                                                                                                                         aaCCACTGCTTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAac                                                                                                          >  2:52318/1‑251 (MQ=255)
                                                                                                                         aaCCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAac                                                                                                          <  2:479838/251‑1 (MQ=255)
                                                                                                                               tGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTCAGCTCATGTTGCACCTGACGAAGGAAacgcacgc                                                                                                    >  1:337708/1‑251 (MQ=255)
                                                                                                                               tGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAacgcacgc                                                                                                    >  2:105545/1‑251 (MQ=255)
                                                                                                                                gATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCAcgcg                                                                                                   <  1:90509/251‑1 (MQ=255)
                                                                                                                                gATTCGACCTGAGATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCGGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAacgcacg                                                                                                     <  1:52318/249‑1 (MQ=255)
                                                                                                                                    cgaacttatatTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAATCTCATGTTGCACCTGACGAAGGAAACGCACGCtttt                                                                                                <  2:402963/243‑1 (MQ=255)
                                                                                                                                    cGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGttt                                                                                                <  1:85615/250‑1 (MQ=255)
                                                                                                                                    cGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGttt                                                                                                <  1:236498/250‑1 (MQ=255)
                                                                                                                                        cTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCa                                                                                            >  1:127864/1‑250 (MQ=255)
                                                                                                                                              ttcttcGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTAt                                                                                     >  2:484890/1‑251 (MQ=255)
                                                                                                                                                    gggCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAAc                                                                                >  2:621694/1‑250 (MQ=255)
                                                                                                                                                           aaaGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGa                                                                         >  1:336148/1‑250 (MQ=255)
                                                                                                                                                                     cATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCGAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCa                                                              <  2:260153/251‑1 (MQ=255)
                                                                                                                                                                                      ccTTTCAGCGCGTAGAAACGGCCTTGCCAACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagag                                               >  2:280270/1‑249 (MQ=255)
                                                                                                                                                                                            aGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCGAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAAt                                       >  1:475829/1‑251 (MQ=255)
                                                                                                                                                                                            aGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGact                                       >  2:486322/1‑249 (MQ=255)
                                                                                                                                                                                                cgTATAAACGGCCGTGCTCCCCAGGAAGATGGTGTGACCAGGTCACCATATCCATCAGATAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCTAGTTTATGCTCATGTTGCACCTGACGATGGAGACGCACGCGTTTATCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCt                                   <  1:591578/251‑1 (MQ=255)
                                                                                                                                                                                                 gTANAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGGGGAATGCCTg                                  <  2:609057/251‑1 (MQ=255)
                                                                                                                                                                                                               gCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCATTCCTGGt                       <  2:649881/248‑1 (MQ=255)
                                                                                                                                                                                                               gCTCACCAGGAAGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGt                       >  1:649881/1‑248 (MQ=255)
                                                                                                                                                                                                                         aaGATGGTGGCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACAt          >  2:300681/1‑251 (MQ=255)
                                                                                                                                                                                                                                tggCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGCGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATaa    >  1:582123/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAGTGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTCCCGACATCGATaaa   <  1:320256/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCACCAGTTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAAc  <  2:336148/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
CTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAAC  >  NC_000913/3923002‑3923477

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