Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,770,539 C→T V281M (GTG→ATG)  yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,5390CT94.5% 195.2 / ‑0.9 74V281M (GTG→ATG) yibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (38/31);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

CCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAC‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCT  >  NC_000913/3770297‑3770784
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
ccGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTcc                                                                                                                                                                                                                                                 <  1:109829/250‑1 (MQ=255)
  gCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGGTAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGc                                                                                                                                                                                                                                               >  2:188933/1‑250 (MQ=255)
       aCCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCCGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCatcatc                                                                                                                                                                                                                                         >  2:303269/1‑251 (MQ=255)
       aCCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCatcatc                                                                                                                                                                                                                                         >  1:94048/1‑251 (MQ=255)
          gAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGGTCCCAGCACA‑CCGGCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGACTGCTAAGAACCGCCTGGCACGACAGGGAAAATACTAGTCAGTTTGCCGGGAAACACCTGCCCAGGGAGCGCGTTAAACAT‑CACTTCAGCCTCATCGcc                                                                                                                                                                                                                                      >  2:551798/1‑251 (MQ=255)
                    tGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTc                                                                                                                                                                                                                               >  1:177039/1‑248 (MQ=255)
                      aaaTGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCAAGACAGGTAAAATACTAGTCAGTTTACCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGACAGGTTTCAGACg                                                                                                                                                                                                                          >  1:42145/1‑251 (MQ=255)
                      aaaTGGTCGGAGTAAACCGCCACCTGGGAGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGATAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCACAGCTAGCGCGTTAAA‑ACTCACTTCCGCATCATCGCCCGGGTTTAGACg                                                                                                                                                                                                                          >  2:634558/1‑251 (MQ=255)
                       aaTGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACCGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGcatcat                                                                                                                                                                                                                                          >  1:588042/1‑234 (MQ=255)
                       aaTGGTCGGAGAAAACGTCCACCGGGGCGTAGATGCCGACGGTGAAGGCATCGATATCATCGTTAGGGTCCAGTTAAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTTCCCAGGTAGCGCGTTAAACAT‑CACTGCCGcatcat                                                                                                                                                                                                                                          <  2:588042/234‑1 (MQ=255)
                        aTGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCTCGACCGTTAATGATTGCAATACCCCCAGCGCCTGCTAAGAACCGCCTGGCACGGCAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCCTCCCAGGAAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAAACgg                                                                                                                                                                                                                         >  1:461855/1‑249 (MQ=255)
                             cGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGc                                                                                                                                                                                                                   >  1:51317/1‑251 (MQ=255)
                              ggAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCg                                                                                                                                                                                                                  <  1:188933/251‑1 (MQ=255)
                                   aaatcGCCCCCTGGGCGAAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCCCA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGttt                                                                                                                                                                                                             <  2:303724/247‑1 (MQ=255)
                                      ccGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGGCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCACGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGc                                                                                                                                                                                                          >  2:511695/1‑251 (MQ=255)
                                                ggCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTggg                                                                                                                                                                                                >  1:435981/1‑251 (MQ=255)
                                                 gCGTGGATGCCGCCGGCAAAGGCATCTATATGATTGATAGGGTCTAGTGCAATGTTTCCCAGCACC‑CCGTCTGTGCCGGTCCCGACCGTTAATGATGTCAATACCCCCTGCGCCTGATAAAAACCGCATGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCACGTTATACAT‑CTCTGCCGCATCATCGCCAGGTTTCAGACGTAACCGCGATTTTTGctcaaattttga                                                                                                                                                                                               <  2:3173/251‑12 (MQ=255)
                                                  ggTAGATGCCGTCGGTTAAGGTATCTATATCCTCGTTAGGGTCCAGTTCAATGGTTCCTAGCACC‑CCGGCCGTGACGGGCACTACCGTTAATGAATTTACTACCCCCTTCCCCTGAGAAGAACCGCTTGGCACGACAGGTAAAATACTAGTCATTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTGCCGCTTCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCg                                                                                                                                                                                              <  2:569181/250‑1 (MQ=255)
                                                                      gCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTc                                                                                                                                                                          >  1:430965/1‑251 (MQ=255)
                                                                          cGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCgggg                                                                                                                                                                      <  1:634558/251‑1 (MQ=255)
                                                                              atcatcGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGttt                                                                                                                                                                                                                                <  1:463048/189‑1 (MQ=255)
                                                                              atcatcGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGAACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGAAGCGCGTTAAACAT‑CACTTCCGAATCATCGCCAGGttt                                                                                                                                                                                                                                >  2:463048/1‑189 (MQ=255)
                                                                                catcGTTAGGGTCCAGTTCAATGGTTCCCAGCACC‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAg                                                                                                                                                                <  1:551798/251‑1 (MQ=255)
                                                                                catcGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAg                                                                                                                                                                <  1:52285/251‑1 (MQ=255)
                                                                                catcGTTAGGGCACAGTTCAATGGGTCCCAG‑ACATCCGTCCGTGTCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGACTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAg                                                                                                                                                                <  2:2601/251‑1 (MQ=255)
                                                                                 atcGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGa                                                                                                                                                                  >  1:510399/1‑248 (MQ=255)
                                                                                      tAGGGTCCAGTTAAATGGTTCCCAGCAAA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACACCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGGAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGAGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATAGGGCGACAATTTGCCGTTTTTGCTCGGGGGTGTAGACCATc                                                                                                                                                          >  2:82309/1‑251 (MQ=255)
                                                                                       agattCCAGTTCAAGGTTTCCTCGCACA‑CCGGCCGTTCCGGTCTCGACCGTTAATGATTGCAATTCCCCCTGCGCCTTAAAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CTTTGCCGCAACATCGCCAGGTGTCAGACGGAACAGCGAGTTTTGCCTAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATc                                                                                                                                                          <  2:558793/246‑1 (MQ=255)
                                                                                       aGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCTGTGAAACACCTGCCCAGGTAGCGCGTTAAACAC‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATc                                                                                                                                                          <  1:316310/250‑1 (MQ=255)
                                                                                          gTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCg                                                                                                                                                      >  1:186319/1‑251 (MQ=255)
                                                                                          gTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTAGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAC‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAAt                                                                                                                                                                                                    >  1:323745/1‑205 (MQ=255)
                                                                                             cAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGAc                                                                                                                                                   >  1:345887/1‑251 (MQ=255)
                                                                                                 tCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGc                                                                                                                                                 >  1:443855/1‑249 (MQ=255)
                                                                                                       gTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCa                                                                                                                                          >  1:92473/1‑250 (MQ=255)
                                                                                                        ttCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCCGCGGCa                                                                                                                                          >  1:630981/1‑249 (MQ=255)
                                                                                                          cccAGCACC‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATCCCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAgg                                                                                                                                       <  1:82309/250‑1 (MQ=255)
                                                                                                          cccAGAACA‑CCGACCGAGCCGGGCACGACCGTTAAAGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAAAAGCGAGTTTTGCCGCAATTGGGCGACAATTTGCAGTTTTTGCTCGAGGATGAAGACCATCACCGGACGAAGCggaagga                                                                                                                                       >  2:601744/1‑245 (MQ=255)
                                                                                                           ccAGCACA‑CCGTCCGTGCCGGGCACGCCCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAACATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAACCAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGAAACAGCGCGTTTTGCCGAAATTGGGCGCCAATTTGCCGGTTTTGCTCGGGGATGAAGACCACCACCGGCCGCAGCGGCAAGGc                                                                                                                                      >  2:514745/1‑250 (MQ=255)
                                                                                                            aaacaca‑CCGTCCGTGCCGGGCACGACCGTTACTGATTGCCATACACCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CCCTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGACGAAATTGGACGACAATTTGCCGTTTTTGATCGGGGCTGAAGACCATCACCGGACGCCGCGGAAAGGc                                                                                                                                      >  2:543040/4‑249 (MQ=255)
                                                                                                                cacccctccctgtccggcctcGACCGTTAATGTTTGCATTACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGGCAGTTTTCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGTGGATTAAGACCATCACCGGGCTCAGCGGTGAGGCAgttg                                                                                                                                 <  2:630981/231‑4 (MQ=255)
                                                                                                                     catccggGCCGGTCCCGACCGTTAATGATTGCAATACCCCCTGCGCCGGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCTGCAAGGCAGCTGCGTAt                                                                                                                            <  2:430965/249‑1 (MQ=255)
                                                                                                                            gCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTg                                                                                                                     <  1:514745/251‑1 (MQ=255)
                                                                                                                                        cGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGt                                                                                                          >  1:108106/1‑250 (MQ=255)
                                                                                                                                         gTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGAc                                                                                                                                                   <  2:509637/208‑1 (MQ=255)
                                                                                                                                         gTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGAc                                                                                                                                                   >  1:509637/1‑208 (MQ=255)
                                                                                                                                               gATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGt                                                                                                  >  1:671264/1‑251 (MQ=255)
                                                                                                                                                aTTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTAGGGACCTGGGCGGATCAGTACCTGAGt                                                                                                  >  2:463761/1‑250 (MQ=255)
                                                                                                                                                  ggCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAAGACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAAc                                                                                               <  2:292829/250‑1 (MQ=255)
                                                                                                                                                      atccccccGGCGCCTTATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑TACTTCCGCATCATCGCCAGGTTTCAGACGGAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGGTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCGGAGTAACggag                                                                                           <  2:51317/248‑4 (MQ=255)
                                                                                                                                                        aCCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGc                                                                                          <  1:219944/250‑1 (MQ=255)
                                                                                                                                                          ccccAGCGCCTGAAAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCATGAAACCCCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCAt                                                                                       <  2:177039/251‑1 (MQ=255)
                                                                                                                                                                  ccTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgc                                                                               <  1:531001/251‑1 (MQ=255)
                                                                                                                                                                       aaaGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGcgcgcg                                                                          <  2:135550/250‑1 (MQ=255)
                                                                                                                                                                        aaGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCCGGTAGCGCGTTCAACAT‑CCCTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGTATGAAGACCATCACCGGACACAGAGGCAAGGCAGCTGCGTATGAACCTGGGCGGGTCAGTACCTGAGTAACGTAGCCATTGGTCGgcgaggg                                                                          >  2:551712/1‑246 (MQ=255)
                                                                                                                                                                        aaGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTAAAACAT‑CACTTCAGCATCATCGCCAGGCTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTACTCGGGGATGAAGACCATCCCCGGACTCAGCGGCTAGGCAGCTGCGTATGCACCTGGGCGGACCAGTCCCTAAGTAACGTCGCCATTGCACAGcgcgcgc                                                                         >  2:357681/1‑250 (MQ=255)
                                                                                                                                                                           aaCCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACCCCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATg                                                                      >  1:45159/1‑251 (MQ=255)
                                                                                                                                                                           aaCCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGc                                                                                                                  <  2:646605/207‑1 (MQ=255)
                                                                                                                                                                           aaCCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGc                                                                                                                  >  1:646605/1‑207 (MQ=255)
                                                                                                                                                                                ccTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑Gt                                                                <  1:163155/251‑1 (MQ=255)
                                                                                                                                                                                 cTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTc                                                               >  1:341751/1‑251 (MQ=255)
                                                                                                                                                                                   gtcactaccggtaaagaattAGTCAGTTGGCCGGGAAACACCGGCCCAGGTAGCGCGATAAACAT‑CACTTCCGCATCATCGCCAGGTTGCAGACGTAACAGCGAGTTGTGCCGAAATTGGGCGACAATTTGCCGGTTTTGCTCGGGGATGAAGACCATCACCGTACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTg                                                             <  2:273903/232‑1 (MQ=255)
                                                                                                                                                                                        gACAGGAAAAATACTAGGCAGTTTGACGTGCAACACCTGCCCAGGTATCGCGTTAAACAT‑CACTTACGCATCATCGCCAGGTTTCAGACGAAACAGCAAGTTTTGCCGAAATTGGGCGACAATTTACCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGAAGGTACATGGGCGGAACATGACCTGAGTAACGTCGCCATTGCTAGACGAGCGCATGACC‑ATCTGCTCa                                                         >  2:125110/1‑250 (MQ=255)
                                                                                                                                                                                            ggTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCTCGTTAAACAG‑CACGTCCGCATCATCGCCAGGTTGCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGTGATGAAGACCATCACCGGACGCCGCGGCAAGGCAGCTGCGTATGTACCTGGGAGGATCAGTACCTGATGAACTTATCCATTGCTCGGCGCGCGAATGACA‑GTCTGCGCAAgggt                                                    <  2:101750/251‑3 (MQ=255)
                                                                                                                                                                                              tAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTtat                                                  >  1:537775/1‑251 (MQ=255)
                                                                                                                                                                                                           aGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGt                                     >  1:387215/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccGTGAAACACCTGCCCTGGTAGCGCGTTAAACAT‑CACTTCCGCATCAGCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGTCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTg                               <  2:186319/251‑1 (MQ=255)
                                                                                                                                                                                                                       aacccatGCCCAGGTAGCGCGTGAAACAT‑CACTTCCGCATCATCGCCATGTTTTCGACGTAACAGCTAGTTTTTCCGACAATGGGCGACAATTTTCCGTTTTTGCTCGTGGATGAAGACCATCACCGGTCGCAGCGGCAAGGCAGCTGCGTATGTTCCTTGGCGTATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATAGTTTGCTTCAGTAAGTTGCGctct                         <  2:540176/245‑1 (MQ=255)
                                                                                                                                                                                                                          cccTGCCCAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTATGACGTAACATCGAGTTGTGCCGAAAGTGGGCGACAAGTTGCCGTTTTTGCTCGGGGATGAAGACCAGCACCGGACGCAGCGGCAAGGCAGCGGCGTAGGTACCTGGGCGGATCAGTGCCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATAGTTTGCTTCAGTAAGTTGCGCTCTTaa                      <  2:94048/250‑1 (MQ=255)
                                                                                                                                                                                                                               cccAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTc                  >  1:271162/1‑250 (MQ=255)
                                                                                                                                                                                                                               cccAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTc                  >  1:497730/1‑250 (MQ=255)
                                                                                                                                                                                                                                ccAGGTAGCGCGTTAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTcac                <  1:570485/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cgTTAAACAC‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGATGTTTGCCGAAATTGGGCGACAAATTGCCGTTGTTGCTCGGGGCTGAAGACCATCTCCGGACGGAGCGGGAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTGCCTGAGTAACGTTGCCATTGCGCGGCGCGCGAATGACA‑GGCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGTGCTCTTAAGCTCGCCATCTg          <  2:306940/247‑1 (MQ=255)
                                                                                                                                                                                                                                            ttAAACAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCCAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAT‑ACAAGTCTGCTCACGGTTATATTTTGCTTCAGTTAGTTGCGCTCTTAAGCTCACAATCTGCGAttt    >  2:664539/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aaaCAT‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATGTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCt  <  1:543040/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
CCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACA‑CCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAC‑CACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACA‑GTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCT  >  NC_000913/3770297‑3770784

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