Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,314,800 A→C 5.3% Q438P (CAG→CCG) ‡ rcsD → phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,314,8000AC 5.3% 897.3 / 30.8 258Q438P (CAG→CCG) ‡rcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base A (130/114);  new base C (6/8);  total (136/122)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 261 total aligned reads displayed.

CGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATCAGGAATCGAACCGCACTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGTTAGCGAACATGCTTGCGGATGATAGCTGGAAAAGTGAGACGGTGCTGTTCTCCGTG  >  NC_000913/2314531‑2315074
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                  
cGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATCAGGAATCGAACCGCCCTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACCAGCTGGATGAGATCCGAATGTTCCGCAGCCAGGTCGCGCCGCGGCCACAAATTTTAATTATTAGCGATCAGGaac                                                                                                                                                                                                                                                                               >  4:520771/1‑273 (MQ=255)
                                 gaTAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATCAGGAATCGAACCGCACTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCCAAACATCACAACCATTGCGGAACAGCATCAGGGGATTATTAAGGCGACGGTCAATAACGCGCTGTCTGAGATCCGCATGTTCCGCATCACGGTCGGGCCGCCCAACCCAATTTTCATTATTAGCTACCAGGCTCGCGCAGTGCTGCGAACCAAACAACTCAAACAgg                                                                                                                                                                                                                                              >  4:75233/1‑275 (MQ=255)
                                     gTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATCAGGAATCGAACCGCACTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCGCCACCATGGCGGAACAGCATCAGAGGATTATTCAGGCGACGATAAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGGTCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGACGCa                                                                                                                                                                                                                                          <  2:309200/275‑1 (MQ=255)
                                               tgcCGCTCGGCCTGCTGGTTCACGATCAGGAATCGAACCGCACTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTa                                                                                                                                                                                                                                 <  2:691565/274‑1 (MQ=255)
                                                                gttacagagaaGGAATCGAACCGCACTGTCATAAGTAACAAAATTGCCGATCAGTTGCTGACGCATTTGAATCTGAAAAACATCACCACCATGGCGGAACAGCAGCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAgg                                                                                                                                                                                                                <  4:192763/264‑1 (MQ=255)
                                                                                       cACTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTCTGAGATCCGCATGTTCCGCAGCCAGGTCGCTCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAAAAAATCAAGCAGGCGCCGCGTCGGTATGAGAAAACCCAACAGGGAGCGATGATCTTTACTaaaaaac                                                                                                                                                                                        >  4:155688/1‑273 (MQ=255)
                                                                                                      cAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcgc                                                                                                                                                                                                                                                                             <  3:146767/175‑1 (MQ=255)
                                                                                                      cAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCCATGCGCTGa                                                                                                                                                                         >  2:183902/1‑275 (MQ=255)
                                                                                                      cAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGGATGAGAAAAACCAGCAGGGGCGGATGATCGTTATGAAAAACATTGGCAGTGCGCTGa                                                                                                                                                                         >  4:361113/1‑275 (MQ=255)
                                                                                                          aTTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaa                                                                                                                                                                                                                                                                <  4:683473/184‑1 (MQ=255)
                                                                                                                                 tttGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcgcg                                                                                                                                                                                                                                                                            >  7:913964/1‑149 (MQ=255)
                                                                                                                                 tttGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATc                                                                                                                                                                                                                                                                               <  8:64810/146‑1 (MQ=255)
                                                                                                                                 tttGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATc                                                                                                                                                                                                                                                                               >  7:64810/1‑146 (MQ=255)
                                                                                                                                   tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaa                                                                                                                                                                                                                                                                          <  8:35560/149‑1 (MQ=255)
                                                                                                                                    gAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAg                                                                                                                                                                                                                                                                         >  5:1708422/1‑149 (MQ=255)
                                                                                                                                      aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTg                                                                                                                                                                                                                                                                       >  5:921821/1‑149 (MQ=255)
                                                                                                                                         tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGCAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAACGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCt                                                                                                                                      >  2:270868/1‑275 (MQ=255)
                                                                                                                                         tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                                                                                                                                                    <  7:1605230/149‑1 (MQ=255)
                                                                                                                                         tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                                                                                                                                                    <  7:345329/149‑1 (MQ=255)
                                                                                                                                          gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGc                                                                                                                                                                                                                                                                      >  2:579030/1‑146 (MQ=255)
                                                                                                                                              aaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaaa                                                                                                                                                                                                                                                               >  5:330944/1‑149 (MQ=255)
                                                                                                                                               aCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAAc                                                                                                                                                                                                                                                              <  5:127336/149‑1 (MQ=255)
                                                                                                                                                  tCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAg                                                                                                                                                                                                                                                           >  7:415053/1‑149 (MQ=255)
                                                                                                                                                    accaccATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaata                                                                                                                                                                                                                                                                        <  7:469497/134‑3 (MQ=255)
                                                                                                                                                     ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                                                                                                                                                        <  8:415053/149‑1 (MQ=255)
                                                                                                                                                     ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                                                                                                                                                        <  6:126090/149‑1 (MQ=255)
                                                                                                                                                        ccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTc                                                                                                                                                                                                                                                     <  5:15744/149‑1 (MQ=255)
                                                                                                                                                          aTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCaa                                                                                                                                                                                                                                                   <  8:581851/149‑1 (MQ=255)
                                                                                                                                                          aTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCaa                                                                                                                                                                                                                                                   <  8:708800/149‑1 (MQ=255)
                                                                                                                                                          aTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGa                                                                                                                     <  2:520775/275‑1 (MQ=255)
                                                                                                                                                                 aaCAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAg                                                                                                                                                                                                                                                                                   >  8:1011022/1‑110 (MQ=255)
                                                                                                                                                                     gCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGc                                                                                                                                                                                                                                        >  8:1317647/1‑149 (MQ=255)
                                                                                                                                                                         cAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAACCTCAAGCAGGCGCAGCGTCt                                                                                                                                                                                                                                    <  6:449045/149‑1 (MQ=255)
                                                                                                                                                                            gggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTa                                                                                                                                                                                                                                 <  7:73404/149‑1 (MQ=255)
                                                                                                                                                                            gggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTa                                                                                                                                                                                                                                 <  5:863489/149‑1 (MQ=255)
                                                                                                                                                                              gATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                                                                                                                                                        <  7:1603333/124‑1 (MQ=255)
                                                                                                                                                                                ttattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                                                                                                                                                    >  7:1030838/1‑110 (MQ=255)
                                                                                                                                                                                ttattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAACGCAAACAACTGGCGAATCAg                                                                                               >  2:672645/1‑275 (MQ=255)
                                                                                                                                                                                 tattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATgaga                                                                                                                                                                                                                            <  5:482521/149‑1 (MQ=255)
                                                                                                                                                                                   ttCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaa                                                                                                                                                                                                                          >  6:1497532/1‑149 (MQ=255)
                                                                                                                                                                                      aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCTGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTAt                                                                                                                                                                                                                                >  6:1583785/1‑140 (MQ=255)
                                                                                                                                                                                      aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                                                                                                                                                       >  7:1529235/1‑149 (MQ=255)
                                                                                                                                                                                        gCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                                                                                                                                                    <  7:553569/102‑1 (MQ=255)
                                                                                                                                                                                        gCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                                                                                                                                                     <  6:1190448/149‑1 (MQ=255)
                                                                                                                                                                                        gCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                                                                                                                                                     >  7:1325048/1‑149 (MQ=255)
                                                                                                                                                                                          gacgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACcagc                                                                                                                                                                                                                   >  8:1343976/1‑149 (MQ=255)
                                                                                                                                                                                              aTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAggg                                                                                                                                                                                                               <  5:274949/149‑1 (MQ=255)
                                                                                                                                                                                              aTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAggg                                                                                                                                                                                                               >  8:602359/1‑149 (MQ=255)
                                                                                                                                                                                                cAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCa                                                                                                                                                                                                                                                                                    >  6:1676906/1‑78 (MQ=255)
                                                                                                                                                                                                cAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGc                                                                                                                                                                                                             <  6:730691/149‑1 (MQ=255)
                                                                                                                                                                                                     aCGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTTAAGCAGGCGCAGCGTCTGTATGAGAAAATCCAGCAGGGGCGgatg                                                                                                                                                                                                        >  8:263443/1‑149 (MQ=255)
                                                                                                                                                                                                       gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACcag                                                                                                                                                                                                                    >  6:486093/1‑135 (MQ=255)
                                                                                                                                                                                                       gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatgat                                                                                                                                                                                                      >  8:739873/1‑149 (MQ=255)
                                                                                                                                                                                                       gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatgat                                                                                                                                                                                                      <  8:1088834/149‑1 (MQ=255)
                                                                                                                                                                                                         gCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTCGTAAACAAGAAACTCAAGCAGGCTCAGCGTCGGTCTGAGCAAAACCAGCAGGGNCGGATGATCTTTATGAAAAACATTGGAGATGCGCCGAAAGAACCCGCAAAGTCCCCGGCGGAGAGCAGGGAGAAACTAAACGCACCGGAAAGCAAACAACTggtgtnannagta                                                                                            >  3:190419/1‑242 (MQ=255)
                                                                                                                                                                                                          cTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCtt                                                                                                                                                                                                   <  7:18980/149‑1 (MQ=255)
                                                                                                                                                                                                           tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAgcgc                                                                                                                                          >  4:23138/1‑205 (MQ=255)
                                                                                                                                                                                                            gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                                                                                                                                                 <  6:984315/149‑1 (MQ=255)
                                                                                                                                                                                                            gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                                                                                                                                                 <  6:1088486/149‑1 (MQ=255)
                                                                                                                                                                                                             tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATTATCTTTAt                                                                                                                                                                                                <  6:558639/149‑1 (MQ=255)
                                                                                                                                                                                                             tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTAt                                                                                                                                                                                                <  6:821168/149‑1 (MQ=255)
                                                                                                                                                                                                                 agaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCttt                                                                                                                                                                                                  >  6:967686/1‑143 (MQ=255)
                                                                                                                                                                                                                 agaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGaaa                                                                                                                                                                                            <  6:921821/149‑1 (MQ=255)
                                                                                                                                                                                                                 agaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCg                                                                                                                                                                                  <  4:654328/159‑1 (MQ=255)
                                                                                                                                                                                                                   aTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                                                                                                                                                     >  5:1367901/1‑138 (MQ=255)
                                                                                                                                                                                                                   aTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                                                                                                                                                     <  6:1367901/138‑1 (MQ=255)
                                                                                                                                                                                                                     ccGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                                                                                                                                                        <  8:1642108/85‑1 (MQ=255)
                                                                                                                                                                                                                     ccGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGGTGAAATACCGt                                                          >  2:229294/1‑275 (MQ=255)
                                                                                                                                                                                                                     ccGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGt                                                          <  2:155690/275‑1 (MQ=255)
                                                                                                                                                                                                                       gCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCt                                                                                                                                                                                                                                                                     >  2:330629/1‑70 (MQ=255)
                                                                                                                                                                                                                       gCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGt                                                                                                                                                                                                                                      >  2:428603/1‑101 (MQ=255)
                                                                                                                                                                                                                       gCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACcagc                                                                                                                                                                                                                   >  1:183671/1‑120 (MQ=255)
                                                                                                                                                                                                                       gCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGt                                                                                                                                                          <  1:277661/177‑1 (MQ=255)
                                                                                                                                                                                                                         aTGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATTAAAAACATTgg                                                                                                                                                                                    <  7:985885/149‑1 (MQ=255)
                                                                                                                                                                                                                          tGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGc                                                                                                                                                                                   <  7:1260284/149‑1 (MQ=255)
                                                                                                                                                                                                                            tcccccAGCCGGGCCGCGCCGAGCAATCAAATTTGCATGATTCGCGATCAGGATCGGGAAGGGCTGGTAAACAAGAAACCCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGTTAGCGaa                                                   <  4:477385/270‑1 (MQ=255)
                                                                                                                                                                                                                            gaccGCAGCCAGAGCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGTTAGCGAACAt                                                <  1:102563/276‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                                                                                                                                                              <  7:1108533/95‑1 (MQ=255)
                                                                                                                                                                                                                                      aGGTCGCGCCGCGCACACAAATTTTCATTATTCTCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatgat                                                                                                                                                                                                      >  5:524475/1‑118 (MQ=255)
                                                                                                                                                                                                                                      aGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGaaa                                                                                                                                                                       >  5:681508/1‑149 (MQ=255)
                                                                                                                                                                                                                                      aGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGaaa                                                                                                                                                                       <  5:1497532/149‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTAt                                                                                                                                                                                                                                >  5:635460/1‑91 (MQ=255)
                                                                                                                                                                                                                                         tCGCGCCGCGCAGACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGaa                                                                                                                                                                    <  5:779917/149‑1 (MQ=255)
                                                                                                                                                                                                                                             gccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                                                                                                                                                                >  5:229405/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  gcccacaAATTTTCATTATTCGCGATCAGGATCGCGACGTCCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTCTGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCAGCACAg                                                                                                                                                           >  5:1625741/4‑149 (MQ=255)
                                                                                                                                                                                                                                                  gcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGAACTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGTTAGCGAACATGCTTGCGGATGATAGCTgg                             >  2:645926/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   cacacaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGt                                                                                                                                                          <  6:668772/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    acacaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTc                                                                                                                                                         <  5:950380/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACTCGCACAGTccc                                                                                                                                                       >  8:1485636/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                                                                                                                                                       >  8:911485/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                                                                                                                                                       >  5:374926/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aaaTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCt                                                                                                                                                      >  7:277738/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        aaaTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg                                                                                                                                                     >  7:811910/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aaTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTgg                                                                                                                                                    <  8:169130/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          aTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGc                                                                                                                                                   >  6:373678/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                attattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGagag                                                                                                                                             >  7:765929/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                attattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTg                                                                          >  2:727660/1‑216 (MQ=255)
                                                                                                                                                                                                                                                                    ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                                                                                                                                                     >  8:1184929/1‑73 (MQ=39)
                                                                                                                                                                                                                                                                    ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                                                                                                                                                     <  7:1184929/73‑1 (MQ=39)
                                                                                                                                                                                                                                                                            cAGGATCGCGAAGGGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGTTAGCGAACATGCTTGCGGATGATAGCTGGAAAAGTGAGACGGTGCTGTTCTCCGt   <  2:656446/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                  
CGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATCAGGAATCGAACCGCACTGTCATAAGTAACAAAATTGCCGATCATTTGCTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAACGCCCCGGAAAGCAAACAACTGGCGAATCAGGCAGATGTGCTGGTGCGGCTGGTTGATGAAATACAGTTAGCGAACATGCTTGCGGATGATAGCTGGAAAAGTGAGACGGTGCTGTTCTCCGTG  >  NC_000913/2314531‑2315074

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: