Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,314,801 G→T 11.3% Q438H (CAG→CAT) ‡ 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,8010GT11.3% 305.6 / 25.4 114Q438H (CAG→CAT) ‡rcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base G (45/56);  new base T (7/6);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

Only 100 of 117 total aligned reads displayed.

CTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAA  >  NC_000913/2314652‑2314950
                                                                                                                                                     |                                                                                                                                                     
ctgcCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAg                                                                                                                                                       >  2:1415950/1‑150 (MQ=255)
    cGCATTTGAATCTGCAAAACCTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAAAAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGGCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATc                                                                                                                                                   <  2:288996/150‑1 (MQ=255)
     gCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcg                                                                                                                                                  >  1:884760/1‑150 (MQ=255)
          tGAATCTTCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaat                                                                                                                                             <  2:324831/150‑2 (MQ=255)
          tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaat                                                                                                                                             <  1:1549402/150‑2 (MQ=255)
          tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaat                                                                                                                                             <  2:1323188/150‑2 (MQ=255)
          tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAg                                                                                                                                             >  1:1718821/1‑150 (MQ=255)
          tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGAGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAg                                                                                                                                             >  1:345359/1‑150 (MQ=255)
          ggAATCTGCAAAACATCACCCCCATGGCGGAACAGCATCAGGGGATTATTCCGGCGACGAGCACTACCGAGCTGTATGAGCTCCGCATGGTCCGCAGCCAGGTCGCGCCGCGCACACAAAGTTTCATTATTCGCGATCAGGATCGCGaat                                                                                                                                             <  1:214622/149‑2 (MQ=255)
             aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTTTTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGc                                                                                                                                          <  1:1375601/150‑1 (MQ=255)
              tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaataatg                                                                                                                                         <  1:1121357/150‑6 (MQ=255)
              tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCt                                                                                                                                         <  1:1676723/150‑1 (MQ=255)
              tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTGCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaataatg                                                                                                                                         <  2:467919/150‑6 (MQ=255)
               cTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGCTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACTCAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                        <  1:1292880/150‑1 (MQ=255)
               cTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACCAGTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                        >  2:210820/1‑150 (MQ=255)
               cTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTg                                                                                                                                        >  2:216092/1‑150 (MQ=255)
                tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTgg                                                                                                                                       <  2:1142659/150‑1 (MQ=255)
                tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTgg                                                                                                                                       >  1:504991/1‑150 (MQ=255)
                 gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGt                                                                                                                                      <  2:1410311/150‑1 (MQ=255)
                 gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGACAGGTCGCGCCGCGAACACAAATATTCATTATTAGCGATCAGGATCGCGAAGTGCTGGt                                                                                                                                      >  1:841122/1‑150 (MQ=255)
                  cAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTa                                                                                                                                     >  1:920368/1‑150 (MQ=255)
                  cAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTa                                                                                                                                     >  2:478934/1‑150 (MQ=255)
                   aaaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaa                                                                                                                                    >  2:29128/1‑150 (MQ=255)
                   aaaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGAGCAGGATCGCGAAGTGCTGGTaa                                                                                                                                    >  1:1006137/1‑150 (MQ=255)
                    aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTgg                                                                                                                                       <  2:1213093/146‑1 (MQ=255)
                    aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTgg                                                                                                                                       >  1:1213093/1‑146 (MQ=255)
                    aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaaa                                                                                                                                   >  2:880619/1‑150 (MQ=255)
                       cATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                                >  2:474994/1‑150 (MQ=255)
                        aTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAg                                                                                                                               <  1:1275736/150‑1 (MQ=255)
                         tCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGa                                                                                                                              <  2:1155671/150‑1 (MQ=255)
                            ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAAc                                                                                                                           <  1:438383/150‑1 (MQ=255)
                               ccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                                        >  2:806115/1‑150 (MQ=255)
                                caTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCaa                                                                                                                       <  2:1389659/150‑1 (MQ=255)
                                  tGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCTCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGc                                                                                                                     <  2:607865/150‑1 (MQ=255)
                                  tGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGc                                                                                                                     >  2:750235/1‑150 (MQ=255)
                                    gCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAg                                                                                                                   >  1:1406349/1‑150 (MQ=255)
                                     cGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                                  >  2:560432/1‑150 (MQ=255)
                                         aCAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCa                                                                                                              >  1:979749/1‑150 (MQ=255)
                                          cAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAg                                                                                                             <  1:1487043/150‑1 (MQ=255)
                                           aGCATCAGGGGATTATTCATGCGACTATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGc                                                                                                            <  1:494992/150‑1 (MQ=255)
                                            gCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACTAATTTTCATTATTCGCGATCATGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAgag                                                                                                           >  1:1489018/1‑148 (MQ=255)
                                            gCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGACACTCAAGCAGGCGCAGCg                                                                                                           >  2:814163/1‑150 (MQ=255)
                                             cATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGt                                                                                                          >  1:1091813/1‑150 (MQ=255)
                                                 aGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGt                                                                                                      <  2:359333/150‑1 (MQ=255)
                                                   gggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTAt                                                                                                    <  1:1532380/150‑1 (MQ=255)
                                                    ggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATg                                                                                                   <  1:660083/150‑1 (MQ=255)
                                                     gATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                                  >  1:1046868/1‑150 (MQ=255)
                                                      attattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATgag                                                                                                 <  2:162040/150‑1 (MQ=255)
                                                      attattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATgag                                                                                                 >  2:675816/1‑150 (MQ=255)
                                                        tattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaa                                                                                               <  1:675816/150‑1 (MQ=255)
                                                           tCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGAGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaaaa                                                                                            >  2:57520/1‑150 (MQ=255)
                                                             aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCATACAGACGCTGCGCCTGCTTgagtttcttg                                                                                          <  2:504991/150‑8 (MQ=255)
                                                              ggCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                         >  2:1555634/1‑150 (MQ=255)
                                                                   cgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAgg                                                                                    <  2:920368/150‑1 (MQ=255)
                                                                   cgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAgg                                                                                    <  1:1307039/150‑1 (MQ=255)
                                                                          tAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgat                                                                             <  1:923287/150‑1 (MQ=255)
                                                                              gAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATAGAGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                         <  1:1493234/150‑1 (MQ=255)
                                                                                  tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTTGGAGAAACACCAGCAGGGGCGGATGATCTTTa                                                                     >  1:74388/1‑150 (MQ=255)
                                                                                  tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                     <  2:1406349/150‑1 (MQ=255)
                                                                                   gTATGAGATCCGCATGTTCCGCCGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAACCCAGCAGGGGCGGATGATCTTTAt                                                                    >  1:650060/1‑150 (MQ=255)
                                                                                   gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCAATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGAGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTAt                                                                    <  2:1664684/150‑1 (MQ=255)
                                                                                     aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGa                                                                  <  2:713428/150‑1 (MQ=255)
                                                                                       gagcTCCGAATGTTCCGCAGCCATGTCGCGCCGCGCACACAAATATTCATTATTCGCTATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCATGCAGGCGCAGCGTCTGTTTGAGAAAAACCAGCTGGGGCTGATGATCTTTATGaaa                                                                <  2:1233464/150‑1 (MQ=255)
                                                                                         gaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGGCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGaaaaa                                                              <  2:345359/150‑1 (MQ=255)
                                                                                          aTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAAc                                                             <  1:1395599/150‑1 (MQ=255)
                                                                                           tCCTCATGTTCCGCAGCCAGGTCGCGCAGCGCACACAAATTTTCATTATTCGCGATCAGGATCGAGAAGTGCTGGTAAACATGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACa                                                            <  2:934459/150‑1 (MQ=255)
                                                                                           tCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACa                                                            >  2:1555226/1‑150 (MQ=255)
                                                                                           tCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACa                                                            <  2:928261/150‑1 (MQ=255)
                                                                                                  gTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGa                                                     >  1:129155/1‑150 (MQ=255)
                                                                                                     ccGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATgc                                                  <  1:1415950/150‑1 (MQ=255)
                                                                                                       gCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATgcgc                                                <  1:560432/150‑1 (MQ=255)
                                                                                                         aGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTg                                              >  1:1118296/1‑150 (MQ=255)
                                                                                                         aGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTg                                              <  1:478934/150‑1 (MQ=255)
                                                                                                         aGCCAGGTCGCGCCGCGCACACAAATTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGa                                             >  1:680604/1‑150 (MQ=255)
                                                                                                              ggTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCTAAGTGCTGGTAAACAAGAAACTCTAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGa                                         <  2:884760/150‑1 (MQ=255)
                                                                                                              ggTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGa                                         >  1:1456734/1‑150 (MQ=255)
                                                                                                                tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAAc                                       <  2:769134/150‑1 (MQ=255)
                                                                                                                 cgcgccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAAcc                                      >  1:1077408/1‑150 (MQ=255)
                                                                                                                  gcgccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAAccc                                     >  1:1410369/1‑150 (MQ=255)
                                                                                                                   cgccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCGGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                                    <  2:1718821/150‑1 (MQ=255)
                                                                                                                   cgccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                                    <  2:1475344/150‑1 (MQ=255)
                                                                                                                      cgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGcac                                 >  1:499126/1‑150 (MQ=255)
                                                                                                                      cgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGcac                                 <  1:1555226/150‑1 (MQ=255)
                                                                                                                          cacacaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTc                             >  2:1314001/1‑150 (MQ=255)
                                                                                                                              caAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTTTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg                         >  2:1003326/1‑150 (MQ=255)
                                                                                                                                    ccagcacttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCg                      >  1:190389/8‑147 (MQ=255)
                                                                                                                                    ccagcacttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCg                      <  2:190389/140‑1 (MQ=255)
                                                                                                                                      cATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGagag                 >  1:773509/1‑150 (MQ=255)
                                                                                                                                          cttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGcgcgcg             >  1:562001/2‑150 (MQ=255)
                                                                                                                                          cttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg             >  1:694196/2‑150 (MQ=255)
                                                                                                                                          cttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGACGGAGAGcgcg             >  2:1141369/2‑150 (MQ=255)
                                                                                                                                               cgATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTaa        <  1:770418/150‑1 (MQ=255)
                                                                                                                                                gATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTaaa       >  2:1168687/1‑150 (MQ=255)
                                                                                                                                                gATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTaaa       >  1:307400/1‑150 (MQ=255)
                                                                                                                                                 aTCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAAc      <  2:499126/150‑1 (MQ=255)
                                                                                                                                                  tCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACt     <  1:1292965/150‑1 (MQ=255)
                                                                                                                                                    aGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGTGGAGGATGATATTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCAAAAcaca   >  2:247652/1‑150 (MQ=255)
                                                                                                                                                    aGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTaagcaca   >  1:861950/1‑145 (MQ=255)
                                                                                                                                                    aGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCa   <  2:566054/150‑1 (MQ=255)
                                                                                                                                                     ggATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCaa  <  1:485259/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAA  >  NC_000913/2314652‑2314950

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: