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.

GAGAGCGATCAGAGCGACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  NC_000913/4447571‑4448116
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
gagaGCGATCAGAGCGACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCgtgt                                                                                                                                                                                                                                                                                   >  2:226970/1‑274 (MQ=255)
                        gCAATGACCTCGATGCTGACTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTCAAATCGGTGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGCAATACGGCGTGGGTATATTTGACTCTTTAGCAACACTCACGTTACGTTAACGCCTTATGCCTAAGCTATTACGGGAAGCCGTTT‑CGGGTGTAAACCAGGCCCTGAAttg                                                                                                                                                                                                                                                         >  4:75509/1‑274 (MQ=255)
                                                        tGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCACAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGCATATTTGACTATATAGAAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGACGTGT‑CAGATGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGAAGTGCGAat                                                                                                                                                                                                                         >  2:536975/1‑275 (MQ=255)
                                                                                tAAAATCGGCGAGACGTGTGGCATAAGCAATTTAGCGCTCGACACCCAAGGAGTAGGCAACTCCTCCCAGGTAGTGGTCAGAGTCTATGTATTACCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAACCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                                                                                                                  <  1:72699/274‑1 (MQ=255)
                                                                                               gTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACa                                                                                                                                                                                   <  2:378567/274‑1 (MQ=255)
                                                                                                                                gggAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  2:278990/1‑206 (MQ=255)
                                                                                                                                   agggggggACTCCTCCCAGGTAGTGGTCGGCGGCTAGGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATAGTTGACTCTATAGCAACACTCACGTTACGTTATAGGCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                                                                                                                              <  4:132224/268‑1 (MQ=255)
                                                                                                                                    gTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTc                                                                                                                                                                                                                  <  4:278989/206‑1 (MQ=255)
                                                                                                                                        gCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTCGCAATGCGAATATTTGTCTACGGCAGTTTCCGCCACAAACCAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                                                                                                                  >  2:132225/1‑266 (MQ=255)
                                                                                                                                         cGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGAc                                                                                                                                                                                                                                                                       >  6:692861/1‑148 (MQ=255)
                                                                                                                                         cGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGAc                                                                                                                                                                                                                                                                       >  6:702346/1‑148 (MQ=255)
                                                                                                                                             tcctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGTAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATAGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGGNTAGACCAGGCACTGGATGTGCTCTATCAGTTCGAGTTTTAGCAATGGGAATATGTGTCTACGGAGGTTTAGGACACAAAANAGGCAACAGTCACTAGATTACCAACNCCAAGTTAATGGGAGCTTTCAGTaaagaa                                                                                                                                    >  3:94941/1‑270 (MQ=255)
                                                                                                                                               ctccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  4:226969/274‑1 (MQ=255)
                                                                                                                                                 cccAGGTAGTGGTCAGCGGCTACGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAAc                                                                                                                                 >  2:298251/1‑274 (MQ=255)
                                                                                                                                                   ctggTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTg                                                                                                                                                                                                                                                                               <  7:572680/128‑1 (MQ=255)
                                                                                                                                                   cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  5:323178/144‑1 (MQ=255)
                                                                                                                                                   cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGGGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                <  1:117083/273‑1 (MQ=255)
                                                                                                                                                    aGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGa                                                                                                                                                                                                                                                            >  8:914379/1‑148 (MQ=255)
                                                                                                                                                       tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGacac                                                                                                                                                                                                                                                                >  8:423507/1‑141 (MQ=255)
                                                                                                                                                          tGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:430550/1‑264 (MQ=255)
                                                                                                                                                           ggTCAGCGGTTATGTATTGCCAGGTCTGCAAGTGAAATCCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      <  2:497033/275‑1 (MQ=255)
                                                                                                                                                            gTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  2:981/1‑163 (MQ=255)
                                                                                                                                                               aGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAg                                                                                                                   >  4:202437/1‑274 (MQ=255)
                                                                                                                                                                gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCa                                                                                                                                                                                                                                 <  4:981/163‑1 (MQ=255)
                                                                                                                                                                gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATTGcc                                                                                                                 >  2:138523/1‑275 (MQ=255)
                                                                                                                                                                 cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAgac                                                                                                                                                                                                                                                                  <  5:698403/129‑3 (MQ=255)
                                                                                                                                                                        gTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                                    <  4:308680/248‑1 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  7:480045/1‑141 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  8:480045/141‑1 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 <  7:1123768/119‑1 (MQ=255)
                                                                                                                                                                            tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGACa                                                                                                                                                                                                                                                                 >  8:1123768/1‑119 (MQ=255)
                                                                                                                                                                                   tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATACCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                              >  2:116152/1‑275 (MQ=255)
                                                                                                                                                                                     tGCAAGTGAAATACGGCGCGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGttt                                                                                                                                                                                                       <  2:216459/168‑1 (MQ=255)
                                                                                                                                                                                          gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCa                                                                                                                                                                                                                                                                     <  8:200026/101‑1 (MQ=255)
                                                                                                                                                                                          gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  7:218250/1‑148 (MQ=255)
                                                                                                                                                                                                 aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                                                                                                                                                                                                                             <  7:908533/118‑1 (MQ=255)
                                                                                                                                                                                                     cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                                                                                                                      >  1:27903/1‑233 (MQ=255)
                                                                                                                                                                                                         ggTTTATTTTACTTTTTAGCAACCCTCACGTTACGTTATCGCCTGATGCCTAAGCTTTATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGNTTTGCTNTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGc                                                                                                                  <  3:27903/233‑1 (MQ=255)
                                                                                                                                                                                                           tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                                                                                                                                                           <  6:844040/126‑1 (MQ=255)
                                                                                                                                                                                                            atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     >  5:956572/1‑115 (MQ=255)
                                                                                                                                                                                                            atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                                                                                                                                    <  5:841085/148‑1 (MQ=255)
                                                                                                                                                                                                               tttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGACTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGt                                                                  <  2:67416/275‑1 (MQ=255)
                                                                                                                                                                                                                 tGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCAc                                                                                                                                                                                                                                                                <  5:728201/83‑1 (MQ=255)
                                                                                                                                                                                                                   aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCa                                                                                                                                                                                                                                                                 <  6:375631/80‑1 (MQ=255)
                                                                                                                                                                                                                   aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                                  >  2:176199/1‑207 (MQ=255)
                                                                                                                                                                                                                   aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGC‑CTGGTCTACACCAGacac                                                                                                                                                                                                                                                                <  6:386811/81‑1 (MQ=38)
                                                                                                                                                                                                                     tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  8:218250/148‑1 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGc                                                                                                                                                                                                                                                                  >  5:912234/1‑76 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                                                                                                                                                           >  6:948350/1‑115 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          <  6:899195/148‑1 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          <  5:48998/148‑1 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                                                                                                                                          >  5:119894/1‑148 (MQ=255)
                                                                                                                                                                                                                      ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGATCTATc                                                                                                                                                                                                                                                 <  6:606869/93‑1 (MQ=255)
                                                                                                                                                                                                                       tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                                >  3:82877/1‑205 (MQ=255)
                                                                                                                                                                                                                       tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                              <  4:176198/207‑1 (MQ=255)
                                                                                                                                                                                                                           gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  8:729935/1‑94 (MQ=255)
                                                                                                                                                                                                                           gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           <  7:729935/94‑1 (MQ=255)
                                                                                                                                                                                                                               cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTg                                                  >  2:515127/1‑275 (MQ=255)
                                                                                                                                                                                                                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGTTCTATCAGCTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTTCGGGGAACGGAACGGTACCCGGTGAGGTTAATGGTATGGAc                                                >  4:470624/1‑275 (MQ=255)
                                                                                                                                                                                                                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                                                                                                                                                              >  8:554530/1‑117 (MQ=255)
                                                                                                                                                                                                                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           <  5:576791/136‑1 (MQ=255)
                                                                                                                                                                                                                                     cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGgtgt                                                                                                                                                                                                                                                                           >  7:986353/1‑52 (MQ=255)
                                                                                                                                                                                                                                     cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                                                                                                                                           <  5:602983/148‑1 (MQ=255)
                                                                                                                                                                                                                                     cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                                                                                                               <  2:270806/160‑1 (MQ=255)
                                                                                                                                                                                                                                       ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                                                                                                                                                              <  5:385274/95‑1 (MQ=255)
                                                                                                                                                                                                                                         acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       <  6:735218/148‑1 (MQ=255)
                                                                                                                                                                                                                                             taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTg                                                                                                                                                                                                                                                        <  6:571992/63‑1 (MQ=255)
                                                                                                                                                                                                                                              aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                                                                                                                                                         >  8:772118/1‑93 (MQ=255)
                                                                                                                                                                                                                                                tgCCTGATGCCTAAGCTATATTTGGAAGCC‑TGTGATGGTGTAGACCAGGTACTGGATTTGCTCTATCAGTTCGAGTTTTGGTAAGGGGAAGATTTGTCTACGGCAGTTTAAGCCACAAAAAAGGCTACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCCGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGTGAACGTAACGGTACACGGTGAAGTTTTTCGTATTGACAACGCCACGCTGGcc                                 <  4:440370/274‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      >  4:548391/1‑94 (MQ=255)
                                                                                                                                                                                                                                                 gCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                                                                                                                                                      <  5:686188/93‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAg                                                                                                                                                                                                                                   <  2:569108/79‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                                                                                                                    <  1:83116/142‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                                                                                                               >  5:862535/1‑129 (MQ=255)
                                                                                                                                                                                                                                                      aTGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACt                            >  4:134835/1‑274 (MQ=255)
                                                                                                                                                                                                                                                       tGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       <  6:82370/134‑1 (MQ=255)
                                                                                                                                                                                                                                                           tAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                                                                                                                                                     <  7:774241/68‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                                                                                                                                                  <  8:801526/134‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                               >  4:523334/1‑232 (MQ=255)
                                                                                                                                                                                                                                                              gCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                                                                                                                                                           >  6:1112369/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                tataTCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  6:130657/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                tataTCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCa                                                                                                                                                                                      >  6:998951/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                 atatCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                                                                                                                                          >  7:1061095/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                 atatCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                                                                                                                                          <  8:1061095/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                   atCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                                                                                                              <  7:422475/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                   atCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATTACCAATGCCCAGTAACTGGGcatt                                                                                                                                                >  7:133464/1‑142 (MQ=255)
                                                                                                                                                                                                                                                                       ggAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    >  8:321116/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                       ggAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                                                                                                                    <  7:321116/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                                                                                                                                                         >  6:31925/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                        gAAGCCGTGT‑CGGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                                                                                                                           >  7:81599/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                         aaGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                                                                                                                                                >  8:586787/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                            ccGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGGACACGGTGAAGTTTATCGTATTGACAACGCCACCCTGGCCGAACTTGTGGCCTTTGCCACCAGGGGCg     >  4:458523/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  5:272109/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gtgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGca                                                                                                                                                  <  5:797122/132‑2 (MQ=255)
                                                                                                                                                                                                                                                                               tgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCTACAAACAAGGCAACAGTCACTGGATGa                                                                                                                                                                       >  7:739957/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                               tgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGTCGATTTCAGTATCGATAAc                                                                                                                                 >  8:860465/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               tgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                                   <  6:593717/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tgt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCAACTATCCAGGCGCCGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACCCCACGCATGCCGAAATTGATGCCTTTCGCCCCAGGGGCGGTg  >  4:344483/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                gt‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCAc             <  2:458524/263‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 t‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                                                                                                                                                       >  7:512124/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
GAGAGCGATCAGAGCGACAGTGCGGCAATGACCTCGATGCTGATTGGTTTGGGGGTTGCGCAAAGTGGCCAGATTGTGGGTAAAATCGGCGAGACGTTTGGCGTAAGCAATTTAGCGCTCGACACCCAGGGAGTAGGCGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGT‑CTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTG  >  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: ‑