Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,203,911 G→T 5.3% A1303A (GCC→GCA PP_3691 ← DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,203,9110GT 5.3% 237.4 / 2.7 75A1303A (GCC→GCAPP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (36/35);  new base T (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

TGCAACAACGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAACAGGTCGAA  >  NC_002947/4203762‑4204060
                                                                                                                                                     |                                                                                                                                                     
tGCAACAACGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACg                                                                                                                                                       >  2:1665313/1‑150 (MQ=255)
tGCAACAACGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACg                                                                                                                                                       <  1:1687903/150‑1 (MQ=255)
   aacaacGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCa                                                                                                                                                    <  1:11681/150‑1 (MQ=255)
   aacaacGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCa                                                                                                                                                    <  1:11374/150‑1 (MQ=255)
        cGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGc                                                                                                                                               >  2:1274474/1‑150 (MQ=255)
          cAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTa                                                                                                                                             >  1:1158365/1‑150 (MQ=255)
                tttCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGcc                                                                                                                                       >  1:399374/1‑150 (MQ=255)
                     cGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTgg                                                                                                                                  >  1:1879607/1‑150 (MQ=255)
                      gCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGa                                                                                                                                 >  1:1861486/1‑150 (MQ=255)
                           ccAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGcc                                                                                                                            <  2:1879607/150‑1 (MQ=255)
                            cAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCca                                                                                                                           >  1:1166924/1‑150 (MQ=255)
                              aCGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCcaca                                                                                                                         >  2:771057/1‑150 (MQ=255)
                                gCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGc                                                                                                                       <  2:777790/150‑1 (MQ=255)
                                  cTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCt                                                                                                                     <  1:556207/150‑1 (MQ=255)
                                     cGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGAc                                                                                                                  <  1:174591/150‑1 (MQ=255)
                                       agagCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGt                                                                                                                <  2:1861486/150‑1 (MQ=255)
                                           cTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGAGCACTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCaa                                                                                                            <  1:1107538/150‑1 (MQ=255)
                                            tATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAAt                                                                                                           <  2:1006526/150‑1 (MQ=255)
                                             aTCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAAtc                                                                                                          <  1:771057/150‑1 (MQ=255)
                                             aTCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAAtc                                                                                                          >  1:1414289/1‑150 (MQ=255)
                                              tCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAAtct                                                                                                         >  2:1591771/1‑150 (MQ=255)
                                              tCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGGTTGCAGGGGAGCCATGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCAACAGCCTGACGTGCAAtct                                                                                                         >  1:1018712/1‑150 (MQ=255)
                                                 cACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCaa                                                                                                      <  2:399374/150‑1 (MQ=255)
                                                      gCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTg                                                                                                 <  1:616564/150‑1 (MQ=255)
                                                          gCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCaa                                                                                                      <  2:394933/141‑1 (MQ=255)
                                                          gCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCaa                                                                                                      >  1:394933/1‑141 (MQ=255)
                                                            tctAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTc                                                                                           <  2:1158365/150‑1 (MQ=255)
                                                               aGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAgg                                                                                        >  1:1414320/1‑150 (MQ=255)
                                                               aGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAgg                                                                                        <  2:1173132/150‑1 (MQ=255)
                                                                gCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTTCCCTGTGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGt                                                                                       <  1:338922/150‑1 (MQ=255)
                                                                gCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGt                                                                                       >  1:565178/1‑150 (MQ=255)
                                                                 cAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCtt                                                                                                                                    <  1:195141/104‑1 (MQ=255)
                                                                 cAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCtt                                                                                                                                    >  2:195141/1‑104 (MQ=255)
                                                                      cGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTcc                                                                                 >  2:1144365/1‑150 (MQ=255)
                                                                          gCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGAc                                                                             <  2:565178/150‑1 (MQ=255)
                                                                            cATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCg                                                                           >  1:522480/1‑150 (MQ=255)
                                                                               cGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTa                                                                        <  1:1274474/150‑1 (MQ=255)
                                                                                 cccAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTAcc                                                                      >  2:1498371/1‑150 (MQ=255)
                                                                                  ccAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCa                                                                     >  2:1300482/1‑150 (MQ=255)
                                                                                  ccAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCa                                                                     >  2:505251/1‑150 (MQ=255)
                                                                                       agCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAActct                                                                <  2:1029298/150‑1 (MQ=255)
                                                                                       agCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAActct                                                                <  1:1300482/150‑1 (MQ=255)
                                                                                          gTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGACGTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGc                                                             >  2:1911255/1‑150 (MQ=255)
                                                                                           tttGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACTGCCGTTGcg                                                                                                                                              <  2:1906112/68‑2 (MQ=255)
                                                                                           tttGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACTGCCGTTGcg                                                                                                                                              >  1:1906112/1‑67 (MQ=255)
                                                                                           tttGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCa                                                            <  1:1665313/150‑1 (MQ=255)
                                                                                              gAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCg                                                                           <  2:1747871/132‑1 (MQ=255)
                                                                                              gAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCg                                                                           >  1:1747871/1‑132 (MQ=255)
                                                                                                aGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCa                                                       <  1:1498371/150‑1 (MQ=255)
                                                                                                 ggggAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAg                                                      >  1:1057139/1‑150 (MQ=255)
                                                                                                  gggAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGc                                                     >  1:32754/1‑150 (MQ=255)
                                                                                                    gAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACTGCCGTTGcgaa                                                                                                                                            <  2:988364/61‑4 (MQ=255)
                                                                                                    gAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACTGCCGTTGcgaa                                                                                                                                            >  1:988364/1‑58 (MQ=255)
                                                                                                    gAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCcaca                                                                                                                         >  1:1033770/1‑80 (MQ=255)
                                                                                                    gAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCcaca                                                                                                                         <  2:1033770/80‑1 (MQ=255)
                                                                                                         ttGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGa                                                                                                                   <  2:923526/81‑1 (MQ=255)
                                                                                                         ttGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGa                                                                                                                   >  1:923526/1‑81 (MQ=255)
                                                                                                               gCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGa                                        >  2:692392/1‑150 (MQ=255)
                                                                                                                 tGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAg                                                                                         >  2:1418482/1‑99 (MQ=255)
                                                                                                                 tGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAg                                                                                         <  1:1418482/99‑1 (MQ=255)
                                                                                                                           tCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAAcc                            >  1:1170674/1‑150 (MQ=255)
                                                                                                                                     gcgcGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAAccgctgcccgct                  >  2:1931169/1‑150 (MQ=255)
                                                                                                                                     gcgcGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAAccg                           >  2:1078147/1‑141 (MQ=255)
                                                                                                                                     gcgcGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAAccg                           <  1:1078147/141‑1 (MQ=255)
                                                                                                                                         gAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGCCGGACAGATCACGTTCAACCGCTGCCCGCTGCtt              <  1:1591771/150‑1 (MQ=255)
                                                                                                                                            cTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTAc                                                                       >  2:20749/1‑90 (MQ=255)
                                                                                                                                            cTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTaa           >  2:1356081/1‑150 (MQ=255)
                                                                                                                                            cTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTaa           <  2:1170674/150‑1 (MQ=255)
                                                                                                                                            cTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCAACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTaa           <  2:1018712/150‑1 (MQ=255)
                                                                                                                                            cTTCGCAACGGCAGTTGCTACCTGCCCGTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTAc                                                                       <  1:20749/90‑1 (MQ=255)
                                                                                                                                             ttCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAAc          >  1:146332/1‑150 (MQ=255)
                                                                                                                                               cGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAACAg        <  2:522480/150‑1 (MQ=255)
                                                                                                                                                gCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCtt              <  1:1288105/143‑1 (MQ=255)
                                                                                                                                                gCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCtt              >  2:1288105/1‑143 (MQ=255)
                                                                                                                                                gCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAACAgg       >  2:1261952/1‑150 (MQ=255)
                                                                                                                                                 cAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAACAGGt      <  1:1911255/150‑1 (MQ=255)
                                                                                                                                                  aaCGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTATGACAGGTCGGACAGATCACGTTCAACAGCTGCCCGCTGCTTTAAGAGGTc     >  1:1087602/1‑150 (MQ=255)
                                                                                                                                                     ggCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAACAGGTCGaa  <  1:1144365/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAACAACGCAAGCCTTTCACGCAGGCCAACGCCTGCGAGAGCTATCGCACCTGCAGGCTCTAGCAAAGCGTTGCCATCGCCCAAGAGCGTTTGAAGGGGAGCCTTGTAAGCTGCAATCTGCTCCGGGTTCAGCGCGAGCTTCGCAACGGCAGTTGCTACCTGCCCTTGGAACGCCACAGCCTGACGTGCAATCTCAAGCTTGGTCTTCAGGTCACTCCAGACCGGTACCAACTCTCGCAGGTCAGCCAGGTCGGACAGATCACGTTCAACCGCTGCCCGCTGCTTTAACAGGTCGAA  >  NC_002947/4203762‑4204060

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