Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,447,844 T→C 7.2% S1245P (TCT→CCT)  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,8440TC 7.2% 835.5 / 37.6 252S1245P (TCT→CCT) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (123/110);  new base C (9/9);  total (132/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-02

Only 100 of 264 total aligned reads displayed.

GAGAGCGATCAGAGCGACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447571‑4448116
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
gagaGCGATCAGAGCGACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCgtgt                                                                                                                                                                                                                                                                                   >  2:226970/1‑274 (MQ=255)
                        gCAATGACCTCGATGCTGACTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTCAAATCGGTGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGCAATACGGCGTGGGTATATTTGACTCTTTAGCAACACTCACGTTACGTTAACGCCTTATGCCTAAGCTATTACGGGAAGCCGTTTCGGGTGTAAACCAGGCCCTGAAttg                                                                                                                                                                                                                                                         >  4:75509/1‑274 (MQ=255)
                                                        tGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCACAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGCATATTTGACTATATAGAAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGACGTGTCAGATGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGAAGTGCGAat                                                                                                                                                                                                                         >  2:536975/1‑275 (MQ=255)
                                                                                tAAAATCGGCGAGACGTGTGGCATAAGCAATTTAGCGCTCGACACCCAAGGAGTAGGCAACTCCTCCCAGGTAGTGGTCAGAGTCTATGTATTACCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAACCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                                                                                                                  <  1:72699/274‑1 (MQ=255)
                                                                                               gTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACa                                                                                                                                                                                   <  2:378567/274‑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)
                                                                                                                                         cGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAc                                                                                                                                                                                                                                                                       >  6:692861/1‑148 (MQ=255)
                                                                                                                                         cGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAc                                                                                                                                                                                                                                                                       >  6:702346/1‑148 (MQ=255)
                                                                                                                                             tcctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGTAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATAGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGGNTAGACCAGGCACTGGATGTGCTCTATCAGTTCGAGTTTTAGCAATGGGAATATGTGTCTACGGAGGTTTAGGACACAAAANAGGCAACAGTCACTAGATTACCAACNCCAAGTTAATGGGAGCTTTCAGTaaagaa                                                                                                                                    >  3:94941/1‑270 (MQ=255)
                                                                                                                                               ctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  4:226969/274‑1 (MQ=255)
                                                                                                                                                 cccAGGTAGTGGTCAGCGGCTACGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAAc                                                                                                                                 >  2:298251/1‑274 (MQ=255)
                                                                                                                                                   ctggTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTg                                                                                                                                                                                                                                                                               <  7:572680/128‑1 (MQ=255)
                                                                                                                                                   cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:323178/144‑1 (MQ=255)
                                                                                                                                                   cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGGGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                <  1:117083/273‑1 (MQ=255)
                                                                                                                                                    aGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGa                                                                                                                                                                                                                                                            >  8:914379/1‑148 (MQ=255)
                                                                                                                                                       tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                                                                                                                                                >  8:423507/1‑141 (MQ=255)
                                                                                                                                                          tGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:430550/1‑264 (MQ=255)
                                                                                                                                                           ggTCAGCGGTTATGTATTGCCAGGTCTGCAAGTGAAATCCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      <  2:497033/275‑1 (MQ=255)
                                                                                                                                                            gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  2:981/1‑163 (MQ=255)
                                                                                                                                                               aGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAg                                                                                                                   >  4:202437/1‑274 (MQ=255)
                                                                                                                                                                gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCa                                                                                                                                                                                                                                 <  4:981/163‑1 (MQ=255)
                                                                                                                                                                gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATTGcc                                                                                                                 >  2:138523/1‑275 (MQ=255)
                                                                                                                                                                 cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAgac                                                                                                                                                                                                                                                                  <  5:698403/129‑3 (MQ=255)
                                                                                                                                                                        gTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                                    <  4:308680/248‑1 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  7:480045/1‑141 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  8:480045/141‑1 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  7:1123768/119‑1 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  8:1123768/1‑119 (MQ=255)
                                                                                                                                                                                   tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATACCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                              >  2:116152/1‑275 (MQ=255)
                                                                                                                                                                                     tGCAAGTGAAATACGGCGCGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGttt                                                                                                                                                                                                       <  2:216459/168‑1 (MQ=255)
                                                                                                                                                                                          gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                                                                                                                                                                                                                                     <  8:200026/101‑1 (MQ=255)
                                                                                                                                                                                          gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  7:218250/1‑148 (MQ=255)
                                                                                                                                                                                                 aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                                                                                                                                                                                                                             <  7:908533/118‑1 (MQ=255)
                                                                                                                                                                                                     cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      >  1:27903/1‑233 (MQ=255)
                                                                                                                                                                                                         ggTTTATTTTACTTTTTAGCAACCCTCACGTTACGTTATCGCCTGATGCCTAAGCTTTATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGNTTTGCTNTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGc                                                                                                                  <  3:27903/233‑1 (MQ=255)
                                                                                                                                                                                                           tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                                                                                                                                                           <  6:844040/126‑1 (MQ=255)
                                                                                                                                                                                                            atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  5:956572/1‑115 (MQ=255)
                                                                                                                                                                                                            atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                                                                                                                                    <  5:841085/148‑1 (MQ=255)
                                                                                                                                                                                                               tttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGACTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGt                                                                  <  2:67416/275‑1 (MQ=255)
                                                                                                                                                                                                                 tGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:728201/83‑1 (MQ=255)
                                                                                                                                                                                                                   aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 <  6:375631/80‑1 (MQ=255)
                                                                                                                                                                                                                   aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:176199/1‑207 (MQ=255)
                                                                                                                                                                                                                   aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                                                                                                                                                <  6:386811/81‑1 (MQ=38)
                                                                                                                                                                                                                     tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  8:218250/148‑1 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGc                                                                                                                                                                                                                                                                  >  5:912234/1‑76 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                                                                                                                                                           >  6:948350/1‑115 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          <  6:899195/148‑1 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          <  5:48998/148‑1 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          >  5:119894/1‑148 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGATCTATc                                                                                                                                                                                                                                                 <  6:606869/93‑1 (MQ=255)
                                                                                                                                                                                                                       tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                >  3:82877/1‑205 (MQ=255)
                                                                                                                                                                                                                       tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  4:176198/207‑1 (MQ=255)
                                                                                                                                                                                                                           gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  8:729935/1‑94 (MQ=255)
                                                                                                                                                                                                                           gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  7:729935/94‑1 (MQ=255)
                                                                                                                                                                                                                               cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTg                                                  >  2:515127/1‑275 (MQ=255)
                                                                                                                                                                                                                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGTTCTATCAGCTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTTCGGGGAACGGAACGGTACCCGGTGAGGTTAATGGTATGGAc                                                >  4:470624/1‑275 (MQ=255)
                                                                                                                                                                                                                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                                                                                                                                                              >  8:554530/1‑117 (MQ=255)
                                                                                                                                                                                                                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  5:576791/136‑1 (MQ=255)
                                                                                                                                                                                                                                     cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGgtgt                                                                                                                                                                                                                                                                           >  7:986353/1‑52 (MQ=255)
                                                                                                                                                                                                                                     cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                                                                                                                                           <  5:602983/148‑1 (MQ=255)
                                                                                                                                                                                                                                     cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                                                                                                               <  2:270806/160‑1 (MQ=255)
                                                                                                                                                                                                                                       ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                                                                                                                                                              <  5:385274/95‑1 (MQ=255)
                                                                                                                                                                                                                                         acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       <  6:735218/148‑1 (MQ=255)
                                                                                                                                                                                                                                             taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTg                                                                                                                                                                                                                                                        <  6:571992/63‑1 (MQ=255)
                                                                                                                                                                                                                                              aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                                                                                                                                                         >  8:772118/1‑93 (MQ=255)
                                                                                                                                                                                                                                                tgCCTGATGCCTAAGCTATATTTGGAAGCCTGTGATGGTGTAGACCAGGTACTGGATTTGCTCTATCAGTTCGAGTTTTGGTAAGGGGAAGATTTGTCTACGGCAGTTTAAGCCACAAAAAAGGCTACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCCGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGTGAACGTAACGGTACACGGTGAAGTTTTTCGTATTGACAACGCCACGCTGGcc                                 <  4:440370/274‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  4:548391/1‑94 (MQ=255)
                                                                                                                                                                                                                                                 gCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      <  5:686188/93‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAg                                                                                                                                                                                                                                   <  2:569108/79‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                                                                                                                    <  1:83116/142‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               >  5:862535/1‑129 (MQ=255)
                                                                                                                                                                                                                                                      aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACt                            >  4:134835/1‑274 (MQ=255)
                                                                                                                                                                                                                                                       tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       <  6:82370/134‑1 (MQ=255)
                                                                                                                                                                                                                                                           tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     <  7:774241/68‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                                                                                                                                                  <  8:801526/134‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                               >  4:523334/1‑232 (MQ=255)
                                                                                                                                                                                                                                                              gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  6:1112369/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:130657/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCa                                                                                                                                                                                      >  6:998951/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                 atatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                                                                                                                                          >  7:1061095/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                 atatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                                                                                                                                          <  8:1061095/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                   atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                                                                                                              <  7:422475/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                   atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATTACCAATGCCCAGTAACTGGGcatt                                                                                                                                                >  7:133464/1‑142 (MQ=255)
                                                                                                                                                                                                                                                                       ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    >  8:321116/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                       ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    <  7:321116/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                                                                                                                                                         >  6:31925/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                        gAAGCCGTGTCGGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           >  7:81599/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                         aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  8:586787/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                            ccGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGGACACGGTGAAGTTTATCGTATTGACAACGCCACCCTGGCCGAACTTGTGGCCTTTGCCACCAGGGGCg     >  4:458523/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  5:272109/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gtgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGca                                                                                                                                                  <  5:797122/132‑2 (MQ=255)
                                                                                                                                                                                                                                                                               tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCTACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  7:739957/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                               tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGTCGATTTCAGTATCGATAAc                                                                                                                                 >  8:860465/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  6:593717/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCAACTATCCAGGCGCCGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACCCCACGCATGCCGAAATTGATGCCTTTCGCCCCAGGGGCGGTg  >  4:344483/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                gtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCAc             <  2:458524/263‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                                                                                                                                                       >  7:512124/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
GAGAGCGATCAGAGCGACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447571‑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: