Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 18,169 T→C T227T (ACT→ACC)  nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,1690TC98.2% 164.5 / ‑4.5 54T227T (ACT→ACC) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/30);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

CGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTG  >  NC_000913/17927‑18414
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
cGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGTCGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTAGCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTGGGGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCgg                                                                                                                                                                                                                                                <  2:310817/250‑1 (MQ=255)
      ccGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCGTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTGGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTaa                                                                                                                                                                                                                                          <  2:471954/250‑1 (MQ=255)
          tAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGGGTTGCAATCGGGGGTTCTCGCAACCCTGGCGGGGGTGATTGTc                                                                                                                                                                                                                                     >  2:198793/1‑251 (MQ=255)
                tGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGGTGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGGATTGAATCTGTGTGGGGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGGAATTGGCgggttt                                                                                                                                                                                                                               >  2:54739/1‑247 (MQ=255)
                 gAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGcttc                                                                                                                                                                                                                               <  1:812/250‑1 (MQ=255)
                 gAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGcttc                                                                                                                                                                                                                               <  1:414017/250‑1 (MQ=255)
                    gATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCtttt                                                                                                                                                                                                                           >  2:373634/1‑250 (MQ=255)
                         ttttGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTc                                                                                                                                                                                                                        <  1:462238/249‑1 (MQ=255)
                          tttGATGGCTCTGGCTATTATCGACGATCTTGGTGCCAGCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGACTCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTGGTGGTGTGGACTGCGGTGTTGGAATCGGGGGTTCACGCAACCCTGGCggggggtatttgagggttctttattactt                                                                                                                                                                                                                     >  2:91719/1‑227 (MQ=255)
                                                gacgaTCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGCAAGAGAAGCATGGGCGttt                                                                                                                                                                                               >  1:49336/1‑250 (MQ=255)
                                                         ggggCCAACATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCagcga                                                                                                                                                                                       <  2:151093/250‑1 (MQ=255)
                                                           ggtcatcatTATCATTGCATTGTTCTACACTAATGACTGATCGATGGCCTTTCTGGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGCCTGCGGTGTTGAAATCGGGGGTTCACGCATCCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGTCGTTCTCCagcgaagc                                                                                                                                                                                    <  1:417213/248‑1 (MQ=255)
                                                                          cGCATTGTTCTGCACTAATGTCTTATCGAGGGCCTCTCTTGGCGTCGCGGCTGGAGCAAGTGCGGTACTCGTGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGGTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGTCTGGAGCAtgtg                                                                                                                                                                      <  2:367871/250‑1 (MQ=255)
                                                                             aTTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTAGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGGGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGGAATTGACAGATTCTTTATTCCTTTGGAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGAATGGAGCATGGGtta                                                                                                                                                                   >  2:449121/1‑249 (MQ=255)
                                                                                   cTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCGCCAGCGAAGCGACTGGGGCATGGGTTACACCCgg                                                                                                                                                            >  1:104020/1‑250 (MQ=255)
                                                                                      cacTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCggggg                                                                                                                                                         >  1:77920/1‑251 (MQ=255)
                                                                                        cTAATGACTTATCGATGGCCTCTATTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGCCTGCGGTGTTGACATCGGGGGATCACGCATCCCTGGCGGGGGGAATTGTCGGCTTATTTTTTCCTTTGAAAGAGAAGCATGGGCGATCTACAGCGAAACGACTGGCACATGTGttaccccatgggggg                                                                                                                                                       >  2:18921/1‑238 (MQ=255)
                                                                                           aTGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCCACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCg                                                                                                                                                    >  1:282436/1‑251 (MQ=255)
                                                                                             gattttgcGATTGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTGCTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACACGTGGGTGGCGTa                                                                                                                                                  <  2:282183/244‑1 (MQ=255)
                                                                                              ccttttcggtggcccttctTGGCGCCGCGGCTGTAGCAATTGCGGTACTCGCTGTAATGAATCTTTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACTCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGGGTTACTCCCGTGGGTGGCGTAt                                                                                                                                                 <  2:439963/235‑1 (MQ=255)
                                                                                                     gATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGAttt                                                                                                                                          <  1:185467/251‑1 (MQ=255)
                                                                                                        ggCCTCTCTTGGCGTCGCGGCTGTAGCCATTGCGGTACTCGCGGTATTGAATCTGGGTGGTGTACGCCGCACGGGCGTCTATATTCTGGGTGGCGTAGGGTTGTGGACTGCGGTGTTGAGATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGATTCTTTTTTTCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGGAACGACTGGGGCGTGTGTTGCACACGTGGGGGGCGGATCTGCTTTTGc                                                                                                                                       >  2:423143/1‑251 (MQ=255)
                                                                                                                     gTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCAtt                                                                                                                          <  1:287907/251‑1 (MQ=255)
                                                                                                                     gTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAAGCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTATTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTTGGTGGCGTATCTGATTTTGCCGCTGTTTGCAtt                                                                                                                          <  1:15718/251‑1 (MQ=255)
                                                                                                                                   aaTTGCGGTACTCGCGTTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGTGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTGGCTAATGCTGGCg                                                                                                            <  2:98053/251‑1 (MQ=255)
                                                                                                                                        cGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCa                                                                                                       <  1:401710/251‑1 (MQ=255)
                                                                                                                                        cGGTACTCGCGGTATTGAAGCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCCGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCa                                                                                                       <  1:115283/251‑1 (MQ=255)
                                                                                                                                             cTCGCGGTATTGAATCTGTGTGGGGTACGCCGCACGGGCGTCTATATTCTGGGTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGCGAAGCATGGGCGGTCTCCAGCGAAGTGTCTGGTGCATGTGGTGCACCCCGGGGGGGCGTATCTGGTTTTGCCGCTATTTGCCTTTGCTAATCCTGGCGTTTCAATGc                                                                                                   >  2:12370/1‑250 (MQ=255)
                                                                                                                                                gcgGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTCTATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGAAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCAAGGGCGTTCTCCAGCGAAGCGGCTGGAGCATGTGTTACACCCGTGGGTGGCGTCTCTGATTTTGCCGCTGTTTGCGTTTGCTAATGCTGGCGTTTCACTGCAAgg                                                                                               >  2:353803/1‑251 (MQ=255)
                                                                                                                                                   gTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGAGTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGTGGGGGGTAATTGTCGGCTTCTGTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAggtgt                                                                                            <  1:458593/251‑4 (MQ=255)
                                                                                                                                                       tGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACg                                                                                        >  1:205500/1‑251 (MQ=255)
                                                                                                                                                                gtggtgTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGAGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATgg                                                                                <  1:300010/250‑1 (MQ=255)
                                                                                                                                                                   gtgtaCGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTg                                                                            >  1:121523/5‑251 (MQ=255)
                                                                                                                                                                         gccgcTCGGGCGTCTATATTCTGGTTTGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTcc                                                                      <  2:284707/251‑1 (MQ=255)
                                                                                                                                                                         gccgcACGGGCGGCTATAGTCTGGTTGGCGTAGGGGTGTGGACTGCGGTGTTGGAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTcc                                                                      <  1:268226/251‑1 (MQ=255)
                                                                                                                                                                                  gCGTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGcc                                                              >  1:482091/1‑250 (MQ=255)
                                                                                                                                                                                    gTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGCTGGCTTGACCTCCATTCTGCCAtt                                                           >  1:280983/1‑251 (MQ=255)
                                                                                                                                                                                    gTCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCAtt                                                           <  1:126221/251‑1 (MQ=255)
                                                                                                                                                                                     tCTATATTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGATAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTg                                                          >  1:18696/1‑251 (MQ=255)
                                                                                                                                                                                       tataTTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTgg                                                         <  1:356597/250‑1 (MQ=255)
                                                                                                                                                                                          atTCTGGTTGGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGat                                                     >  1:274726/1‑251 (MQ=255)
                                                                                                                                                                                             ctggttgGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGAACTCCATTCTGCCATTGGGGatcat                                                  >  1:484627/4‑251 (MQ=255)
                                                                                                                                                                                               ggttgGCGTAGTGTTGTGGACTGCGGTTGTGACATCGGGGGTTCACGCAACCCTGGTGGGGGTAATTGGCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGGTCTCCAGCGAAGCGACTGGTGCATGTGTTTCACCCGTGGGTGGCGTATCTGATGTGGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGACACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGatcat                                                  <  1:205900/248‑1 (MQ=255)
                                                                                                                                                                                               ggttgGCGTAGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGatcatc                                                 <  2:234406/249‑1 (MQ=255)
                                                                                                                                                                                                             tgtgGACTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGAtt                                  >  1:82616/1‑251 (MQ=255)
                                                                                                                                                                                                                   cTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGGGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCCCTGTTTGCAt                                                                                                                           <  2:349254/156‑1 (MQ=255)
                                                                                                                                                                                                                   cTGCGGTGTTGAAATCGGGGGTTCACGCAACCCTGGCGGGGGGAATTGTCGGCTTCTTTTTTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAACGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCAt                                                                                                                           >  1:349254/1‑156 (MQ=255)
                                                                                                                                                                                                                           ttGAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTggg                    >  1:49480/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAAATCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGa                  <  1:373634/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaTCGGGGGTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGAtt                <  1:241389/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggggTTCACGCAACCCTGGCGGGGGTAATTGTCTGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTTCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCt           <  1:433046/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggggTTCACGCAACCCTGGCGGGGGTAATTGTCGGCTTCTTTATTACTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCt           <  1:115876/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gggTTTACGCAACCCTGGCGGGGGTAATTGTCGGCTGCGTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTg          <  2:26166/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gCAACCCTGGCGGGGGGAATTGGCGGCTTCTTTATTCCTTTGAAAGCGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCACACCGCTGGGGGTTAGTCTGTTctgatg  >  2:362755/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
CGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGACTGCGGTGTTGAAATCGGGGGTTCACGCAACTCTGGCGGGGGTAATTGTCGGCTTCTTTATTCCTTTGAAAGAGAAGCATGGGCGTTCTCCAGCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTG  >  NC_000913/17927‑18414

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