Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,314,801 G→T 5.6% 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,8010GT 5.6% 420.3 / 8.4 144Q438H (CAG→CATrcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base G (65/70);  new base T (5/4);  total (70/74)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

GCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAA  >  NC_000913/2314654‑2314950
                                                                                                                                                   |                                                                                                                                                     
gcCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGa                                                                                                                                                     >  2:639180/1‑150 (MQ=255)
   gCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcg                                                                                                                                                  <  2:1789508/150‑1 (MQ=255)
   gCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcg                                                                                                                                                  >  1:361927/1‑150 (MQ=255)
           aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaataat                                                                                                                                          <  1:626131/150‑5 (MQ=255)
            tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCt                                                                                                                                         <  1:1355706/150‑1 (MQ=255)
            tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCt                                                                                                                                         <  1:540050/150‑1 (MQ=255)
            tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCt                                                                                                                                         <  1:36043/150‑1 (MQ=255)
              tGCAAAACATCCCCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTAATCGCGACCCTGATCGCGaataatgaa                                                                                                                                       <  1:1138881/150‑8 (MQ=255)
               gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGt                                                                                                                                      <  1:1046361/150‑1 (MQ=255)
               gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGt                                                                                                                                      >  1:1079226/1‑150 (MQ=255)
                cAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTa                                                                                                                                     >  2:483127/1‑150 (MQ=255)
                    aCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGAAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCTAGATCAGGATCGCGAAGTGCTGCTAAACa                                                                                                                                 >  1:80758/1‑150 (MQ=255)
                     cATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                                <  2:414090/150‑1 (MQ=255)
                      aTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAg                                                                                                                               <  1:18557/150‑1 (MQ=255)
                       tCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGa                                                                                                                              <  1:980082/150‑1 (MQ=255)
                         accaccATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                            >  1:1392140/1‑150 (MQ=255)
                          ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCTCGATCAGGATCGCGAAGTTCTGGTAAACAAGCAAc                                                                                                                           >  2:159427/1‑150 (MQ=255)
                          ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAAc                                                                                                                           >  2:803068/1‑150 (MQ=255)
                          ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGACAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAAc                                                                                                                           <  1:111012/150‑1 (MQ=255)
                            accaTTGCGGAACAGCATCAGTGTATTATTCAGGCGACGATCAATAACGATATGTATGAGATCCGGATGTTCCGCAGCCAGGTCGCGCCTCGCTCACAAATTATCATTATTCGCGATCAGGATCGCGAAGTGCTGGTTAACAAGAAACTc                                                                                                                         <  1:307230/150‑1 (MQ=255)
                            accaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTc                                                                                                                         >  2:315701/1‑150 (MQ=255)
                              caTGGCGGTACACCATGAGGGGTTTATTCTGGCGTCGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCTAAGTGCTGGTAAACAAGAAtctcta                                                                                                                       <  2:335981/150‑3 (MQ=255)
                               aTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAg                                                                                                                      <  1:1479848/150‑1 (MQ=255)
                                tGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGc                                                                                                                     <  1:1628264/150‑1 (MQ=255)
                                tGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGc                                                                                                                     <  1:1287212/150‑1 (MQ=255)
                                   cGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                                  >  1:564038/1‑150 (MQ=255)
                                   cGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGATGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCCCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                                  >  1:404592/1‑150 (MQ=255)
                                            ttaaGGGGATTATTCAGGAGACGATCATTAACGAGCTGTCTGAGATCCGCATGTTCAGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCACGCAGGCGCAGCGTc                                                                                                         <  1:1150976/147‑1 (MQ=255)
                                            aTCAGGTGATTATTCAGGCGACGAGCAATAACGAGCTGTATGAGCTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTc                                                                                                         <  2:153106/150‑1 (MQ=255)
                                                  ggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATg                                                                                                   <  1:685717/150‑1 (MQ=255)
                                                      tattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaa                                                                                               <  1:655520/150‑1 (MQ=255)
                                                      tattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaa                                                                                               >  2:1113975/1‑150 (MQ=255)
                                                        ttCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCCGCCAGGTCGCGCCGCGCACACAAATTTTCATTTTTCGCGAGCAGGATGGGGAAGTGCGGGTAAGCAAGAACCGCAAGCAGGCGCAGCGTCTGTAGGAGaaaa                                                                                             >  2:278576/1‑150 (MQ=255)
                                                         tCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaaaa                                                                                            >  2:749129/1‑150 (MQ=255)
                                                          cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCTCAGCCAGGTCGCGCCGAGCACACAAATTTTCATTATTCGCGATCAGGATCTCGAAGTGCTGGTAATCAAGAAACTCAATCAGGCGCAGCGGCTGTATGAGaaaaaa                                                                                           >  1:1748413/1‑149 (MQ=255)
                                                          cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                           >  2:1471082/1‑150 (MQ=255)
                                                          cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                           >  1:507736/1‑150 (MQ=255)
                                                          cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAAATTTCATTATTCGCGATCAGGATCGCGAAGTGCTGCTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAc                                                                                           >  2:1076145/1‑150 (MQ=255)
                                                           aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAcc                                                                                          <  2:1534449/150‑1 (MQ=255)
                                                           aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAcc                                                                                          >  2:1061881/1‑150 (MQ=255)
                                                            ggCGACTATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                         <  1:1740853/150‑1 (MQ=255)
                                                                cctaTCAATAACTCGCTGATTGTCACCCGCAGGTGCCGCAGCCAGGCCGCGCCGCGCACAAAAATTTTCATAATTCGCCATCAGCAGCCCGAAGTGCTGGTAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAgg                                                                                    <  2:1400588/147‑1 (MQ=255)
                                                                 cgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAgg                                                                                    >  1:144170/1‑150 (MQ=255)
                                                                      aaTAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCgg                                                                               >  2:911180/1‑150 (MQ=255)
                                                                      aaTAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGAGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCgg                                                                               >  1:983823/1‑150 (MQ=255)
                                                                        tAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgat                                                                             <  2:534076/150‑1 (MQ=255)
                                                                           cGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatgat                                                                          <  1:918353/150‑1 (MQ=255)
                                                                           cGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatgat                                                                          <  1:1465225/150‑1 (MQ=255)
                                                                               cTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCttt                                                                      >  2:933420/1‑150 (MQ=255)
                                                                                 gTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTAt                                                                    <  1:494401/150‑1 (MQ=255)
                                                                                  tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATg                                                                   <  1:1496084/150‑1 (MQ=255)
                                                                                    tGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGCGGCGGATGATCTTTATGaa                                                                 >  2:1104723/1‑150 (MQ=255)
                                                                                     gagaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGaaa                                                                >  1:154940/1‑150 (MQ=255)
                                                                                       gaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGaaaaa                                                              >  2:616147/1‑150 (MQ=255)
                                                                                          ccGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACAt                                                           >  2:1215562/1‑150 (MQ=255)
                                                                                            gCATGTTCCGCAGCCAGGTCGCGCGGCGCACACAATTTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTg                                                         <  2:431843/150‑1 (MQ=255)
                                                                                             cATGTTCCGCAGCCAGGACGCGCCGCGCACACAAATTTTCATTATTCGCGAACAGGATAGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGTCGCAGCGTCTGTATGAGAAAAACCAGCAGCGGCGGATGATCTTTATGAAAAACATTgg                                                        >  1:624496/1‑150 (MQ=255)
                                                                                              aTGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGc                                                       <  1:1808774/150‑1 (MQ=255)
                                                                                               ttatCCGCAGCCAGGTAGCGCCGCGCACAAAAATTTTCATTATACGCGATAAGGATCGCGAAGAGCTGGTAAACAAGAAACTAAAGCAGGCGCAGCCTCTGAATGAGAAAAACCAGCAGGGGCAGAAGAGATTTATGAAAAACATTGGCg                                                      >  2:1673721/4‑150 (MQ=255)
                                                                                               tGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCg                                                      >  1:862666/1‑150 (MQ=255)
                                                                                               tGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCg                                                      >  2:489246/1‑150 (MQ=255)
                                                                                                gTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGa                                                     >  1:1515311/1‑150 (MQ=255)
                                                                                                gTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGa                                                     >  2:1000456/1‑150 (MQ=255)
                                                                                                  tCCGCAGCCAGGTAGCGCCGCGCACACAAATTTTCATTATTCGCGATCATGATCGCGAAGTGCTGGTAAACAAGAAACTCAATCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATg                                                   <  1:1042/150‑1 (MQ=255)
                                                                                                   ccGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGTTCGCGATGTGCTGGTAAGCAAGCAACTCAAGCAGGCGCAGCGTCTGTATGAGAATAACCAGCAGGGGCGGATGAGCTTTATGAAAAACATTGGCGATgc                                                  >  2:1306211/1‑150 (MQ=255)
                                                                                                   ccGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGATGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATgc                                                  >  2:1069054/1‑150 (MQ=255)
                                                                                                       cgCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTg                                              <  2:1687465/149‑1 (MQ=255)
                                                                                                       aGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTg                                              <  1:1602479/150‑1 (MQ=255)
                                                                                                          cAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGaaa                                           >  1:800557/1‑150 (MQ=255)
                                                                                                           aGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAATAACATTGGCGATTCGCTGaaat                                          >  1:648122/1‑149 (MQ=255)
                                                                                                           aGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAg                                          <  1:1000456/150‑1 (MQ=255)
                                                                                                            ggTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAACCATTGGCGATGCGCTGAAAGa                                         >  1:1668180/1‑150 (MQ=255)
                                                                                                               cgcgccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAAcc                                      <  2:571054/150‑1 (MQ=255)
                                                                                                                 ccccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGCAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                                    <  1:1201953/148‑1 (MQ=255)
                                                                                                                  gccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGc                                   <  2:1804411/150‑1 (MQ=255)
                                                                                                                  gccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGc                                   >  2:1313227/1‑150 (MQ=255)
                                                                                                                  gccgcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAAAAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGc                                   <  2:361927/150‑1 (MQ=255)
                                                                                                                     gcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGcaca                                <  1:913873/150‑1 (MQ=255)
                                                                                                                      cgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTTGTAAACAAGAAACTCAAGCAGGCGCAGCGGCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAg                               <  2:1400989/150‑1 (MQ=255)
                                                                                                                      cgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAg                               <  2:676173/150‑1 (MQ=255)
                                                                                                                         acacaTATTATCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTcc                            <  1:336567/150‑1 (MQ=255)
                                                                                                                         acacaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGACACTCAAGCAGGCGCAGCGTCTGTATGCGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGCTGCGCTGAAAGAACCCGCACAGTcc                            >  1:1385742/1‑150 (MQ=255)
                                                                                                                          cacaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                           >  1:735294/1‑150 (MQ=255)
                                                                                                                          cacaAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTccc                           >  1:996922/1‑150 (MQ=255)
                                                                                                                           acaAATTTTCATTATTCGCGAGCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCt                          <  1:1654839/150‑1 (MQ=255)
                                                                                                                            caAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg                         <  1:1259490/150‑1 (MQ=255)
                                                                                                                                    cATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGagag                 <  1:1724981/150‑1 (MQ=255)
                                                                                                                                     attattCGCGATCAGGATCGCGAAGTGCTGGTAACCAAGAAACTCACGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCCCAGTCCCTGGCGGAGAgc                >  2:1462877/1‑150 (MQ=255)
                                                                                                                                      ttattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAgcg               >  1:204235/1‑150 (MQ=255)
                                                                                                                                      cacttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAACACAGCAGGGGCGGAAGATCATTATGAAAAACATTGTCGATTCTCTGAAAGAACCCGCAAAGTCCCTGGCGGAGAgcg               >  1:1220393/4‑150 (MQ=255)
                                                                                                                                       tattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAgcgc              >  1:555326/1‑150 (MQ=255)
                                                                                                                                        attCGCGATCAGGATCGCGAAGTGCTGGTAACCAAGAAACTCAAGCCCGCGCAGCGGCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg             <  1:450385/150‑1 (MQ=255)
                                                                                                                                        attCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg             <  2:618574/150‑1 (MQ=255)
                                                                                                                                        attCGCGATCAGGATCGCGAAGTGCTCGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCTGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg             <  2:1177864/150‑1 (MQ=255)
                                                                                                                                         ttCGCGATCCTGCTCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGCTGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCgg            >  1:563/1‑150 (MQ=255)
                                                                                                                                              gATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTaaa       >  2:293316/1‑150 (MQ=255)
                                                                                                                                                tCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACt     <  1:1471082/150‑1 (MQ=255)
                                                                                                                                                   ggATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCaa  >  2:1117360/1‑150 (MQ=255)
                                                                                                                                                   ggATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCaa  >  2:1350079/1‑150 (MQ=255)
                                                                                                                                                   ggATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCaa  >  2:1597968/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACTCAA  >  NC_000913/2314654‑2314950

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