Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,314,800 A→C 8.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 8.3% 623.6 / 38.6 188Q438P (CAG→CCG) ‡rcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base A (80/92);  new base C (9/7);  total (89/99)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 188 total aligned reads displayed.

TGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAAC  >  NC_000913/2314653‑2314946
                                                                                                                                                   |                                                                                                                                                  
tgcCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAg                                                                                                                                                   <  2:12437/149‑1 (MQ=255)
    gCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATc                                                                                                                                               <  3:1350998/149‑1 (MQ=255)
      aTTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGACCAGGATcgc                                                                                                                                             <  2:45320/149‑1 (MQ=255)
       tttGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCGGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcgcg                                                                                                                                            <  4:772745/149‑1 (MQ=255)
         tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                          <  4:264044/149‑1 (MQ=255)
         tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                          <  4:652753/149‑1 (MQ=255)
         tGAATCTGCAAAACATCACCACCATGGCGCAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaa                                                                                                                                          <  1:44199/149‑1 (MQ=255)
            aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAGTTTTCATTATTCGCGATCCTGATCGCGaataa                                                                                                                                       <  3:1006839/149‑4 (MQ=255)
            aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTg                                                                                                                                       <  1:1073775/149‑1 (MQ=255)
               tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATTCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACATAACTGTTCATTATTTGCCATCCGGATCGCTaattcctt                                                                                                                                    <  2:361479/149‑7 (MQ=255)
               tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCCTTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                    <  2:1481767/149‑1 (MQ=255)
                gCAAAACATCGCCACCATGGCGGAACAGCATGAGGGGATTATTCAGGCGACGATCAATAACGAGCTATATGAGATACGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAATTTTTCATTATTCGCGATCAGGATCGCGGAGTGCTgg                                                                                                                                   <  4:343151/149‑1 (MQ=255)
                gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTgg                                                                                                                                   <  1:84712/149‑1 (MQ=255)
                  aaaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTa                                                                                                                                 >  4:1545303/1‑149 (MQ=255)
                   aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                          <  3:124298/139‑1 (MQ=255)
                   aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAAATTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                          >  4:124298/1‑139 (MQ=255)
                    aaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaaa                                                                                                                               <  2:1100359/149‑1 (MQ=255)
                       aTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                            >  1:1136746/1‑149 (MQ=255)
                        tCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAg                                                                                                                           <  3:1542431/149‑1 (MQ=255)
                        tCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAg                                                                                                                           >  1:317964/1‑149 (MQ=255)
                         caccacCATGGTGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGa                                                                                                                          <  3:1017585/149‑1 (MQ=255)
                          accaccATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaa                                                                                                                         >  4:321488/1‑149 (MQ=255)
                                aTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                          >  4:484779/1‑126 (MQ=255)
                                aTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCaa                                                                                                                   >  4:728258/1‑149 (MQ=255)
                                      gAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGgc                                                                                                             <  4:502122/149‑1 (MQ=255)
                                        aCAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGgcgc                                                                                                           >  2:140634/1‑149 (MQ=255)
                                        aCAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCACATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGgcgc                                                                                                           >  2:626493/1‑149 (MQ=255)
                                             aTCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCTGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGt                                                                                                      <  2:1136746/149‑1 (MQ=255)
                                              tCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                              <  1:407613/140‑1 (MQ=255)
                                              tCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                              >  2:407613/1‑140 (MQ=255)
                                              tCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTc                                                                                                     >  4:856385/1‑149 (MQ=255)
                                               cAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCt                                                                                                    <  4:690621/149‑1 (MQ=255)
                                                 ggggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGt                                                                                                  >  4:1292182/1‑149 (MQ=255)
                                                 ggggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGt                                                                                                  >  1:1426066/1‑149 (MQ=255)
                                                 ggggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGt                                                                                                  <  4:374530/149‑1 (MQ=255)
                                                 ggggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGt                                                                                                  <  2:15998/149‑1 (MQ=255)
                                                    gATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGa                                                                                                                          <  1:192417/122‑1 (MQ=255)
                                                     attattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                              <  4:97269/149‑1 (MQ=255)
                                                      ttattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATgag                                                                                             >  3:630247/1‑149 (MQ=255)
                                                      ttattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGACGCTGCGCCTGCTTgag                                                                                             <  1:857188/149‑1 (MQ=255)
                                                         ttCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaa                                                                                          >  4:397381/1‑149 (MQ=255)
                                                            aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                       >  2:1102441/1‑149 (MQ=255)
                                                            aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                       <  1:1229210/149‑1 (MQ=255)
                                                            aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                       <  2:859474/149‑1 (MQ=255)
                                                                 acgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACcagca                                                                                  <  3:397381/149‑1 (MQ=255)
                                                                   gaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAgg                                                                                <  2:945338/149‑1 (MQ=255)
                                                                    aTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAggg                                                                               >  3:348314/1‑149 (MQ=255)
                                                                         tAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGga                                                                          >  3:52330/1‑149 (MQ=255)
                                                                         tAACGAGCTGTATGAGATCCACATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGga                                                                          <  1:626493/149‑1 (MQ=255)
                                                                          aaCGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgat                                                                         <  1:909208/149‑1 (MQ=255)
                                                                             gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaat                                                                                                                                         >  1:965039/1‑81 (MQ=255)
                                                                             gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaat                                                                                                                                         <  2:965039/82‑2 (MQ=255)
                                                                             gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatgat                                                                      <  4:697853/149‑1 (MQ=255)
                                                                                  gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                 <  4:151406/149‑1 (MQ=255)
                                                                                  gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                 <  3:856385/149‑1 (MQ=255)
                                                                                  gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                 <  2:1014936/149‑1 (MQ=255)
                                                                                   tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTAt                                                                <  3:1314/149‑1 (MQ=255)
                                                                                   tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTAt                                                                >  3:307350/1‑149 (MQ=255)
                                                                                    aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAg                                                                                                               <  1:219085/101‑1 (MQ=255)
                                                                                    aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAg                                                                                                               >  2:219085/1‑101 (MQ=255)
                                                                                        gaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                              >  1:586943/1‑114 (MQ=255)
                                                                                        gaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                              <  2:586943/114‑1 (MQ=255)
                                                                                         aTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACATAAATTTGCATACTTCGCGATCAGGATCGCGAAATGCTGGTAAACAAGAAACTCAAGCAGGCTGAGCGTCTGTATGACAAAGACCTGCAGTCCCGGATGATCTCTAttccaat                                                          <  2:124674/149‑7 (MQ=255)
                                                                                         aTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACGCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGaaaaa                                                          <  2:165288/149‑1 (MQ=255)
                                                                                         aTCCGCATGTTCCGCAGCAAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGaaaaa                                                          >  4:953305/1‑149 (MQ=255)
                                                                                             gCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACAtt                                                      >  1:1115816/1‑149 (MQ=255)
                                                                                               aTGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTgg                                                    >  4:931064/1‑149 (MQ=255)
                                                                                                 gTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTa                                                                                                 >  2:1034860/1‑102 (MQ=255)
                                                                                                  ttCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGa                                                 <  2:268196/149‑1 (MQ=255)
                                                                                                  ttCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGa                                                 <  1:1101444/149‑1 (MQ=255)
                                                                                                  ttCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGa                                                 <  1:134720/149‑1 (MQ=255)
                                                                                                    ccGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATg                                               >  1:1449359/1‑149 (MQ=255)
                                                                                                      gCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATgcg                                             >  4:545129/1‑149 (MQ=255)
                                                                                                             ggTCGCGCCGCGCACACAAATTTTCATGATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAg                                      <  2:1325882/149‑1 (MQ=255)
                                                                                                               tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGaa                                    <  2:533930/149‑1 (MQ=255)
                                                                                                               tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGaa                                    <  1:953998/149‑1 (MQ=255)
                                                                                                                   gccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                                >  1:1015155/1‑149 (MQ=255)
                                                                                                                            acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                       <  4:1481539/149‑1 (MQ=255)
                                                                                                                            acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                       >  4:955772/1‑149 (MQ=255)
                                                                                                                            acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                       >  4:280893/1‑149 (MQ=255)
                                                                                                                            acaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                       >  1:484761/1‑149 (MQ=255)
                                                                                                                              aaaTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg                     >  3:1169609/1‑149 (MQ=255)
                                                                                                                              aaaTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg                     <  4:1062571/149‑1 (MQ=255)
                                                                                                                                aTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGca                              >  1:1572756/1‑138 (MQ=255)
                                                                                                                                 ttttCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCg                  >  3:986727/1‑149 (MQ=255)
                                                                                                                                  tttCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCgg                 >  3:559164/1‑149 (MQ=255)
                                                                                                                                     cATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGgaga              <  2:19972/149‑1 (MQ=255)
                                                                                                                                     cATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGgaga              <  2:447239/149‑1 (MQ=255)
                                                                                                                                       ttattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGgaga              <  2:1190708/147‑1 (MQ=255)
                                                                                                                                       ttattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGgaga              >  1:1190708/1‑147 (MQ=255)
                                                                                                                                       ttattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAgc            >  1:1400107/1‑149 (MQ=255)
                                                                                                                                          ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCt                                           <  2:1539344/115‑1 (MQ=255)
                                                                                                                                          ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCt                                           >  1:1539344/1‑115 (MQ=255)
                                                                                                                                          ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg         >  2:1573415/1‑149 (MQ=255)
                                                                                                                                          ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg         >  1:1000808/1‑149 (MQ=255)
                                                                                                                                          ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg         >  2:537421/1‑149 (MQ=255)
                                                                                                                                          ttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg         >  4:1261095/1‑149 (MQ=255)
                                                                                                                                          ttCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg         >  3:741890/1‑149 (MQ=255)
                                                                                                                                                aTCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTaaa   <  2:807498/149‑1 (MQ=255)
                                                                                                                                                 tCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAAc  >  2:190644/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAAC  >  NC_000913/2314653‑2314946

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: