Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,681,341 C→A 8.5% N123K (AAC→AAA PP_1476 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,681,3410CA 8.5% 211.5 / 6.3 71N123K (AAC→AAAPP_1476hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base A (3/3);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CTGGACGATCGCGTGCTGCACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCA  >  NC_002947/1681192‑1681490
                                                                                                                                                      |                                                                                                                                                     
cTGGACGATCGCGTGCTGCACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAAc                                                                                                                                                       >  1:742691/1‑150 (MQ=255)
 tGGACGATCGCGTGCTGCACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAAcc                                                                                                                                                      >  2:192888/1‑150 (MQ=255)
  ggACGATCGCGTGCTGCTCCGCAGTTAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGCGCTCGCGCCGTAACGACCGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCt                                                                                                                                                     <  1:1533347/150‑1 (MQ=255)
    aCGATCGCGTGCTGCACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGa                                                                                                                                                   <  1:456974/150‑1 (MQ=255)
             tgctgcACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGa                                                                                                                                          <  2:1124444/150‑1 (MQ=255)
                  cACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGACGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCAtt                                                                                                                                     <  2:1037606/150‑1 (MQ=255)
                   aCTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTg                                                                                                                                    >  1:1108746/1‑150 (MQ=255)
                     tGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGcc                                                                                                                                  <  1:573103/150‑1 (MQ=255)
                          tGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCtgt                                                                                                                             <  2:1557829/150‑1 (MQ=255)
                                 ggCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCaa                                                                                                                      <  1:192888/150‑1 (MQ=255)
                                   ccAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAaga                                                                                                                    >  1:1737292/1‑150 (MQ=255)
                                   ccAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAaga                                                                                                                    >  1:342638/1‑150 (MQ=255)
                                   ccAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAaga                                                                                                                    >  2:831318/1‑150 (MQ=255)
                                    cAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACAGTTATGGTCAAGGTGTTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAagag                                                                                                                   <  2:1483463/150‑1 (MQ=255)
                                     aGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAgg                                                                                                                  >  1:1229508/1‑150 (MQ=255)
                                          gATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTCCGACTGTTGTTGGGTAAAGTTTTGATCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTg                                                                                                             <  1:231640/150‑1 (MQ=255)
                                           aTGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTgg                                                                                                            >  1:108985/1‑150 (MQ=255)
                                                 aTCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAAccc                                                                                                      >  1:1379595/1‑150 (MQ=255)
                                                   cGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTg                                                                                                    <  2:460247/150‑1 (MQ=255)
                                                    gCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTgg                                                                                                   <  1:1052545/150‑1 (MQ=255)
                                                      cTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcg                                                                                                 >  1:1859051/1‑150 (MQ=255)
                                                        gACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcgac                                                                                               >  2:1615580/1‑150 (MQ=255)
                                                         aCCTGGTGCTCGCGCCGTACCGAATGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcgacc                                                                                              >  2:1282906/1‑150 (MQ=255)
                                                           cTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcgaccga                                                                                            >  2:1673951/1‑150 (MQ=255)
                                                           cTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcgaccga                                                                                            >  2:1477615/1‑150 (MQ=255)
                                                              gTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACAt                                                                                         <  2:1713603/150‑1 (MQ=255)
                                                              gTGCTCGCGCCGGACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACAt                                                                                         <  1:1529260/150‑1 (MQ=255)
                                                                 cTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGcc                                                                                      <  1:1477615/150‑1 (MQ=255)
                                                                  tCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATg                                                                                        >  2:1724244/1‑147 (MQ=255)
                                                                  tCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATg                                                                                        <  1:1724244/147‑1 (MQ=255)
                                                                   cgcgCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCcgc                                                                                    <  1:1777794/150‑1 (MQ=255)
                                                                   cgcgCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGAGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGTCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCcgc                                                                                    <  2:573946/150‑1 (MQ=255)
                                                                            ccGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCa                                                                           >  1:1884721/1‑150 (MQ=255)
                                                                                 aCGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTg                                                                       <  1:622273/150‑1 (MQ=255)
                                                                                  cGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTgg                                                                      <  1:1673951/150‑1 (MQ=255)
                                                                                   gCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTggg                                                                     >  1:1750111/1‑150 (MQ=255)
                                                                                     aCCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGAc                                                                   >  2:341423/1‑150 (MQ=255)
                                                                                         gCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtga                                                               >  2:1300486/1‑150 (MQ=255)
                                                                                              cTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATc                                                          >  1:493682/1‑150 (MQ=255)
                                                                                                ggCTTACGGCAGCCAGCGAGTGTGGGAGGGACAGCAGGTCAGGGCTGTAGCACAGCTCAAAAACCGAAGCAGTTCCGCGGTGCCTCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAg                                                        <  2:1859051/150‑1 (MQ=255)
                                                                                                       ggCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGt                                                 >  2:1666561/1‑150 (MQ=255)
                                                                                                          aGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACg                                              <  2:436195/150‑1 (MQ=255)
                                                                                                            ccTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGt                                            >  2:223737/1‑150 (MQ=255)
                                                                                                              tGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCt                                                                                                     <  1:1792128/91‑1 (MQ=255)
                                                                                                              tGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCt                                                                                                     >  2:1792128/1‑91 (MQ=255)
                                                                                                                  aGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGt                                      <  1:1702089/150‑1 (MQ=255)
                                                                                                                   gtgtGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGtg                                     >  2:474695/1‑150 (MQ=255)
                                                                                                                        ggAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcgaccga                                                                                            >  2:1839548/1‑90 (MQ=255)
                                                                                                                        ggAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGcgaccga                                                                                            <  1:1839548/90‑1 (MQ=255)
                                                                                                                             gTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACg                                              >  1:1549463/1‑131 (MQ=255)
                                                                                                                             gTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACg                                              <  2:1549463/131‑1 (MQ=255)
                                                                                                                              tAAAGTGGTGAGCGCGGTCAAGAAACTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtg                                                                <  2:646254/112‑1 (MQ=255)
                                                                                                                              tAAAGTGGTGAGCGCGGTCAAGAAACTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtg                                                                <  2:647013/112‑1 (MQ=255)
                                                                                                                              tAAAGTGGTGAGCGCGGTCAAGAAACTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtg                                                                <  2:653012/112‑1 (MQ=255)
                                                                                                                              tAAAGTGGTGAGCGCGGTCAAGAAACTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtg                                                                >  1:653012/1‑112 (MQ=255)
                                                                                                                              tAAAGTGGTGAGCGCGGTCAAGAAACTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtg                                                                >  1:647013/1‑112 (MQ=255)
                                                                                                                              tAAAGTGGTGAGCGCGGTCAAGAAACTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGtg                                                                >  1:646254/1‑112 (MQ=255)
                                                                                                                                aaGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTg                        >  1:875874/1‑150 (MQ=255)
                                                                                                                                 aGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTgg                       <  2:108985/150‑1 (MQ=255)
                                                                                                                                 aGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTgg                       <  1:460151/150‑1 (MQ=255)
                                                                                                                                 aGTGGTGAGCGCGGTCAAGAACCTGAAAAACCGGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAGACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTgg                       <  1:679894/150‑1 (MQ=255)
                                                                                                                                      tgAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAaga                                                                                                                    <  1:1340854/52‑1 (MQ=255)
                                                                                                                                      tgAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAaga                                                                                                                    >  2:1340854/1‑52 (MQ=255)
                                                                                                                                      tgAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGc                  <  1:1666561/150‑1 (MQ=255)
                                                                                                                                       gAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAg                   >  1:610109/1‑148 (MQ=255)
                                                                                                                                       gAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAg                   <  2:610109/148‑1 (MQ=255)
                                                                                                                                       gAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCa                 <  1:474695/150‑1 (MQ=255)
                                                                                                                                            cgGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCaa                                                                                                                      <  2:1138961/44‑1 (MQ=255)
                                                                                                                                            cgGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCaa                                                                                                                      >  1:1138961/1‑44 (MQ=255)
                                                                                                                                             ggTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGt           <  1:1300486/150‑1 (MQ=255)
                                                                                                                                                  aGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTTGTGACTGACGACCGTGGCCAGCATGAAGTGCAAc      <  2:342638/150‑1 (MQ=255)
                                                                                                                                                     aCCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGtga                                    <  2:376887/117‑1 (MQ=255)
                                                                                                                                                     aCCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGtga                                    >  1:376887/1‑117 (MQ=255)
                                                                                                                                                      ccTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAg                                                   >  2:630026/1‑101 (MQ=255)
                                                                                                                                                      ccTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAg                                                   <  1:630026/101‑1 (MQ=255)
                                                                                                                                                      ccTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCa  <  2:1480495/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGGACGATCGCGTGCTGCACTGGATTGAAGTTGGCCAGCCCGATGCCGATCGCCTGACCTGGTGCTCGCGCCGTACCGAACGCACCAGCTTGCTGGCTTACGGCAGCCTGCGAGTGTGGGAGGGTAAAGTGGTGAGCGCGGTCAAGAACCTGAAAAACCTGAGCATTGCCGCTGTGCCGCAAGAGGTGCTGGAAACCCTGGCGACCGACATGCCGCGTAGCATCAAGTGGGACGTGATGATCAGTGAAGGTACGGTATTCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCA  >  NC_002947/1681192‑1681490

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: