Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,446,380 A→C 10.3% *159G (TGA→GGA)  sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,446,3800AC10.3% 490.7 / 35.0 162*159G (TGA→GGA) sulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base A (69/76);  new base C (10/7);  total (79/83)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

AGTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGC  >  NC_002947/2446231‑2446529
                                                                                                                                                     |                                                                                                                                                     
aGTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCa                                                                                                                                                       <  2:108808/150‑1 (MQ=255)
aGTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCa                                                                                                                                                       <  1:3064344/150‑1 (MQ=255)
    gCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCg                                                                                                                                                   >  2:4369043/1‑150 (MQ=255)
      cAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAg                                                                                                                                                 >  2:2722973/1‑150 (MQ=255)
      cAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCGCTTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAg                                                                                                                                                 >  2:1564317/1‑150 (MQ=255)
        aGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGc                                                                                                                                               <  2:2490369/150‑1 (MQ=255)
         ggCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCg                                                                                                                                              >  2:950228/1‑150 (MQ=255)
         ggCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCg                                                                                                                                              <  1:4709517/150‑1 (MQ=255)
         ggCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCg                                                                                                                                              <  1:2016363/150‑1 (MQ=255)
           ccATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGa                                                                                                                                            <  1:1265899/150‑1 (MQ=255)
           ccATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGa                                                                                                                                            <  1:1756028/150‑1 (MQ=255)
            cATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGAt                                                                                                                                           <  1:350249/150‑1 (MQ=255)
                cTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGGACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTc                                                                                                                                       <  2:4376765/150‑1 (MQ=255)
                cTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTc                                                                                                                                       <  1:984256/150‑1 (MQ=255)
                 tCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                      <  1:4116801/150‑1 (MQ=255)
                 tCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCa                                                                                                                                      <  2:1577395/150‑1 (MQ=255)
                 tCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGACAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                      <  2:2070455/150‑1 (MQ=255)
                     tgtgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCTGCAGCGTATCAGGCTATTTGATGTGACAGGCAATGCCTGAACATCAGCCGAGGCGGATGTTCAGGCt                                                                                                                                  >  2:4426659/1‑150 (MQ=255)
                     tgtgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCt                                                                                                                                  >  1:4415494/1‑150 (MQ=255)
                      gtgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCtt                                                                                                                                 <  2:592203/150‑1 (MQ=255)
                       tgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCttt                                                                                                                                >  1:574744/1‑150 (MQ=255)
                         ggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTtgt                                                                                                                              <  1:2400933/150‑1 (MQ=255)
                                aGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTc                                                                                                                       <  2:891059/150‑1 (MQ=255)
                                aGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTc                                                                                                                       >  1:50546/1‑150 (MQ=255)
                                aGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTGTGTGCATTTc                                                                                                                       >  2:745923/1‑150 (MQ=255)
                                 gATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTcc                                                                                                                      <  2:2130687/150‑1 (MQ=255)
                                 gATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTcc                                                                                                                      >  2:350290/1‑150 (MQ=255)
                                 gATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTcc                                                                                                                      <  1:1049625/150‑1 (MQ=255)
                                    ggCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGc                                                                                                                   >  2:4076934/1‑150 (MQ=255)
                                       cGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACt                                                                                                                >  1:3561993/1‑150 (MQ=255)
                                          tCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGc                                                                                                             <  1:1045766/150‑1 (MQ=255)
                                            gCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTg                                                                                                           >  2:2830713/1‑150 (MQ=255)
                                             ccTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGc                                                                                                          <  2:1335609/150‑1 (MQ=255)
                                               tGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCAc                                                                                                        >  2:1933556/1‑150 (MQ=255)
                                                gatggatgGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACg                                                                                                       >  2:1414369/1‑150 (MQ=255)
                                                  tggatggCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTa                                                                                                     <  1:94376/150‑1 (MQ=255)
                                                     atggCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCa                                                                                                  <  1:2031244/150‑1 (MQ=255)
                                                     atggCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCa                                                                                                  <  2:3014026/150‑1 (MQ=255)
                                                            gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGttgttg                                                                                           <  1:809579/150‑1 (MQ=255)
                                                            gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCATTGTGCATTTCCGGCACTGGCTGCACGTAGCAGttgttg                                                                                           <  2:1811347/150‑1 (MQ=255)
                                                              tCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGc                                                                                         >  1:2652760/1‑150 (MQ=255)
                                                                      gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCg                                                                                 >  1:83336/1‑150 (MQ=255)
                                                                      gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCg                                                                                 <  2:3261924/150‑1 (MQ=255)
                                                                      gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCg                                                                                 <  2:696633/150‑1 (MQ=255)
                                                                      gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTACCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCg                                                                                 >  2:783135/1‑150 (MQ=255)
                                                                        tgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGc                                                                               >  2:619552/1‑150 (MQ=255)
                                                                         gaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCa                                                                              >  2:2269789/1‑150 (MQ=255)
                                                                          aaGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCAt                                                                             <  1:4880600/150‑1 (MQ=255)
                                                                           aGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCAtt                                                                            >  2:600188/1‑150 (MQ=255)
                                                                           aGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTGTGTGCATTTCCGGCACTGGCTGCACGTAGCGGGTGTTGGCGTGCAGCGGCAtt                                                                            >  2:2798101/1‑150 (MQ=255)
                                                                            gggCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAATGCCTGAACATCAGCCGAGGCGGATGTTCAGGCGGTGTGCATTTCCGGCACTGGCTGCACGGAGCAGGGGTTGGCGTGCAGCGGCAGTg                                                                           >  2:4660786/1‑150 (MQ=255)
                                                                              gCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCTCGTAGCAGTTGTTGGCGTGCAGCGGCATTGAc                                                                         <  2:2351706/150‑1 (MQ=255)
                                                                                     gtgCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGccc                                                                  >  1:4384345/1‑150 (MQ=255)
                                                                                      tgctgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCcca                                                                 <  2:4886540/150‑1 (MQ=255)
                                                                                        ctgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCccagc                                                               <  1:719658/150‑1 (MQ=255)
                                                                                         tgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCCGGCTTTGTGCCTGTCCGGCACTGGCTGCACGTAGCAGTTGGTGGCGTGCAGCGGCAGTGACGTTGCccagcc                                                              >  2:3157449/1‑150 (MQ=255)
                                                                                          gcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCccagcca                                                             >  1:4709628/1‑150 (MQ=255)
                                                                                                     gcagcaGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACgg                                                  <  1:3785617/150‑1 (MQ=255)
                                                                                                      cagcagCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGgt                                                 <  1:391417/150‑1 (MQ=255)
                                                                                                       agcagcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGgtg                                                <  2:405274/150‑1 (MQ=255)
                                                                                                         cagcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGtgtg                                              >  2:3817165/1‑150 (MQ=255)
                                                                                                         cagcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCGTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGtgtg                                              >  2:2818798/1‑150 (MQ=255)
                                                                                                          agcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTgg                                             <  2:478591/150‑1 (MQ=255)
                                                                                                          agcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTgg                                             <  2:1326894/150‑1 (MQ=255)
                                                                                                              tATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGa                                         <  2:1930619/150‑1 (MQ=255)
                                                                                                              tATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGa                                         <  2:2481215/150‑1 (MQ=255)
                                                                                                               aTCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAg                                        >  1:2472731/1‑150 (MQ=255)
                                                                                                                tCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAgg                                       >  1:3118263/1‑150 (MQ=255)
                                                                                                                  aGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTCGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                     >  2:1444302/1‑150 (MQ=255)
                                                                                                                    gCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCa                                   >  2:3525249/1‑150 (MQ=255)
                                                                                                                     cTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAg                                  <  1:624585/150‑1 (MQ=255)
                                                                                                                     cTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCACTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAg                                  <  1:275092/150‑1 (MQ=255)
                                                                                                                       aTTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGt                                <  2:899564/150‑1 (MQ=255)
                                                                                                                       aTTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGt                                <  2:899348/150‑1 (MQ=255)
                                                                                                                           gAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaa                            <  1:4369043/150‑1 (MQ=255)
                                                                                                                                 atgcacAACGCCTGAACATCCGCCTCGGCTGTTGTTCTGGATTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTc                      <  1:1900869/145‑1 (MQ=255)
                                                                                                                                 aCAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTc                      <  1:1705059/150‑1 (MQ=255)
                                                                                                                                   aGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTcac                    >  2:3419506/1‑150 (MQ=255)
                                                                                                                                   aGGCAAAGCCTGAACATCAGCCAAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCTGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTcac                    >  1:505251/1‑150 (MQ=255)
                                                                                                                                   aGGCAAAGCCTGAACATCAGCCAAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTcac                    >  2:3001593/1‑150 (MQ=255)
                                                                                                                                    ggCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGTCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTcaca                   <  1:3825294/150‑1 (MQ=255)
                                                                                                                                    ggCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTcaca                   <  1:2269789/150‑1 (MQ=255)
                                                                                                                                     gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                  >  2:5022254/1‑150 (MQ=255)
                                                                                                                                     gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                  >  1:305771/1‑150 (MQ=255)
                                                                                                                                     gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                  >  2:1612659/1‑150 (MQ=255)
                                                                                                                                     gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                  >  2:3982124/1‑150 (MQ=255)
                                                                                                                                     acAAAGCCTGAACATCCGCCTCGGCTGATGTTCAGGCGTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                  >  2:4326987/2‑150 (MQ=39)
                                                                                                                                         aGCCTGAACATCAGCCGAGGCGGATGTTCAGGCATTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCa              <  1:1817974/150‑1 (MQ=255)
                                                                                                                                          gCCTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGCGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAg             >  1:2545744/1‑150 (MQ=39)
                                                                                                                                           ccTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGt            >  1:3949586/1‑150 (MQ=255)
                                                                                                                                           ccTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGt            >  1:3949823/1‑150 (MQ=255)
                                                                                                                                           ccTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGt            >  2:3526513/1‑150 (MQ=255)
                                                                                                                                            cTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGGTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGtt           >  1:4064918/1‑150 (MQ=255)
                                                                                                                                             tGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTg          >  1:2047644/1‑150 (MQ=255)
                                                                                                                                               aaCATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCa        >  2:4711510/1‑150 (MQ=255)
                                                                                                                                               aaCATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCa        <  1:68005/150‑1 (MQ=255)
                                                                                                                                                 cATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgc      <  2:1132277/150‑1 (MQ=255)
                                                                                                                                                  aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGtt           <  2:4127181/144‑1 (MQ=255)
                                                                                                                                                  aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgcg     >  1:2612885/1‑150 (MQ=255)
                                                                                                                                                     aGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgcgcgc  >  2:3906389/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGC  >  NC_002947/2446231‑2446529

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