Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,447,850 G→C 7.0% V1247L (GTA→CTA)  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,8500GC 7.0% 913.8 / 47.6 260V1247L (GTA→CTA) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (125/114);  new base C (9/9);  total (136/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-04

Only 100 of 266 total aligned reads displayed.

GACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447586‑4448116
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                          
                                         tGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCCGGCACTGGATTTGCTCTATCAGTTCGGGTTTTAGCAATGCGaa                                                                                                                                                                                                                          >  4:465715/1‑274 (MQ=255)
                                                                 tAAAATCGGCGAGACGTGTGGCATAAGCAATTTAGCGCTCGACACCCAAGGAGTAGGCAACTCCTCCCAGGTAGTGGTCAGAGTCTATGTATTACCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAACCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                                                                                                                  <  1:72699/274‑1 (MQ=255)
                                                                              gggtttgtgtgaaGCAATTTAGTGCTCGNNNCCCGGGGGGTAGGCGACTCCTCCAAGGTTGTGTTGAGCGGTTATGTATTGCCAGGTCTGGAATTGAAATACGGTGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCGATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATGTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                                                                                                                     <  3:153082/263‑1 (MQ=255)
                                                                                gTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACa                                                                                                                                                                                   <  2:378567/274‑1 (MQ=255)
                                                                                                                 gggAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  2:278990/1‑206 (MQ=255)
                                                                                                                         gCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTCGCAATGCGAATATTTGTCTACGGCAGTTTCCGCCACAAACCAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                                                                                                                  >  2:132225/1‑266 (MQ=255)
                                                                                                                              tcctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGTAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATAGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGGNTAGACCAGGCACTGGATGTGCTCTATCAGTTCGAGTTTTAGCAATGGGAATATGTGTCTACGGAGGTTTAGGACACAAAANAGGCAACAGTCACTAGATTACCAACNCCAAGTTAATGGGAGCTTTCAGTaaagaa                                                                                                                                    >  3:94941/1‑270 (MQ=255)
                                                                                                                                 tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:692861/148‑1 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  6:323178/1‑144 (MQ=255)
                                                                                                                                    cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGGGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                <  1:117083/273‑1 (MQ=255)
                                                                                                                                    aaGGTAGTGGTCAGCGGCTGTGTATTGCCATGTCTGCAAGTGAAATACGGCGTGGGGATATTTGACTCTATAGCAACACCCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  4:430549/274‑1 (MQ=255)
                                                                                                                                     aGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGa                                                                                                                                                                                                                                                            >  8:914379/1‑148 (MQ=255)
                                                                                                                                        tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                                                                                                                                                <  7:423507/141‑1 (MQ=255)
                                                                                                                                            ggTCAGCGGTTATGTATTGCCAGGTCTGCAAGTGAAATCCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      <  2:497033/275‑1 (MQ=255)
                                                                                                                                             gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  2:981/1‑163 (MQ=255)
                                                                                                                                             gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGttt                                                                                                                                                                                                                                      >  2:465290/1‑162 (MQ=255)
                                                                                                                                                 gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATTGcc                                                                                                                 >  2:138523/1‑275 (MQ=255)
                                                                                                                                                  cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAgac                                                                                                                                                                                                                                                                  >  6:698403/1‑127 (MQ=255)
                                                                                                                                                     cTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                        >  2:308681/1‑248 (MQ=255)
                                                                                                                                                         gTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                                    <  4:308680/248‑1 (MQ=255)
                                                                                                                                                             tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  8:480045/141‑1 (MQ=255)
                                                                                                                                                             tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  8:1123768/1‑119 (MQ=255)
                                                                                                                                                               ccAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAAAGTCACTGGATGACCAGTGCCAAGTTACTGGGCGATTTCAGTACGCATAACTACCAGTTTGATAGCCTGGGCCACTAACCa                                                                                                   >  4:85642/1‑275 (MQ=255)
                                                                                                                                                                  ggTCTGCAAGTGAAATACGGCGCGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCa                                                                                                                                                                                                           >  4:216458/1‑168 (MQ=255)
                                                                                                                                                                    tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATACCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                              >  2:116152/1‑275 (MQ=255)
                                                                                                                                                                        cAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                                                                                                               <  2:449357/206‑1 (MQ=255)
                                                                                                                                                                           gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                                                                                                                                                                                                                                     <  8:200026/101‑1 (MQ=255)
                                                                                                                                                                                 tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:841085/1‑148 (MQ=255)
                                                                                                                                                                                  aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                                                                                                                                                                                                                             <  7:908533/118‑1 (MQ=255)
                                                                                                                                                                                      cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      >  1:27903/1‑233 (MQ=255)
                                                                                                                                                                                       atGGGTATATTTGACTCTATGGCAACACTCACGTTCCGTTATCTCCTGATGCCTAAGCTATATCTTGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTGTACGCCACAAACAAGGACACAGTCACTCGATGACCATTGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGTCCACTATCCAGGCGCAGTTCCGGGgaacggcacg                                                                           <  4:332727/274‑1 (MQ=255)
                                                                                                                                                                                            tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACAAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                                                                                                                                                    >  7:1055405/1‑117 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     <  6:956572/115‑1 (MQ=255)
                                                                                                                                                                                             atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:84665/87‑1 (MQ=255)
                                                                                                                                                                                                tttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGACTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGt                                                                  <  2:67416/275‑1 (MQ=255)
                                                                                                                                                                                                 ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                                                                                                                                                                                                                                   >  7:481736/1‑97 (MQ=255)
                                                                                                                                                                                                 ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                                                                                                                                                              >  7:214155/1‑118 (MQ=255)
                                                                                                                                                                                                 ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                                                                                                                                                >  6:602983/1‑148 (MQ=255)
                                                                                                                                                                                                 ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                                                                                                                                                >  6:48998/1‑148 (MQ=255)
                                                                                                                                                                                                  tGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                >  6:728201/1‑83 (MQ=255)
                                                                                                                                                                                                   gACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaa                                                                                                                                                                                              <  8:736843/148‑1 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 >  5:375631/1‑80 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:176199/1‑207 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAAAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  4:116489/1‑207 (MQ=255)
                                                                                                                                                                                                    aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                                                                                                                                                >  5:386811/1‑81 (MQ=38)
                                                                                                                                                                                                       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)
                                                                                                                                                                                                            gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  7:729935/94‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGGCTAAGCTATATCTGGCAGCCGTGTCTGGTGGAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGATTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               >  7:112268/1‑148 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                                                                                                                                                              <  7:554530/117‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           >  6:576791/1‑136 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               <  6:119894/148‑1 (MQ=255)
                                                                                                                                                                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                                                                                                                   >  4:270805/1‑160 (MQ=255)
                                                                                                                                                                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACt                                                                                                                                                                             >  5:735218/1‑148 (MQ=255)
                                                                                                                                                                                                                      cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGgtgt                                                                                                                                                                                                                                                                           >  7:986353/1‑52 (MQ=255)
                                                                                                                                                                                                                      cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                                                                                                                                         <  6:944788/134‑1 (MQ=255)
                                                                                                                                                                                                                      cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                                                                                                               <  2:270806/160‑1 (MQ=255)
                                                                                                                                                                                                                        ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                                                                                                                                                              >  6:385274/1‑95 (MQ=255)
                                                                                                                                                                                                                              taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTg                                                                                                                                                                                                                                                        >  5:571992/1‑63 (MQ=255)
                                                                                                                                                                                                                               aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtt                                                                                                                                                                                                                                       >  4:569107/1‑79 (MQ=255)
                                                                                                                                                                                                                               aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                                                                                                                                                         >  8:772118/1‑93 (MQ=255)
                                                                                                                                                                                                                                 cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  4:548391/1‑94 (MQ=255)
                                                                                                                                                                                                                                  gCCGGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  6:686188/1‑93 (MQ=255)
                                                                                                                                                                                                                                   ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  7:300599/71‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                                                                                                                    <  1:83116/142‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               <  6:862535/129‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACt                            >  4:134835/1‑274 (MQ=255)
                                                                                                                                                                                                                                       aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGAAACGGAACGGTACACGGTGACGTTTATCGTATTGACAACGCAACGCTGGccagact                            >  2:23501/1‑269 (MQ=255)
                                                                                                                                                                                                                                        tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  5:82370/1‑134 (MQ=255)
                                                                                                                                                                                                                                            tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     <  7:774241/68‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      <  5:711779/82‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                                                                                                                                                  <  8:801526/134‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                                    >  6:272109/1‑148 (MQ=255)
                                                                                                                                                                                                                                               gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  5:1112369/59‑1 (MQ=255)
                                                                                                                                                                                                                                               gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                                                                                                                                                    >  8:277368/1‑82 (MQ=255)
                                                                                                                                                                                                                                                 tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                <  5:130657/84‑1 (MQ=255)
                                                                                                                                                                                                                                                 tataTCTGGAAGCCGTGTCTAGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCa                                                                                                                                                                                      <  5:998951/110‑1 (MQ=255)
                                                                                                                                                                                                                                                  atatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                                                                                                                                          <  8:1061095/121‑1 (MQ=255)
                                                                                                                                                                                                                                                  atatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGttt                                                           <  2:523335/232‑1 (MQ=255)
                                                                                                                                                                                                                                                   tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGATGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGTCGCCTTgccc               >  4:51939/1‑273 (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)
                                                                                                                                                                                                                                                         gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                                                                                                                                                         <  5:31925/67‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAGCCGTGTCGGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           >  7:81599/1‑97 (MQ=255)
                                                                                                                                                                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  8:586787/1‑75 (MQ=255)
                                                                                                                                                                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGGACACGGTGAAGTTTATCGTATTGACAACGCCACCCTGGCCGAACTTGTGGCCTTTGCCACCAGGGGCg     >  4:458523/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               gtgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGca                                                                                                                                                  >  6:797122/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCTACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  7:739957/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGTCGATTTCAGTATCGATAAc                                                                                                                                 >  8:860465/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   >  5:593717/1‑146 (MQ=255)
                                                                                                                                                                                                                                                                  tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                                                                                                                                                       <  8:512124/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTc                                                                                                                                             >  8:59737/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                    tGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                                                                                                                                     <  7:330547/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    >  7:593575/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                        gtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                                                                                                                           <  2:565092/146‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: