Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,560,900 T→A 5.3% E231D (GAA→GAT)  PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,560,9000TA 5.3% 296.4 / 4.8 94E231D (GAA→GAT) PP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/45);  new base A (2/3);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACA  >  NC_002947/4560753‑4561048
                                                                                                                                                   |                                                                                                                                                    
gCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTc                                                                                                                                                    <  1:2114810/150‑1 (MQ=255)
gCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTc                                                                                                                                                    <  1:842752/150‑1 (MQ=255)
  cAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                  >  1:1691192/1‑150 (MQ=255)
  cAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                  <  1:905726/150‑1 (MQ=255)
  cAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                  >  1:985383/1‑150 (MQ=255)
  cAGAGGTCACCCGCCCAGCACTCATGGCCATGGGATCCGGATTTTGAGAATCCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                  >  1:1004017/1‑150 (MQ=255)
  cAGAGGTCACCCGCCCAGCACTAATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                  >  2:1656622/1‑150 (MQ=255)
   agagGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCg                                                                                                                                                 >  2:2264343/1‑150 (MQ=255)
      ggTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGCat                                                                                                                                              <  1:605292/150‑2 (MQ=255)
       gTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAg                                                                                                                                             <  1:1439688/150‑1 (MQ=255)
        tCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGCatgc                                                                                                                                            <  2:1384491/150‑4 (MQ=255)
          aCCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCa                                                                                                                                          <  1:1707836/150‑1 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                        >  1:794121/1‑150 (MQ=255)
            ccGACCAGCACTCATGGCAAGGGGCTCCGGCTGTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                        <  2:380936/150‑1 (MQ=255)
             cGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGc                                                                                                                                       <  1:2059862/150‑1 (MQ=255)
                   gCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGc                                                                                                                                 >  1:174585/1‑150 (MQ=255)
                    cACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCa                                                                                                                                >  1:579885/1‑150 (MQ=255)
                       tCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGACCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACgg                                                                                                                             >  1:1323098/1‑150 (MQ=255)
                             ccATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGa                                                                                                                       >  2:1100050/1‑150 (MQ=255)
                             ccATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGa                                                                                                                       <  1:1502682/150‑1 (MQ=255)
                              cATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGaa                                                                                                                      <  2:1998746/150‑1 (MQ=255)
                              cATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGaa                                                                                                                      <  1:2145514/150‑1 (MQ=255)
                                  ggCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgc                                                                                                                  <  2:794121/150‑1 (MQ=255)
                                  ggCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgc                                                                                                                  <  2:1141499/150‑1 (MQ=255)
                                     tCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCt                                                                                                               >  1:2286221/1‑150 (MQ=255)
                                      ccGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCtt                                                                                                              >  1:1885131/1‑150 (MQ=255)
                                       cGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTc                                                                                                             >  2:1353953/1‑150 (MQ=255)
                                       cGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTc                                                                                                             <  2:1691192/150‑1 (MQ=255)
                                               gcgcATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCt                                                                                                     >  1:455331/1‑150 (MQ=255)
                                                     gCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAAcc                                                                                               <  2:455331/150‑1 (MQ=255)
                                                     gCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAAcc                                                                                               <  1:495014/150‑1 (MQ=255)
                                                       aCCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGc                                                                                             >  2:397333/1‑150 (MQ=255)
                                                        ccccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCa                                                                                            <  1:1100050/150‑1 (MQ=255)
                                                          ccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAcc                                                                                          <  2:1383308/150‑1 (MQ=255)
                                                           cGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCa                                                                                         >  2:161630/1‑150 (MQ=255)
                                                           cGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCa                                                                                         >  1:286711/1‑150 (MQ=255)
                                                            ggCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCaa                                                                                        <  2:174585/150‑1 (MQ=255)
                                                                tGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAgg                                                                                    >  1:1534367/1‑150 (MQ=255)
                                                                     aaCCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGa                                                                               <  1:1935136/150‑1 (MQ=255)
                                                                      aCCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAg                                                                              >  1:1290593/1‑150 (MQ=255)
                                                                        cTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGgc                                                                            <  2:579885/150‑1 (MQ=255)
                                                                             cgGGCAGGTCGCTGGCCTGCTTCAGCACCTTTGCCAGCCTGGCATCGAGCTGCATCAGCGAATGCGGGATTTCCCGCAAGGCATGCATCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCt                                                                       <  1:1306739/150‑1 (MQ=255)
                                                                                    gTCGCTTGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGaa                                                                <  1:1698636/150‑1 (MQ=255)
                                                                                    gTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGaa                                                                <  1:2270866/150‑1 (MQ=255)
                                                                                    gTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGaa                                                                >  2:2043043/1‑150 (MQ=255)
                                                                                     tCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAAc                                                               <  2:124500/150‑1 (MQ=255)
                                                                                        cTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGt                                                                                                                          <  1:2386358/88‑1 (MQ=255)
                                                                                        cTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGt                                                                                                                          >  2:2386358/1‑88 (MQ=255)
                                                                                        cTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                            <  2:1290593/150‑1 (MQ=255)
                                                                                        cTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                            <  1:2104132/150‑1 (MQ=255)
                                                                                          ggCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCa                                                          <  2:1534367/150‑1 (MQ=255)
                                                                                            ccTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGt                                                        >  2:336917/1‑150 (MQ=255)
                                                                                            ccTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCACCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGt                                                        >  1:321138/1‑150 (MQ=255)
                                                                                             cTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgc                                                                                                                  <  1:1074724/91‑1 (MQ=255)
                                                                                             cTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgc                                                                                                                  >  2:1074724/1‑91 (MQ=255)
                                                                                             cTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGtt                                                       >  1:1375706/1‑150 (MQ=255)
                                                                                              tGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                            <  2:1641091/144‑1 (MQ=255)
                                                                                              tGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                            >  1:1641091/1‑144 (MQ=255)
                                                                                              tGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTg                                                      >  1:934308/1‑150 (MQ=255)
                                                                                                cTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGcc                                                    <  1:1656622/150‑1 (MQ=255)
                                                                                                 tCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGccc                                                   >  2:1044050/1‑150 (MQ=255)
                                                                                                 tCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGccc                                                   >  1:2312197/1‑150 (MQ=255)
                                                                                                  ccAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCa                                                  <  1:1073201/150‑1 (MQ=255)
                                                                                                  ccAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCa                                                  <  2:2312197/150‑1 (MQ=255)
                                                                                                   cAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAg                                                 >  2:2186023/1‑150 (MQ=255)
                                                                                                     gCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGt                                               <  2:537715/150‑1 (MQ=255)
                                                                                                         cTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGt                                               <  2:2073906/146‑1 (MQ=255)
                                                                                                         cTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGt                                               >  1:2073906/1‑146 (MQ=255)
                                                                                                            ggCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCg                                        >  1:1954546/1‑150 (MQ=255)
                                                                                                              ccAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCAt                                                                                                          <  2:2378028/82‑1 (MQ=255)
                                                                                                              ccAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGaa                                      <  1:2186023/150‑1 (MQ=255)
                                                                                                              ccAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTGCCAt                                                                                                          >  1:2378028/1‑82 (MQ=255)
                                                                                                               cAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAg                                     >  2:517596/1‑150 (MQ=255)
                                                                                                                  ccTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCtt                                  >  1:1952026/1‑150 (MQ=255)
                                                                                                                    tGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttg                                >  2:2276278/1‑150 (MQ=255)
                                                                                                                         tCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACAGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgtt                           <  2:1430263/150‑1 (MQ=255)
                                                                                                                         tCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACAGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgtt                           >  1:1430263/1‑150 (MQ=255)
                                                                                                                           gAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAG‑TGCCCAGGTAGCTGCGAAGCTttgccgttgcc                        >  1:2426836/1‑150 (MQ=255)
                                                                                                                               tGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCcgcg                     >  1:1103286/1‑150 (MQ=255)
                                                                                                                                gCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCga                    <  2:1885131/150‑1 (MQ=255)
                                                                                                                                gCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCga                    >  1:106816/1‑150 (MQ=255)
                                                                                                                                 cATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCgag                   >  1:1815831/1‑150 (MQ=255)
                                                                                                                                  aTCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCgaga                  <  2:1004017/150‑1 (MQ=255)
                                                                                                                                     atgccttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAAcc               >  1:697796/7‑150 (MQ=255)
                                                                                                                                     aGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAAcc               <  1:1573865/150‑1 (MQ=255)
                                                                                                                                        gcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg            <  2:2286221/150‑1 (MQ=255)
                                                                                                                                        gcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg            <  2:1103286/150‑1 (MQ=255)
                                                                                                                                        gcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg            >  1:1132397/1‑150 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg            >  2:1769237/3‑149 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg            <  1:1769237/147‑1 (MQ=255)
                                                                                                                                         cATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGTCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTgg           <  1:786847/150‑1 (MQ=255)
                                                                                                                                          aTGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGt          <  1:1044050/150‑1 (MQ=255)
                                                                                                                                             cGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCaa       >  1:817499/1‑150 (MQ=255)
                                                                                                                                              gggATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCaaa      >  2:46358/1‑150 (MQ=255)
                                                                                                                                                  tttCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACa  <  2:286711/150‑1 (MQ=255)
                                                                                                                                                  tttCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACa  >  1:2057574/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACA  >  NC_002947/4560753‑4561048

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