Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,447,853 G→C 6.5% D1248H (GAC→CAC)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,447,8530GC 6.5% 939.5 / 44.0 266D1248H (GAC→CAC) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (130/118);  new base C (8/9);  total (139/127)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-03

Only 100 of 268 total aligned reads displayed.

GACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447586‑4448116
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                       
gACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTTTAGACCAGGc                                                                                                                                                                                                                                                                  >  2:152812/1‑275 (MQ=255)
                                           cgcAAATTGGCCAGATTGTGGGTAAAATCGGCGGGACTTTTGGCGTAAGCAATTTGGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAatat                                                                                                                                                                                                                       <  2:504396/275‑1 (MQ=255)
                                                      cggttgttgGGTAAAATCGGGCGGACGTTTGGCGTAAACAATTTGGCGCTCGGCACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg                                                                                                                                                                                                             <  4:14187/267‑1 (MQ=255)
                                                                 tAAAATCGGCGAGACGTGTGGCATAAGCAATTTAGCGCTCGACACCCAAGGAGTAGGCAACTCCTCCCAGGTAGTGGTCAGAGTCTATGTATTACCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAACCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                                                                                                                  <  1:72699/274‑1 (MQ=255)
                                                                                    tgCGGAAGCAATTTAGCGTGTGACAACCAGGGAGTAGGCGAGTCCTCCCAGGTAGGGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATGTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                                                                                                              <  2:271341/274‑1 (MQ=255)
                                                                                                tgaGCGCTCGCCACCCAGGGGGTAGGCGAACCCTCCCAGGTAGTGGTCAGCGGCTATGTTTTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTGTACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAGTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                                                                                                                   <  4:222817/272‑1 (MQ=255)
                                                                                                                 gggAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  2:278990/1‑206 (MQ=255)
                                                                                                                    agggggggACTCCTCCCAGGTAGTGGTCGGCGGCTAGGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATAGTTGACTCTATAGCAACACTCACGTTACGTTATAGGCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                                                                                                                              <  4:132224/268‑1 (MQ=255)
                                                                                                                     gTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTc                                                                                                                                                                                                                  <  4:278989/206‑1 (MQ=255)
                                                                                                                         gCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTCGCAATGCGAATATTTGTCTACGGCAGTTTCCGCCACAAACCAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                                                                                                                  >  2:132225/1‑266 (MQ=255)
                                                                                                                            acgcctccAAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTTTAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        <  1:94941/270‑1 (MQ=255)
                                                                                                                             ctcctcCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                                                                                                                                     <  2:418659/275‑1 (MQ=255)
                                                                                                                              tcctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGTAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATAGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGGNTAGACCAGGCACTGGATGTGCTCTATCAGTTCGAGTTTTAGCAATGGGAATATGTGTCTACGGAGGTTTAGGACACAAAANAGGCAACAGTCACTAGATTACCAACNCCAAGTTAATGGGAGCTTTCAGTaaagaa                                                                                                                                    >  3:94941/1‑270 (MQ=255)
                                                                                                                                 tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:702346/148‑1 (MQ=255)
                                                                                                                                 tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:692861/148‑1 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGacacg                                                                                                                                                                                                                                                               >  7:878056/1‑145 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:323178/144‑1 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGGGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                <  1:117083/273‑1 (MQ=255)
                                                                                                                                        tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGacac                                                                                                                                                                                                                                                                <  7:423507/141‑1 (MQ=255)
                                                                                                                                             gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  2:981/1‑163 (MQ=255)
                                                                                                                                             gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGttt                                                                                                                                                                                                                                      >  2:465290/1‑162 (MQ=255)
                                                                                                                                                 gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCa                                                                                                                                                                                                                                 <  4:981/163‑1 (MQ=255)
                                                                                                                                                  cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAgac                                                                                                                                                                                                                                                                  >  6:698403/1‑127 (MQ=255)
                                                                                                                                                     cTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        >  2:308681/1‑248 (MQ=255)
                                                                                                                                                         gTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                                    <  4:308680/248‑1 (MQ=255)
                                                                                                                                                             tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  7:480045/1‑141 (MQ=255)
                                                                                                                                                             tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  8:1123768/1‑119 (MQ=255)
                                                                                                                                                                    tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATACCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                              >  2:116152/1‑275 (MQ=255)
                                                                                                                                                                        cAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                                                                                                               <  2:449357/206‑1 (MQ=255)
                                                                                                                                                                          aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTAGACCAGACa                                                                                                                                                                                                                                                                 >  6:681916/1‑106 (MQ=255)
                                                                                                                                                                           gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCa                                                                                                                                                                                                                                                                     >  7:200026/1‑101 (MQ=255)
                                                                                                                                                                           gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCa                                                                                                                                                                                                                                                                     <  8:200026/101‑1 (MQ=255)
                                                                                                                                                                                 tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  5:403517/121‑1 (MQ=255)
                                                                                                                                                                                 tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:841085/1‑148 (MQ=255)
                                                                                                                                                                                  aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                                                                                                                                                                                                                             >  8:908533/1‑118 (MQ=255)
                                                                                                                                                                                      cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      >  1:27903/1‑233 (MQ=255)
                                                                                                                                                                                        tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGt                                                                                                                                                                                                         >  7:736843/1‑148 (MQ=255)
                                                                                                                                                                                          ggTTTATTTTACTTTTTAGCAACCCTCACGTTACGTTATCGCCTGATGCCTAAGCTTTATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGNTTTGCTNTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGc                                                                                                                  <  3:27903/233‑1 (MQ=255)
                                                                                                                                                                                            tataTTTGACTCTATAGCAACACTCACGTTACGTTAACGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                                                                                                                                                    <  8:1055405/117‑1 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  7:953931/1‑115 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     <  6:956572/115‑1 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGacacggc                                                                                                                                                                                                                                                             <  6:1040807/91‑4 (MQ=38)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  6:84665/1‑87 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:84665/87‑1 (MQ=255)
                                                                                                                                                                                               atTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACt                                                                                                                                                                                                                                                               >  5:590754/1‑87 (MQ=255)
                                                                                                                                                                                                tttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGACTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGt                                                                  <  2:67416/275‑1 (MQ=255)
                                                                                                                                                                                                 ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                                                                                                                                                >  6:48998/1‑148 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 >  8:423734/1‑80 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 >  5:375631/1‑80 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                                                                                                                                     >  4:569578/1‑172 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:176199/1‑207 (MQ=255)
                                                                                                                                                                                                      tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  8:218250/148‑1 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGc                                                                                                                                                                                                                                                                  <  6:912234/76‑1 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATc                                                                                                                                                                                                                                                 >  5:606869/1‑93 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                                                                                                                                                           <  5:948350/115‑1 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          >  5:899195/1‑148 (MQ=255)
                                                                                                                                                                                                        tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                >  3:82877/1‑205 (MQ=255)
                                                                                                                                                                                                        tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  4:176198/207‑1 (MQ=255)
                                                                                                                                                                                                         ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                                                                                                                                                         >  5:207944/1‑115 (MQ=255)
                                                                                                                                                                                                            gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  8:729935/1‑94 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGTTCTATCAGCTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTTCGGGGAACGGAACGGTACCCGGTGAGGTTAATGGTATGGAc                                                >  4:470624/1‑275 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                                                                                                                                                              <  7:554530/117‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           >  6:576791/1‑136 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               <  6:119894/148‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               >  7:1049410/1‑148 (MQ=255)
                                                                                                                                                                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACt                                                                                                                                                                             >  5:735218/1‑148 (MQ=255)
                                                                                                                                                                                                                        ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                                                                                                                                                              >  6:385274/1‑95 (MQ=255)
                                                                                                                                                                                                                              taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTg                                                                                                                                                                                                                                                        >  5:571992/1‑63 (MQ=255)
                                                                                                                                                                                                                               aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                                                                                                                                                         <  7:772118/93‑1 (MQ=255)
                                                                                                                                                                                                                                 tgCCTGATGCCTAAGCTATATTTGGAAGCC‑TGTGATGGTGTAGACCAGGTACTGGATTTGCTCTATCAGTTCGAGTTTTGGTAAGGGGAAGATTTGTCTACGGCAGTTTAAGCCACAAAAAAGGCTACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCCGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGTGAACGTAACGGTACACGGTGAAGTTTTTCGTATTGACAACGCCACGCTGGcc                                 <  4:440370/274‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  4:548391/1‑94 (MQ=255)
                                                                                                                                                                                                                                   ccTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAg                                                                                                                                                                                                                                   <  2:569108/79‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGAAACGGAACGGTACACGGTGACGTTTATCGTATTGACAACGCAACGCTGGccagact                            >  2:23501/1‑269 (MQ=255)
                                                                                                                                                                                                                                        tGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  5:82370/1‑134 (MQ=255)
                                                                                                                                                                                                                                            tAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  8:774241/1‑68 (MQ=255)
                                                                                                                                                                                                                                            tAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGc                                                                                                                                                                                                            >  7:1066234/1‑93 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  6:711779/1‑82 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                                    >  6:272109/1‑148 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                                    <  8:594819/148‑1 (MQ=255)
                                                                                                                                                                                                                                              aGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  6:328877/1‑128 (MQ=255)
                                                                                                                                                                                                                                              aGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                               >  4:523334/1‑232 (MQ=255)
                                                                                                                                                                                                                                               gCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  5:1112369/59‑1 (MQ=255)
                                                                                                                                                                                                                                               gCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                                                                                                                                                    <  7:277368/82‑1 (MQ=255)
                                                                                                                                                                                                                                                  atatCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGttt                                                           <  2:523335/232‑1 (MQ=255)
                                                                                                                                                                                                                                                    atCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                                                                                                              >  8:422475/1‑115 (MQ=255)
                                                                                                                                                                                                                                                      cTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  7:120906/100‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    >  8:321116/1‑105 (MQ=255)
                                                                                                                                                                                                                                                        ggAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  2:540561/159‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  8:81599/97‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                                                                                                                  <  7:778222/138‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        >  8:76422/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          aaGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:757407/1‑75 (MQ=255)
                                                                                                                                                                                                                                                          aaGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  8:586787/1‑75 (MQ=255)
                                                                                                                                                                                                                                                             ccGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGGACACGGTGAAGTTTATCGTATTGACAACGCCACCCTGGCCGAACTTGTGGCCTTTGCCACCAGGGGCg     >  4:458523/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                tgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  6:593717/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                tgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCAACTATCCAGGCGCCGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACCCCACGCATGCCGAAATTGATGCCTTTCGCCCCAGGGGCGGTg  >  4:344483/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  t‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                                                                                                                                                       <  8:512124/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTc                                                                                                                                             <  7:59737/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                                                                                                                                     <  7:330547/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    >  7:593575/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                       
GACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447586‑4448116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑