Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,728,699 C→T 7.1% intergenic (+51/+72) 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,6990CT 7.1% 277.9 / 7.6 95intergenic (+51/+72)cspA‑I/PP_1523major cold shock protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/51);  new base T (3/4);  total (38/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCC  >  NC_002947/1728553‑1728845
                                                                                                                                                    |                                                                                                                                                  
gAGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGt                                                                                                                                                 <  2:395969/150‑1 (MQ=255)
 aGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTg                                                                                                                                                <  1:321083/150‑1 (MQ=255)
    ggCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCATCCAAGCCGAAAAAATAGTGGGCCTCTAACCCGCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCt                                                                                                                                             <  1:1434199/150‑1 (MQ=255)
      cATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGc                                                                                                                                           <  2:1531401/150‑1 (MQ=255)
      cATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGc                                                                                                                                           <  2:682011/150‑1 (MQ=255)
         cGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGAAGGGCCAGAAAGGCTTGCAAGCCGAAGACGTAGTGGGCCTATGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAGGCAGCCACCGTGGTGGCTGcttc                                                                                                                                        <  1:656335/150‑4 (MQ=255)
           tACCCTGGTCGAAGGGCTCAGGGTGCAGTCCTGCGGCCTGCAGGGCAAGAAAGGCCTGCAAGGCTAGCACGGTAGGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGTG‑‑GCAGCCACCGTGGTGGCTGcttcgcgg                                                                                                                                    <  1:443735/150‑8 (MQ=12)
           tACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACgg                                                                                                                                      <  1:446724/150‑1 (MQ=255)
            aCCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACggg                                                                                                                                     >  2:655756/1‑150 (MQ=255)
                  gTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGc                                                                                                                               >  1:1169124/1‑150 (MQ=255)
                    cGAAGGGCAGCGGGTGGACTACCCCTGCGTGCAGGGCCAGAAAGCCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTCGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGccc                                                                                                                             <  2:691650/150‑1 (MQ=255)
                                 gTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTa                                                                                                                <  2:925068‑M2/150‑2 (MQ=255)
                                 gTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTa                                                                                                                <  1:1453917‑M2/150‑2 (MQ=255)
                                   ggAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaaa                                                                                                              <  1:376137/150‑1 (MQ=255)
                                   ggAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaaa                                                                                                              <  2:985014/150‑1 (MQ=255)
                                    gAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  1:1050154/150‑1 (MQ=255)
                                    gAGTACGCCTGCGTGCAGGGCCAGAAAGGCATGCAGGCGCTGCAGGAGCTGGACCGCTGAGCCTCTGGTTTCCCGCGGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  1:1600029/150‑1 (MQ=255)
                                    gAGTACGCCTGCGTGCAGGGACAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  2:648513/150‑1 (MQ=255)
                                        aCGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCAt                                                                                                         <  1:1048375/150‑1 (MQ=255)
                                         cGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGccc                                                                                                                             <  2:1216754/129‑1 (MQ=255)
                                         cGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGccc                                                                                                                             >  1:1216754/1‑129 (MQ=255)
                                         cGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCAtt                                                                                                        <  2:1441435/150‑1 (MQ=255)
                                          gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTc                                                                                                       <  1:1054124/150‑1 (MQ=255)
                                          gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCAACGGCATTGATGCGCACTGAGCCGCTGGTTTCAAGAACAAGGTTCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTc                                                                                                       <  1:365534/150‑1 (MQ=255)
                                            cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTCGGAGCGGGCTTTCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                     >  2:1004824/1‑150 (MQ=255)
                                            cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTTCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                     >  2:497841/1‑150 (MQ=255)
                                            cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTcgc                                                                                                                     <  1:828297/134‑1 (MQ=255)
                                            cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTcgc                                                                                                                     >  2:828297/1‑134 (MQ=255)
                                            cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                     >  1:955839/1‑150 (MQ=255)
                                               cGTGCAGGGCCAGAAAGGCCTGCAAGGGGAAGACGTAGTGGGGCTCTGAGCCTCTGGTTTCAAGCAGTAGGAGCGAGTTTGATTAGGAAGAAGCTTCCTAGGTGGTTGTCACGGGCATTGACCGTATTCGCGGGTAAACCCATTCCTACa                                                                                                  <  2:297957/150‑1 (MQ=255)
                                                 tGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGa                                                                                                >  1:459177/1‑150 (MQ=255)
                                                    aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGt                                                                                                                                                 <  2:881406/98‑1 (MQ=255)
                                                    aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGt                                                                                                                                                 >  1:881406/1‑98 (MQ=255)
                                                      ggCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGa                                                                                           <  2:459177/150‑1 (MQ=255)
                                                      ggCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGa                                                                                           <  2:1453917/150‑1 (MQ=255)
                                                              aGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACgg                                                                                                                                      <  2:492857/99‑1 (MQ=255)
                                                              aGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACgg                                                                                                                                      >  1:492857/1‑99 (MQ=255)
                                                                 ccTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGaaa                                                                                              >  1:548433/1‑136 (MQ=255)
                                                                 ccTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGaaa                                                                                              <  2:548433/136‑1 (MQ=255)
                                                                  cTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCg                                                                               >  2:1066080/1‑150 (MQ=255)
                                                                   tGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGctt                                                                                                                                         <  2:1159241/91‑3 (MQ=39)
                                                                   tGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGctt                                                                                                                                         >  1:1159241/1‑89 (MQ=39)
                                                                   tGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGt                                                                              >  2:309790/1‑150 (MQ=255)
                                                                   tGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGt                                                                              >  1:1397886/1‑150 (MQ=255)
                                                                    gCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGtt                                                                             >  2:907701/1‑150 (MQ=255)
                                                                    gCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGtt                                                                             <  1:655756/150‑1 (MQ=255)
                                                                           gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCATCCAAGCAGCCACCATGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGCGCAAGCCTGCCGATTccttgt                                                                      >  1:611661/1‑145 (MQ=255)
                                                                               aCGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTAc                                                                  <  1:497841/150‑1 (MQ=255)
                                                                                  gTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                     <  1:813834/129‑1 (MQ=255)
                                                                                  gTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                     >  2:813834/1‑129 (MQ=255)
                                                                                        gggCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           <  1:1034406/150‑1 (MQ=255)
                                                                                        ggcccTCTCTGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGACACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           <  1:160983/147‑1 (MQ=255)
                                                                                        gatccTCTAAGCTTCGGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           <  2:1169124/147‑1 (MQ=255)
                                                                                         ggCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             >  1:665852/1‑99 (MQ=255)
                                                                                         ggCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  2:665852/99‑1 (MQ=255)
                                                                                              ctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc                                                     <  1:907701/150‑1 (MQ=255)
                                                                                                   ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             >  1:1146765/1‑89 (MQ=255)
                                                                                                   ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  2:1146765/89‑1 (MQ=255)
                                                                                                     tctGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAg                                                                                      >  2:579614/1‑110 (MQ=255)
                                                                                                     tctGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAg                                                                                      <  1:579614/110‑1 (MQ=255)
                                                                                                       tGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAAc                                                                                             <  1:891068/101‑1 (MQ=255)
                                                                                                       tGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAAc                                                                                             >  2:891068/1‑101 (MQ=255)
                                                                                                          tttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGt                                          >  2:548742/1‑149 (MQ=255)
                                                                                                          tttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGt                                          <  1:548742/149‑1 (MQ=255)
                                                                                                               aGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGGTGCCACGGGCACAGCCCGTATGCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCCgcgc                                    >  2:214152/1‑150 (MQ=255)
                                                                                                               aGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTAc                                                                                                   <  2:401440/87‑1 (MQ=255)
                                                                                                               aGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTAc                                                                                                   >  1:401440/1‑87 (MQ=255)
                                                                                                               aGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgcgc                                    >  2:63879/1‑150 (MQ=255)
                                                                                                                gCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCt                                   <  2:1128408/150‑1 (MQ=255)
                                                                                                                  cGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgat                                                                      >  1:846389/1‑113 (MQ=255)
                                                                                                                  cGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgat                                                                      <  2:846389/113‑1 (MQ=255)
                                                                                                                   gcaGGAGCGGGCTTGCCCGCGAAGCAGCCACCAAGCTGGCTGCCACGGTCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGAt                                <  2:611661/148‑1 (MQ=255)
                                                                                                                   gTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGttt                                                                            <  1:997237/106‑1 (MQ=255)
                                                                                                                   gTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGttt                                                                            >  2:997237/1‑106 (MQ=255)
                                                                                                                      ggAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCAtt                                                                                                        <  1:1067092/75‑1 (MQ=255)
                                                                                                                      ggAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCAtt                                                                                                        >  2:1067092/1‑75 (MQ=255)
                                                                                                                      ggAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCctc                                                  >  2:250574/1‑129 (MQ=255)
                                                                                                                      ggAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCctc                                                  <  1:250574/129‑1 (MQ=255)
                                                                                                                      ggAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCg                             >  2:1422167/1‑150 (MQ=255)
                                                                                                                       gAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             >  1:1527256/1‑69 (MQ=255)
                                                                                                                       gAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  2:1527256/69‑1 (MQ=255)
                                                                                                                        aGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGa                           <  1:1066080/150‑1 (MQ=255)
                                                                                                                         gCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgatg                                                                     >  2:929496/1‑107 (MQ=255)
                                                                                                                         gCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgatg                                                                     <  1:929496/107‑1 (MQ=255)
                                                                                                                         gCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGAc                          >  2:977286/1‑150 (MQ=255)
                                                                                                                            ggCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           <  1:341038/114‑1 (MQ=255)
                                                                                                                            ggCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           >  2:341038/1‑114 (MQ=255)
                                                                                                                              cTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaagcccgctcctacggcttgaaaccagaggctcagaggcccactacgtcttcggcttgcaggcctttctggccctgcacgcaggcgtactcc                     <  1:361671‑M2/150‑93 (MQ=255)
                                                                                                                               ttGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaagcccgctcctacggcttgaaaccagaggctcagaggcccactacgtcttcggcttgcaggcctttctggccctgcacgcaggcgtactcca                    <  2:575669‑M2/150‑94 (MQ=255)
                                                                                                                                  cccGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTc                 <  1:1004824/150‑1 (MQ=255)
                                                                                                                                   ccGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa                >  1:155030/1‑150 (MQ=255)
                                                                                                                                       gtggcagcCACCGTGGTGGCTGCCACACCCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAt            >  1:863049/4‑150 (MQ=255)
                                                                                                                                         agcagcCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtct          <  2:1397886/150‑1 (MQ=255)
                                                                                                                                           cagcCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtcttc        <  1:1422167/150‑1 (MQ=255)
                                                                                                                                              ccaccaCGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc     >  1:437520/1‑150 (MQ=255)
                                                                                                                                               caccacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaagcccgctcctacggcttgaaaccagaggctcagaggcccactacgtcttcggcttgcaggcctttctggccctgcacgcaggcgtactccacccgctgcccttcgac    >  1:575669‑M2/1‑41 (MQ=255)
                                                                                                                                               caccacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAGTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCCTCTTCTTcc    >  2:183681/1‑150 (MQ=255)
                                                                                                                                                accacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                     <  2:880304/68‑1 (MQ=255)
                                                                                                                                                accacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                     >  1:880304/1‑68 (MQ=255)
                                                                                                                                                 ccacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcccc  <  1:63879/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCC  >  NC_002947/1728553‑1728845

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