Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 762,013 T→G 5.7% A164A (GCT→GCG sucB → dihydrolipoyltranssuccinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913762,0130TG 5.7% 263.2 / 2.3 77A164A (GCT→GCGsucBdihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base G (2/4);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACG  >  NC_000913/761763‑762281
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                            
aCAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGAAACGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAACAACATCTGGCGAAAGCCCCGGCGAACACGTCTGTTCCGGAAGCGGGTGTGCCGGAGGGGCGAACGGggctggct                                                                                                                                                                                                                                                      >  1:20326/1‑275 (MQ=255)
                                           ttCCACCCCGGCGCAACGCCAGCAGTCGGCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGc                                                                                                                                                                                                           <  2:402429/274‑1 (MQ=255)
                                                cTCCGGCGCAAGCCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGTTCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCgtgt                                                                                                                                                                                                       <  2:79822/274‑1 (MQ=255)
                                                   cGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGc                                                                                                                                                                                                    <  2:116384/274‑1 (MQ=255)
                                                                                gagCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCa                                                                                                                                                                       <  2:384474/274‑1 (MQ=255)
                                                                                             aTGCGTTGAGCCCGGCGAGCCGTCGCCTGCTGGCTGAAACCAATCTCGACGCCAGCGCCATTAAAGGCACCGGGTTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCGGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGAcc                                                                                                                                                         <  2:479178/275‑1 (MQ=255)
                                                                                                cGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGATGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCGCGTAGTGAAAAACGTGTCCCGATGACCCGCCTGCGTACGCGTGTGGCAGAGCGTCTGCCGGAAGCGAAAAACTCCACCGCCATGCTGACCACg                                                                                                                                                      >  2:93878/1‑275 (MQ=255)
                                                                                                     ggCCCGGGGATCCGTGGCCGGCTGGCTGAAGAAAATCTCGGCGCCAGCGCCAATAAAGGCACCGGTGTGGGGGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCaa                                                                                                                                                 <  4:263310/274‑1 (MQ=255)
                                                                                                      gCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGgctgct                                                                                                                                                                                                                                                                              >  7:156168/1‑149 (MQ=255)
                                                                                                          ggCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCCGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGaa                                                                                                                                             >  2:361481/1‑274 (MQ=255)
                                                                                                             gATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCTGGCAGCGGCTGCTCcggcgg                                                                                                                                                                                                                                                                       >  6:365115/1‑149 (MQ=255)
                                                                                                                  gTCGCCTGCTGGCTGAACACATTCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATg                                                                                                                                    <  1:56879/275‑1 (MQ=255)
                                                                                                                              ctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCt                                                                                                                                                                                                                                                      >  5:329948/1‑149 (MQ=255)
                                                                                                                                      aaTCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAg                                                                                                                                                                                                                                              <  6:213595/149‑1 (MQ=255)
                                                                                                                                      aaTCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAg                                                                                                                                                                                                                                              >  8:68354/1‑149 (MQ=255)
                                                                                                                                         ctcGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGa                                                                                                                                                                                                                                           >  5:62476/1‑149 (MQ=255)
                                                                                                                                             aCGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGAGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGc                                                                                                                                                                                                                                                    >  2:109007/1‑136 (MQ=255)
                                                                                                                                              cGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAAc                                                                                                                                                                                                                                      <  8:7221/149‑1 (MQ=255)
                                                                                                                                                ccAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACgt                                                                                                                                                                                                                                    >  6:200251/1‑149 (MQ=255)
                                                                                                                                                ccAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACCTTAAGCCGATTATGGATCTGCGTAACCAGTAc                                                                                                      >  4:336881/1‑275 (MQ=255)
                                                                                                                                                 cAGCGCCATTACAGGCACCGGAGTGGGTGGTCGTATGACTCGTGAAGATGTGGAAAAACATATGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCGCCGGCGGCGCAACCGGCTCTGGCTGCACGt                                                                                                                                                                                                                                                <  4:109006/136‑1 (MQ=255)
                                                                                                                                                  aGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACgtgt                                                                                                                                                                                                                                  <  7:459053/149‑1 (MQ=255)
                                                                                                                                                     ggcaggaagcgaaCCGGGGGGGGTGGTCGGCGAACTCGTGGGCATGTGGAAAAGCAGCGGGCGAAAGCCGCGGCGAAAAAGTCTGCTCCGGCAGCGGCTGATCCGGCGGCGCAACCGGCTCGGGCTGCGCGTAGTGAAAAACGTGTACCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCGATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTg                                                                                                  <  4:93877/262‑1 (MQ=255)
                                                                                                                                                         ttAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATg                                                                                                                                                                                                                           >  7:75211/1‑149 (MQ=255)
                                                                                                                                                          tAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGa                                                                                                                                                                                                                          <  5:275230/149‑1 (MQ=255)
                                                                                                                                                                 aCCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCACGCTGACTCGCATGCGTAAGCGTTTGGCAGAGCGTCTGCTGGCAGCGAACAACTCCACCGCAATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTTTGGCTATCCGTAAGCATTCCGGtacaggcattaaaa                                                                                      >  4:309732/1‑261 (MQ=255)
                                                                                                                                                                   cGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGCCTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAACCGTGTGaaaaa                                                                                    >  2:234994/1‑274 (MQ=255)
                                                                                                                                                                   cGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGc                                                                                                                                                                                                                 >  2:365221/1‑149 (MQ=255)
                                                                                                                                                                     gtgtGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGt                                                                                                                                                                                                               <  7:68354/149‑1 (MQ=255)
                                                                                                                                                                     gtggGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATg                                                                                                                                    <  2:492213/224‑1 (MQ=255)
                                                                                                                                                                       gtgggtggTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTaa                                                                                                                                                                                                             <  4:365216/149‑1 (MQ=255)
                                                                                                                                                                               tcgtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTgg                                                                                                                                                                                                     <  8:37891/149‑1 (MQ=255)
                                                                                                                                                                                 gtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAAc                                                                                                                                                                                                                                      >  6:131290/1‑114 (MQ=255)
                                                                                                                                                                                 gtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAAc                                                                                                                                                                                                                                      <  5:131290/114‑1 (MQ=255)
                                                                                                                                                                                 gtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGAt                                                                                                                                                                                                                            <  8:413783/124‑1 (MQ=255)
                                                                                                                                                                                 gtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGAt                                                                                                                                                                                                                            >  7:413783/1‑124 (MQ=255)
                                                                                                                                                                                 gtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCGGCACGAAGGGAAAAACTTGCCCCGATGACTCGCCTGCGTAAGCGGGTggaa                                                                                                                                                                                                   >  8:84446/1‑147 (MQ=255)
                                                                                                                                                                                  tcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCag                                                                                                                                                                                                  >  5:369002/1‑149 (MQ=255)
                                                                                                                                                                                    ggACGCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCGGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCt                                                                   <  2:309736/273‑1 (MQ=255)
                                                                                                                                                                                        tCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTc                                                                                                                                                                                            >  8:131876/1‑149 (MQ=255)
                                                                                                                                                                                        tCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTc                                                                                                                                                                                            >  5:299905/1‑149 (MQ=255)
                                                                                                                                                                                        tCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTc                                                                                                                                                                                            <  6:299905/149‑1 (MQ=255)
                                                                                                                                                                                           tGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTctgc                                                                                                                                                                                         >  8:54194/1‑149 (MQ=255)
                                                                                                                                                                                            gAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTctgct                                                                                                                                                                                        >  7:50959/1‑149 (MQ=255)
                                                                                                                                                                                            gAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTctgct                                                                                                                                                                                        <  6:62476/149‑1 (MQ=255)
                                                                                                                                                                                                    ggAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGa                                                                                                                                                                                >  7:469778/1‑149 (MQ=255)
                                                                                                                                                                                                    ggAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGa                                                                                                                                                                                >  8:134169/1‑149 (MQ=255)
                                                                                                                                                                                                     gAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGaa                                                                                                                                                                               <  6:329948/149‑1 (MQ=255)
                                                                                                                                                                                                     gAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGGCATGCTGACCACGTTCAACGAGGTCAACATGAAACCGATTCTGGATCTGCGTa                                                                                                              >  4:206655/1‑214 (MQ=255)
                                                                                                                                                                                                     gAAAAACATAGAGAGAAAGCCCCGGGCAAAGAGGCTGCGCCGGCAGGGGCGGCGCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCAGGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGACATGCTGACCACGTTCAACGAAGGCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTa                                                 <  1:86079/275‑1 (MQ=255)
                                                                                                                                                                                                         aaCATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAAc                                                                                                                                                                           <  8:75211/149‑1 (MQ=255)
                                                                                                                                                                                                         aaCATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCGCCGGCGGCGCAACCGGCTATGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCa                                                                                                          <  2:206658/214‑1 (MQ=255)
                                                                                                                                                                                                             tCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCa                                                                                                                                                                       >  8:311344/1‑149 (MQ=255)
                                                                                                                                                                                                             tCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCa                                                                                                                                                                       <  8:156168/149‑1 (MQ=255)
                                                                                                                                                                                                                 gccaaaaCCCCGGCCAAAGAGTTCTCCCCGGGAGGCGGTGCTCCGGCGGGGCAAACGGGTGTTGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGc                                                                                                                                                                   <  7:84446/146‑1 (MQ=255)
                                                                                                                                                                                                                 gCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGc                                                                                                                                                                   >  6:299861/1‑149 (MQ=255)
                                                                                                                                                                                                                    aaaGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCAt                                                                                                                                                                >  8:278959/1‑149 (MQ=255)
                                                                                                                                                                                                                      aGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCg                                                                                                                >  2:44912/1‑195 (MQ=255)
                                                                                                                                                                                                                      aGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCCGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGTAAAACGTGTCCCGATGAATCGCCCGCGTAAGCGTGTGGTAGAGCGTTTGCTGGAAGCGAAAAACTcc                                                                                                                                                                        >  2:332048/1‑139 (MQ=255)
                                                                                                                                                                                                                       gCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGa                               <  2:336885/275‑1 (MQ=255)
                                                                                                                                                                                                                          ccGGCGAAAGAGTCTGCTCTGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGAc                                                                                                                                                          <  5:365115/149‑1 (MQ=255)
                                                                                                                                                                                                                           cGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGAGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCCCGTTCAACGAAGTCAACATGAAGCCGATTATGGATCGGCGTAAGCAGTACGGTGAAGCGTTGGAAAAACGCAAAGGCATCCGTCTGGGCCTTATGTCCTTCCACGGTAAACCgggtgttctgtcc                           >  4:271761/1‑264 (MQ=255)
                                                                                                                                                                                                                             gCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGAt                                                                                                                                                                                                                            >  6:87733/1‑80 (MQ=255)
                                                                                                                                                                                                                             gCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGAt                                                                                                                                                                                                                            <  5:87733/80‑1 (MQ=255)
                                                                                                                                                                                                                               gAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGa                                            <  2:271765/254‑1 (MQ=255)
                                                                                                                                                                                                                                   gagTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACg                   <  1:64153/275‑1 (MQ=255)
                                                                                                                                                                                                                                      tCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGtt                 >  2:422744/1‑274 (MQ=255)
                                                                                                                                                                                                                                       cTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGaa                                                                                                                                             <  6:369002/149‑1 (MQ=255)
                                                                                                                                                                                                                                          cTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTc                                                                                                                                          <  6:396938/149‑1 (MQ=255)
                                                                                                                                                                                                                                             cGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAAc                                                                                                                                       <  7:134169/149‑1 (MQ=255)
                                                                                                                                                                                                                                             cGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAAc                                                                                                                                       <  6:112944/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACAt                                                                                                                                     <  7:278959/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACAt                                                                                                                                     <  5:200251/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGt       >  2:452924/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 gccGGTTGCGCCGCCGGAGCAGCCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCa                                                                                                                                                                       >  8:119232/3‑113 (MQ=255)
                                                                                                                                                                                                                                                 gccGGTTGCGCCGCCGGAGCAGCCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCa                                                                                                                                                                       <  7:119232/111‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggctctgctCTCGGGGCGCAACCGGCCCTGGCTGACCGTAGTGAAAACCGGGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGaa    <  2:85191/271‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGaa    >  2:363236/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    ggCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGAAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCCTTATGTCCTTCTACGTGAAAGCGGTGCTTGAAGCCCTTACACGTTCCCCGGAAGTGaaag  >  2:382077/1‑273 (MQ=255)
                                                                                                                                                                                                                                                     gctgctCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCgtg                                                                                                                                                                                                        <  7:398062/76‑1 (MQ=255)
                                                                                                                                                                                                                                                     gctgctCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCgtg                                                                                                                                                                                                        >  8:398062/1‑76 (MQ=255)
                                                                                                                                                                                                                                                      ctgctCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCg                                                                                                                              >  7:418225/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gctCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGAt                                                                                                                            >  6:241648/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTa                                                                                                                          >  5:94099/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                            
ACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAAAGCGTCCACTCCGGCGCAACGCCAGCAGGCGTCTCTGGAAGAGCAAAACAACGATGCGTTAAGCCCGGCGATCCGTCGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCGAAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAAGCGAAAAACTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACG  >  NC_000913/761763‑762281

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: