Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,125 C→T 5.4% T162T (ACC→ACT ycdY → redox enzyme maturation protein (REMP) chaperone for YcdX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,1250CT 5.4% 737.1 / 21.4 220T162T (ACC→ACTycdYredox enzyme maturation protein (REMP) chaperone for YcdX
Reads supporting (aligned to +/- strand):  ref base C (95/113);  new base T (6/6);  total (101/119)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 220 total aligned reads displayed.

GTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCC  >  NC_000913/1099862‑1100368
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
 tgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGACGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGACGCTTAGAGCGAAGCACTGGAAACCCTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAACGTGGGGGGCCATTCAACCACGCCTTTCTGGCCCACCCTTGCAccccca                                                                                                                                                                                                                                         >  4:241992/1‑271 (MQ=255)
                   cATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCCGTGCGTATCCGTTCCCCTGGTGTGCTGCGTTCCTTTGCAACGTGGAGGCCCCTGCAACCCCGCCTTTCTGGCGCACCATGCCACCGCTAACCCGCGATCCCCtt                                                                                                                                                                                                                          >  4:129352/1‑272 (MQ=255)
                   cATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGt                                                                                                                                                                                                                       >  2:192639/1‑275 (MQ=255)
                     tGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTTGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCCCCATGGCCCCGCAAACCCGCGTAGCCATTAGTg                                                                                                                                                                                                                      >  4:320132/1‑274 (MQ=255)
                         ttggggggggcccggaagcggagaCGCGCCCCTGTCTTTCTGAGGGGTGGATGCCATTAGCGGTTAGGCCAGCCGATCACATCGGCACATTGCTTCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACAGTTCAGGGAGTATCTGTTACCCTAGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAAtg                                                                                                                                                                                                                 <  4:347522/263‑1 (MQ=255)
                                                             tttctGGGGTAGGGATGCAATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCCCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgt                                                                                                                                                                              <  4:532945/274‑1 (MQ=255)
                                                                              ccATTAGCGGATAAGGAGGCCGATACAATGGGGAAATTGGTGCTCGCAGCTTCCTGGCGGGAAGGTCAGGAAACGGAAGATGAGGGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCAt                                                                                                                                                            <  4:7488/275‑1 (MQ=255)
                                                                                      gggtaagCCAGCCGATGACACCGGCCCATTGCTGCTCGCAGCTTCCTAGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCCCTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATACAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCa                                                                                                                                                     <  4:411479/268‑1 (MQ=255)
                                                                                                      cacaTCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgt                                                                                                                                     <  1:2072/274‑1 (MQ=255)
                                                                                                      cacaTAGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgt                                                                                                                                     <  2:484759/274‑1 (MQ=255)
                                                                                                                      tgctCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCt                                                                                                                    <  2:320133/275‑1 (MQ=255)
                                                                                                                       gctCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAACTGTTTGCt                                                                                                                    >  2:398752/1‑274 (MQ=255)
                                                                                                                            cAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCAcc                                                                                                                                                                                                                                             >  5:957064/1‑148 (MQ=255)
                                                                                                                             aGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCg                                                                                                                                                                                                                                            <  6:887567/148‑1 (MQ=255)
                                                                                                                                    tggctggAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAAccc                                                                                                                                                                                                                                     >  5:677267/1‑148 (MQ=255)
                                                                                                                                    tagcagaaaacgcggccatgaatgagatcgagcacacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACCATGGCACCGCTAAccc                                                                                                                                                                                                                                     <  5:145924‑M1/21‑1 (MQ=255)
                                                                                                                                     ggctggAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCcg                                                                                                                                                                                                                                    >  8:562754/1‑148 (MQ=255)
                                                                                                                                           aaGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGc                                                                                                                                                                                                                              >  6:602403/1‑148 (MQ=255)
                                                                                                                                               tCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          <  7:960116/148‑1 (MQ=255)
                                                                                                                                                  gTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCa                                                                                                                                                                                                                                                     >  8:305883/1‑119 (MQ=255)
                                                                                                                                                    cAACGGAAGATGAGAGCGAAGCACTGGAAACAGTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                                     >  7:721171/1‑148 (MQ=255)
                                                                                                                                                    cAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                                     >  5:232968/1‑148 (MQ=255)
                                                                                                                                                       cGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCAGCGATGCCATTAgatcgga                                                                                                                                                                                                                  >  5:376569/1‑142 (MQ=255)
                                                                                                                                                            gATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAATGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgcgc                                                                              <  2:115577/275‑1 (MQ=255)
                                                                                                                                                              atcgagcacacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATg                                                                                                                                                                                                                               <  5:668743‑M1/28‑1 (MQ=255)
                                                                                                                                                               gagaGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTggg                                                                                                                                                                                                             <  5:398412/145‑1 (MQ=255)
                                                                                                                                                                     gAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAAt                                                                                                                                                                                                                  <  2:565454/134‑1 (MQ=255)
                                                                                                                                                                      cacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCATTAg                                                                                                                                                                                                                        >  8:1129154‑M1/94‑128 (MQ=255)
                                                                                                                                                                      cacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCATTAg                                                                                                                                                                                                                        <  7:1129154‑M1/35‑1 (MQ=255)
                                                                                                                                                                      aaGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGga                                                                                                                                                                                                   <  7:562754/148‑1 (MQ=255)
                                                                                                                                                                         cACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAAcc                                                                                                                                                                                                                                      <  6:1051675/110‑1 (MQ=255)
                                                                                                                                                                          ttttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACCATGGCACCGCTAACCCGCGATGc                                                                                                                                                                                                                              <  8:351795‑M1/28‑1 (MQ=255)
                                                                                                                                                                          aCTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGgaagaa                                                                                                                                                                                               <  7:782168/148‑1 (MQ=255)
                                                                                                                                                                           tttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATg                                                                                                                                                                                                                               >  7:676895‑M1/89‑116 (MQ=255)
                                                                                                                                                                            tGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAagaaga                                                                                                                                                                                             <  5:525895/148‑1 (MQ=255)
                                                                                                                                                                              gAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCa                                                                                                                                                                                                                                               <  8:244235/96‑1 (MQ=255)
                                                                                                                                                                               aaaCACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGga                                                                                                                                                                                                   <  7:299041/139‑1 (MQ=255)
                                                                                                                                                                               aaaCACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTc                                                                                                                                                                                          >  5:115877/1‑148 (MQ=255)
                                                                                                                                                                                 acacTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAAtg                                                                                                                                                                                                                 <  7:1024675/123‑1 (MQ=255)
                                                                                                                                                                                 acacTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGAt                                                                                                                                                                                            >  5:792569/1‑144 (MQ=255)
                                                                                                                                                                                   acTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAAt                                                                                                                                                                                  >  4:562413/1‑152 (MQ=255)
                                                                                                                                                                                      acaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCAt                                                                                                                                                                                                                           >  6:785351‑M1/78‑109 (MQ=255)
                                                                                                                                                                                       ttCAGTGAGTATCTGTTAACCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATa                                                                                                                               <  2:108718/199‑1 (MQ=255)
                                                                                                                                                                                          aGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCt                                                                                                                    <  4:112528/207‑1 (MQ=255)
                                                                                                                                                                                           gTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgt                                                                                                                                                                              >  6:1065663/1‑148 (MQ=255)
                                                                                                                                                                                           gTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgt                                                                                                                                                                              >  8:847740/1‑148 (MQ=255)
                                                                                                                                                                                            tatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACCATGGCACCGCTAACCCGCGATg                                                                                                                                                                                                                               <  5:191821‑M1/27‑1 (MQ=255)
                                                                                                                                                                                            tGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          >  6:979456/1‑103 (MQ=255)
                                                                                                                                                                                            tGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgcgc                                                                              <  1:116570/243‑1 (MQ=255)
                                                                                                                                                                                               ccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTATGGCACCGCTAACCCGCGAt                                                                                                                                                                                                                                <  5:1031902‑M1/26‑1 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTgg                                                                                                                                                                                                    >  6:323318/1‑120 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGGTTCTga                                                                                                                                                                                       >  8:884002/1‑133 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  5:1061176/133‑1 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  7:884002/133‑1 (MQ=255)
                                                                                                                                                                                                   cTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      <  6:313345/148‑1 (MQ=255)
                                                                                                                                                                                                   cTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      >  7:552320/1‑148 (MQ=255)
                                                                                                                                                                                                     gTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCcgcg                                                                                                                                                                                                                                  <  6:188973/86‑1 (MQ=255)
                                                                                                                                                                                                     gTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGgaagaa                                                                                                                                                                                               <  6:463730/121‑1 (MQ=255)
                                                                                                                                                                                                      ttACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          >  5:561126/1‑93 (MQ=255)
                                                                                                                                                                                                      ttACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGga                                                                                                                                                                                                   <  5:315557/116‑1 (MQ=255)
                                                                                                                                                                                                      ttACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  5:87729/128‑1 (MQ=255)
                                                                                                                                                                                                       tACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTa                                                                                                                                                                                                                         <  7:1051564/93‑1 (MQ=255)
                                                                                                                                                                                                       tACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      >  6:719965/1‑144 (MQ=255)
                                                                                                                                                                                                       tACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      <  5:719965/144‑1 (MQ=255)
                                                                                                                                                                                                           cTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAAtg                                                                                                                                                                                                                 >  7:233833/1‑97 (MQ=255)
                                                                                                                                                                                                            tggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTc                                                                                                                                                                                          <  6:673163/119‑1 (MQ=255)
                                                                                                                                                                                                            tggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       >  6:1010941/1‑122 (MQ=255)
                                                                                                                                                                                                             ggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGATGAtt                                                                                                                                                                                           >  8:741284/1‑117 (MQ=255)
                                                                                                                                                                                                             ggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATa                                                                                                                                                                                 <  8:1049619/127‑1 (MQ=255)
                                                                                                                                                                                                              gtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTg                                                                                                                                                       <  1:124331/152‑1 (MQ=255)
                                                                                                                                                                                                                        ttCCTTGGCAAAGTGGGGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGAAGAGCTGGAAGACGANTCTGAAGAATAAGTGTGATCTACNNCACCAATANCTGCAACGGATAATttg                                                                                                                                         >  3:30548/1‑156 (MQ=255)
                                                                                                                                                                                                                        ttCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACgg                                                                                                                                                 >  7:900899/1‑148 (MQ=255)
                                                                                                                                                                                                                         cacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          <  6:513743‑M1/33‑1 (MQ=255)
                                                                                                                                                                                                                            ttGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  6:118844/106‑1 (MQ=255)
                                                                                                                                                                                                                              ggCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGAGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACCTACTTTCCTCTATTGCAATCAACACAGGCATTCTCTCCCGTTACTGGGGTTGGGTg            >  4:99721/1‑275 (MQ=255)
                                                                                                                                                                                                                               gCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAtt                                                                                                                                           <  8:291203/147‑1 (MQ=255)
                                                                                                                                                                                                                               gCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAtt                                                                                                                                           >  7:291203/1‑147 (MQ=255)
                                                                                                                                                                                                                                 aaaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgt                                                                                                                                        <  8:552320/148‑1 (MQ=255)
                                                                                                                                                                                                                                  aaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGa                                                                                                                                                >  5:33999/1‑139 (MQ=255)
                                                                                                                                                                                                                                  aaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgtt                                                                                                                                       <  6:850421/148‑1 (MQ=255)
                                                                                                                                                                                                                                  aaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGTTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTgggg                   <  2:99722/264‑1 (MQ=255)
                                                                                                                                                                                                                                        aGGCCCATGCAACCACGCCTTTCTGGCGCACCATTGCACCGCTAACCCCCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTccc  <  2:535913/275‑1 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGc                                                                                                                                                                                                       <  6:55677/73‑1 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGAtt                                                                                                                                                                                           <  6:99022/85‑1 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttg                                                                                                                                      >  5:261260/1‑138 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttg                                                                                                                                      <  6:261260/138‑1 (MQ=255)
                                                                                                                                                                                                                                                gCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgtg                                                                                                                                                                             >  7:336475/1‑96 (MQ=255)
                                                                                                                                                                                                                                                gCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTc                                                                                                                                                                  >  5:479630/1‑107 (MQ=255)
                                                                                                                                                                                                                                                 cAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCTGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtgtg                                                                                                                        <  8:900899/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGa                                                                                                                                                                          >  7:858687/1‑95 (MQ=255)
                                                                                                                                                                                                                                                       cGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTg                                                                                                                                                                                        <  7:616595/78‑1 (MQ=255)
                                                                                                                                                                                                                                                       cGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTg                                                                                                                                                                                        >  8:616595/1‑78 (MQ=255)
                                                                                                                                                                                                                                                       cGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCg                                                                                                                  <  6:232968/148‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATa                                                                                                                                              <  5:883667/118‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgttg                                                                                                                                   <  6:359927/129‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  6:877432/76‑1 (MQ=255)
                                                                                                                                                                                                                                                           tttCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCa                                                                                                                                 <  8:174873/129‑1 (MQ=255)
                                                                                                                                                                                                                                                            ttCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATaaa                                                                                                                             >  8:677509/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                    CACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgt                                                                                                                                        >  7:274918‑M1/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                      CCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgt                                                                                                                                     <  5:986275‑M1/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCC  >  NC_000913/1099862‑1100368

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