Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,728,698 A→G 9.1% intergenic (+50/+73) cspA‑I → / ← PP_1523 major cold shock protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,728,6980AG 9.1% 166.3 / 5.5 63intergenic (+50/+73)cspA‑I/PP_1523major cold shock protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base G (4/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

TGAGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCC  >  NC_002947/1728552‑1728844
                                                                                                                                                  |                                                                                                                                                  
tGAGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACgg                                                                                                                                                 <  2:718749/150‑1 (MQ=255)
 taGGTGCATCGTACCCTGGTCGAAGTGCAGCGGGTGCATTACGCCTGCGTCCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGCCCTCTGATCCTCTGGTTTCAAGCCGTAGGATCGGTCTTGCCCGCGAAGCAGCCACCACGGt                                                                                                                                                <  2:471406/149‑1 (MQ=255)
  aGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTg                                                                                                                                               <  1:250608/150‑1 (MQ=255)
     ggCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCt                                                                                                                                            <  1:90396/150‑1 (MQ=255)
              caccggTCGAACAGCAGCGGATCGAGTAAGCATCCGCGCAGGTCCAGGAAGGCCTGCAGGGCGAAGGCGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCGCGCCCGCGAAGCAGCGACGGCCGGGGCTGGCTCGGGc                                                                                                                                   >  2:15885/5‑150 (MQ=255)
               ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCTCGCGAAACAGCCACCACGGTGGCTGCCACGGGCa                                                                                                                                  >  2:648856/1‑150 (MQ=255)
                  ggTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCaagc                                                                                                                               <  1:842617/150‑4 (MQ=255)
                      gAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCg                                                                                                                           >  1:688623/1‑150 (MQ=255)
                            cAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCACCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTcg                                                                                                                     >  1:359839/1‑150 (MQ=255)
                             aGCGGGTGGAGTACGCCTGCGTGCAGGGCAAGAAAGGCCTGCAAGCCGAAGGCGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTcgc                                                                                                                    <  2:866979/150‑1 (MQ=255)
                               cGGGTGGAGTACGCCTGCGTGCAGGTCCAGCACGGCTTGCAATCCGAAGAGGTAGGGGGCCTCTGAGCCTCTGGTTCCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCgg                                                                                                                  <  2:474057/150‑1 (MQ=255)
                                     gAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATACGCGGGTAAAc                                                                                                            >  1:1287793/1‑150 (MQ=255)
                                         aCGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCAt                                                                                                        >  2:1306224/1‑150 (MQ=255)
                                         aCGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCAt                                                                                                        <  2:28611/150‑1 (MQ=255)
                                             cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACgg                                                                                                                                     <  1:393499/117‑1 (MQ=255)
                                             cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACgg                                                                                                                                     >  2:393499/1‑117 (MQ=255)
                                             cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                    >  1:555213/1‑150 (MQ=255)
                                              tGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGCCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTa                                                                                                   <  1:15885/150‑1 (MQ=255)
                                                     aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCCAAGCAGCCACCACGGTGGCTGCCCCGGGCATAGCCCGTACTCCCGGGTAAACCCATTCCTACAGaaaa                                                                                            >  2:35819/1‑149 (MQ=255)
                                                           aGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACggcagg                                                                                      >  2:1204137/1‑145 (MQ=255)
                                                           aGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCACGTATTCGCGGTGAAACCCATTCCTACAGAAACGAGCaa                                                                                      <  1:629894/150‑1 (MQ=255)
                                                             aaaGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                    <  2:688623/150‑1 (MQ=255)
                                                             aaaGGCCTGCAAGACGAAGACGTAGTGGGCCTCTGGGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCAAGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                    <  1:648856/150‑1 (MQ=255)
                                                                 gCCTGCAAGCCGCAGACGTAGTGAGCCTCTGCGCCTCCGGTTTCAATCAGTATTTCCCGACTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCGTACAGAAACGAGCCGGCCttc                                                                                <  2:87908/150‑3 (MQ=255)
                                                                   cTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCgg                                                                                                                  >  1:761967/1‑114 (MQ=255)
                                                                   cTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCgg                                                                                                                  <  2:761967/114‑1 (MQ=255)
                                                                     gCAAGCCGAAGACGTAGTGGGCCTCTGAGACTCTGATTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATATCCCGTATTCTCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGtt                                                                            >  1:8575/1‑150 (MQ=255)
                                                                       aaGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTc                                                                          >  1:537107/1‑150 (MQ=255)
                                                                        aGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCa                                                                         >  2:1305121/1‑150 (MQ=255)
                                                                         gCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCat                                                                        <  2:436351/150‑1 (MQ=255)
                                                                          ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCTCGGGTAAAc                                                                                                            >  2:1327700/1‑113 (MQ=255)
                                                                          ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                            <  1:1327700/113‑1 (MQ=255)
                                                                            gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGc                                                                                                                              >  2:1342296/1‑93 (MQ=255)
                                                                            gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGc                                                                                                                              <  1:1342296/93‑1 (MQ=255)
                                                                            gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgat                                                                     <  1:1305121/150‑1 (MQ=255)
                                                                            gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgat                                                                     >  2:544630/1‑150 (MQ=255)
                                                                             aaGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgatg                                                                    <  2:1202283/150‑1 (MQ=255)
                                                                             aaGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgatg                                                                    <  2:1202761/150‑1 (MQ=255)
                                                                              aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                            >  1:1079141/1‑109 (MQ=255)
                                                                              aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                            <  2:1079141/109‑1 (MQ=255)
                                                                                  gTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGc                                                               >  2:879143/1‑150 (MQ=255)
                                                                                       gggCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                          <  1:64611/150‑1 (MQ=255)
                                                                                         gCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGa                                                        >  1:521645/1‑150 (MQ=255)
                                                                                          ctcctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAt                                                       <  2:931253/147‑1 (MQ=255)
                                                                                           ctctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAtt                                                      >  1:330401/1‑150 (MQ=255)
                                                                                            tctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAttt                                                     >  1:741890/1‑150 (MQ=255)
                                                                                             gtGAGCCTCGGGTTCCACGCCGTCGGAGAGAGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc                                                    <  2:299647/149‑1 (MQ=255)
                                                                                                aGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCctc                                                 <  2:330401/150‑1 (MQ=255)
                                                                                                        gTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGcc                                                              <  1:699393/129‑1 (MQ=255)
                                                                                                        gTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGcc                                                              >  2:699393/1‑129 (MQ=255)
                                                                                                               gCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCt                                  >  1:821383/1‑150 (MQ=255)
                                                                                                                   tAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATg                              <  1:544630/150‑1 (MQ=255)
                                                                                                                         cGGGCTTACCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCAGCATAATCAGGCAAACTCAGGCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGAcc                        <  1:879143/150‑1 (MQ=255)
                                                                                                                                  ccGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCGGCAATGAAGGAGCATGAGGTGCGCCAGGTGATTTGATATTCACCGACGTCTCTGATGCGGTCCCAGCGGTCa               >  1:1096165/1‑150 (MQ=255)
                                                                                                                                  ccGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa               >  2:571118/1‑150 (MQ=255)
                                                                                                                                  ccGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa               >  1:1195272/1‑150 (MQ=255)
                                                                                                                                    ccgtggcagcCACCGTGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCa                                                             <  2:233900/97‑1 (MQ=255)
                                                                                                                                    ccgtggcagcCACCGTGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCa                                                             >  1:233900/6‑102 (MQ=255)
                                                                                                                                        ggcagcCACCGTGGTGGCTGCTACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAtt                                                      <  2:756769/104‑1 (MQ=255)
                                                                                                                                        ggcagcCACCGTGGTGGCTGCTACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAtt                                                      >  1:756769/2‑105 (MQ=255)
                                                                                                                                         gcagcCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtctt        <  2:821383/150‑1 (MQ=255)
                                                                                                                                          cagcCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtcttc       <  1:35819/150‑1 (MQ=255)
                                                                                                                                           agcCACCGTGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtcttct      >  2:1184887/1‑150 (MQ=255)
                                                                                                                                              caccacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcc   <  2:521645/150‑1 (MQ=255)
                                                                                                                                               accacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTccc  >  2:641135/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGAGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCC  >  NC_002947/1728552‑1728844

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