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.

GACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447586‑4448116
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                       
gACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTTTAGACCAGGc                                                                                                                                                                                                                                                                  >  2:152812/1‑275 (MQ=255)
                                           cgcAAATTGGCCAGATTGTGGGTAAAATCGGCGGGACTTTTGGCGTAAGCAATTTGGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAatat                                                                                                                                                                                                                       <  2:504396/275‑1 (MQ=255)
                                                      cggttgttgGGTAAAATCGGGCGGACGTTTGGCGTAAACAATTTGGCGCTCGGCACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg                                                                                                                                                                                                             <  4:14187/267‑1 (MQ=255)
                                                                 tAAAATCGGCGAGACGTGTGGCATAAGCAATTTAGCGCTCGACACCCAAGGAGTAGGCAACTCCTCCCAGGTAGTGGTCAGAGTCTATGTATTACCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAACCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                                                                                                                  <  1:72699/274‑1 (MQ=255)
                                                                                    tgCGGAAGCAATTTAGCGTGTGACAACCAGGGAGTAGGCGAGTCCTCCCAGGTAGGGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATGTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                                                                                                              <  2:271341/274‑1 (MQ=255)
                                                                                                tgaGCGCTCGCCACCCAGGGGGTAGGCGAACCCTCCCAGGTAGTGGTCAGCGGCTATGTTTTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTGTACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAGTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                                                                                                                   <  4:222817/272‑1 (MQ=255)
                                                                                                                 gggAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  2:278990/1‑206 (MQ=255)
                                                                                                                    agggggggACTCCTCCCAGGTAGTGGTCGGCGGCTAGGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATAGTTGACTCTATAGCAACACTCACGTTACGTTATAGGCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                                                                                                                              <  4:132224/268‑1 (MQ=255)
                                                                                                                     gTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTc                                                                                                                                                                                                                  <  4:278989/206‑1 (MQ=255)
                                                                                                                         gCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTCGCAATGCGAATATTTGTCTACGGCAGTTTCCGCCACAAACCAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                                                                                                                  >  2:132225/1‑266 (MQ=255)
                                                                                                                            acgcctccAAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTTTAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        <  1:94941/270‑1 (MQ=255)
                                                                                                                             ctcctcCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                                                                                                                                     <  2:418659/275‑1 (MQ=255)
                                                                                                                              tcctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGTAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATAGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGGNTAGACCAGGCACTGGATGTGCTCTATCAGTTCGAGTTTTAGCAATGGGAATATGTGTCTACGGAGGTTTAGGACACAAAANAGGCAACAGTCACTAGATTACCAACNCCAAGTTAATGGGAGCTTTCAGTaaagaa                                                                                                                                    >  3:94941/1‑270 (MQ=255)
                                                                                                                                 tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:702346/148‑1 (MQ=255)
                                                                                                                                 tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:692861/148‑1 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacg                                                                                                                                                                                                                                                               >  7:878056/1‑145 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:323178/144‑1 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGGGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                <  1:117083/273‑1 (MQ=255)
                                                                                                                                        tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                                                                                                                                                <  7:423507/141‑1 (MQ=255)
                                                                                                                                             gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  2:981/1‑163 (MQ=255)
                                                                                                                                             gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGttt                                                                                                                                                                                                                                      >  2:465290/1‑162 (MQ=255)
                                                                                                                                                 gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCa                                                                                                                                                                                                                                 <  4:981/163‑1 (MQ=255)
                                                                                                                                                  cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAgac                                                                                                                                                                                                                                                                  >  6:698403/1‑127 (MQ=255)
                                                                                                                                                     cTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        >  2:308681/1‑248 (MQ=255)
                                                                                                                                                         gTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                                    <  4:308680/248‑1 (MQ=255)
                                                                                                                                                             tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  7:480045/1‑141 (MQ=255)
                                                                                                                                                             tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  8:1123768/1‑119 (MQ=255)
                                                                                                                                                                    tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATACCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                              >  2:116152/1‑275 (MQ=255)
                                                                                                                                                                        cAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                                                                                                               <  2:449357/206‑1 (MQ=255)
                                                                                                                                                                          aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTAGACCAGACa                                                                                                                                                                                                                                                                 >  6:681916/1‑106 (MQ=255)
                                                                                                                                                                           gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                                                                                                                                                                                                                                     >  7:200026/1‑101 (MQ=255)
                                                                                                                                                                           gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                                                                                                                                                                                                                                     <  8:200026/101‑1 (MQ=255)
                                                                                                                                                                                 tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  5:403517/121‑1 (MQ=255)
                                                                                                                                                                                 tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:841085/1‑148 (MQ=255)
                                                                                                                                                                                  aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                                                                                                                                                                                                                             >  8:908533/1‑118 (MQ=255)
                                                                                                                                                                                      cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      >  1:27903/1‑233 (MQ=255)
                                                                                                                                                                                        tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGt                                                                                                                                                                                                         >  7:736843/1‑148 (MQ=255)
                                                                                                                                                                                          ggTTTATTTTACTTTTTAGCAACCCTCACGTTACGTTATCGCCTGATGCCTAAGCTTTATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGNTTTGCTNTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGc                                                                                                                  <  3:27903/233‑1 (MQ=255)
                                                                                                                                                                                            tataTTTGACTCTATAGCAACACTCACGTTACGTTAACGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                                                                                                                                                    <  8:1055405/117‑1 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  7:953931/1‑115 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     <  6:956572/115‑1 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacggc                                                                                                                                                                                                                                                             <  6:1040807/91‑4 (MQ=38)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  6:84665/1‑87 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:84665/87‑1 (MQ=255)
                                                                                                                                                                                               atTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACt                                                                                                                                                                                                                                                               >  5:590754/1‑87 (MQ=255)
                                                                                                                                                                                                tttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGACTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGt                                                                  <  2:67416/275‑1 (MQ=255)
                                                                                                                                                                                                 ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                                                                                                                                                >  6:48998/1‑148 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 >  8:423734/1‑80 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 >  5:375631/1‑80 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                                                                                                                                     >  4:569578/1‑172 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:176199/1‑207 (MQ=255)
                                                                                                                                                                                                      tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  8:218250/148‑1 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGc                                                                                                                                                                                                                                                                  <  6:912234/76‑1 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATc                                                                                                                                                                                                                                                 >  5:606869/1‑93 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                                                                                                                                                           <  5:948350/115‑1 (MQ=255)
                                                                                                                                                                                                       ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          >  5:899195/1‑148 (MQ=255)
                                                                                                                                                                                                        tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                >  3:82877/1‑205 (MQ=255)
                                                                                                                                                                                                        tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  4:176198/207‑1 (MQ=255)
                                                                                                                                                                                                         ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                                                                                                                                                         >  5:207944/1‑115 (MQ=255)
                                                                                                                                                                                                            gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  8:729935/1‑94 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGTTCTATCAGCTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTTCGGGGAACGGAACGGTACCCGGTGAGGTTAATGGTATGGAc                                                >  4:470624/1‑275 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                                                                                                                                                              <  7:554530/117‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           >  6:576791/1‑136 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               <  6:119894/148‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               >  7:1049410/1‑148 (MQ=255)
                                                                                                                                                                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACt                                                                                                                                                                             >  5:735218/1‑148 (MQ=255)
                                                                                                                                                                                                                        ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                                                                                                                                                              >  6:385274/1‑95 (MQ=255)
                                                                                                                                                                                                                              taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTg                                                                                                                                                                                                                                                        >  5:571992/1‑63 (MQ=255)
                                                                                                                                                                                                                               aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                                                                                                                                                         <  7:772118/93‑1 (MQ=255)
                                                                                                                                                                                                                                 tgCCTGATGCCTAAGCTATATTTGGAAGCCTGTGATGGTGTAGACCAGGTACTGGATTTGCTCTATCAGTTCGAGTTTTGGTAAGGGGAAGATTTGTCTACGGCAGTTTAAGCCACAAAAAAGGCTACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCCGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGTGAACGTAACGGTACACGGTGAAGTTTTTCGTATTGACAACGCCACGCTGGcc                                 <  4:440370/274‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  4:548391/1‑94 (MQ=255)
                                                                                                                                                                                                                                   ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAg                                                                                                                                                                                                                                   <  2:569108/79‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGAAACGGAACGGTACACGGTGACGTTTATCGTATTGACAACGCAACGCTGGccagact                            >  2:23501/1‑269 (MQ=255)
                                                                                                                                                                                                                                        tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  5:82370/1‑134 (MQ=255)
                                                                                                                                                                                                                                            tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  8:774241/1‑68 (MQ=255)
                                                                                                                                                                                                                                            tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGc                                                                                                                                                                                                            >  7:1066234/1‑93 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  6:711779/1‑82 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                                    >  6:272109/1‑148 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                                    <  8:594819/148‑1 (MQ=255)
                                                                                                                                                                                                                                              aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  6:328877/1‑128 (MQ=255)
                                                                                                                                                                                                                                              aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                               >  4:523334/1‑232 (MQ=255)
                                                                                                                                                                                                                                               gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  5:1112369/59‑1 (MQ=255)
                                                                                                                                                                                                                                               gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                                                                                                                                                    <  7:277368/82‑1 (MQ=255)
                                                                                                                                                                                                                                                  atatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGttt                                                           <  2:523335/232‑1 (MQ=255)
                                                                                                                                                                                                                                                    atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                                                                                                              >  8:422475/1‑115 (MQ=255)
                                                                                                                                                                                                                                                      cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  7:120906/100‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    >  8:321116/1‑105 (MQ=255)
                                                                                                                                                                                                                                                        ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  2:540561/159‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  8:81599/97‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                                                                                                                  <  7:778222/138‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        >  8:76422/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:757407/1‑75 (MQ=255)
                                                                                                                                                                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  8:586787/1‑75 (MQ=255)
                                                                                                                                                                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGGACACGGTGAAGTTTATCGTATTGACAACGCCACCCTGGCCGAACTTGTGGCCTTTGCCACCAGGGGCg     >  4:458523/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  6:593717/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCAACTATCCAGGCGCCGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACCCCACGCATGCCGAAATTGATGCCTTTCGCCCCAGGGGCGGTg  >  4:344483/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                                                                                                                                                       <  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)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                       
GACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  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: