Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,716 A→C 6.0% C60G (TGC→GGC)  PP_1395 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,591,7160AC 6.0% 211.1 / 3.9 67C60G (TGC→GGC) PP_1395AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base C (2/2);  total (31/36)
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.86e-01

CGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAG  >  NC_002947/1591568‑1591863
                                                                                                                                                     |                                                                                                                                                   
cgcAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGca                                                                                                                                                     <  2:184572/149‑1 (MQ=255)
 gcAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGca                                                                                                                                                     <  6:280807/148‑1 (MQ=255)
    ggCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAgg                                                                                                                                                 <  5:281728/149‑1 (MQ=255)
       ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAATGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                                                >  1:214045/1‑147 (MQ=255)
       ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGTGCATCCTTGCCCAGGGCa                                                                                                                                              <  5:120471/149‑1 (MQ=255)
       ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGTGCATCCTTGCCCAGGGCa                                                                                                                                              <  1:250301/149‑1 (MQ=255)
       ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                                                >  2:210263/1‑147 (MQ=255)
        cAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAg                                                                                                                                             <  3:350789/149‑1 (MQ=255)
          ggCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGc                                                                                                                                           <  2:311238/149‑1 (MQ=255)
            cTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTg                                                                                                                                         <  3:79591/149‑1 (MQ=255)
             ttGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCTAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTgc                                                                                                                                        <  8:303916/149‑1 (MQ=255)
                  ttCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAg                                                                                                                                   <  6:128821/149‑1 (MQ=255)
                  ttCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAg                                                                                                                                   <  2:36911/149‑1 (MQ=255)
                       ggggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtgg                                                                                                                              <  5:188666/149‑1 (MQ=255)
                         ggCTGGCATCCAGCCTCATGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtggg                                                                                                                             >  3:319858/1‑148 (MQ=255)
                         ggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGCCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtgggt                                                                                                                            <  1:116729/149‑1 (MQ=255)
                                aTCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtgg                                                                                                                              <  3:327911/140‑1 (MQ=255)
                                aTCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAAc                                                                                                                     >  8:216165/1‑149 (MQ=255)
                                aTCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAAc                                                                                                                     >  6:54836/1‑149 (MQ=255)
                                aTCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAAc                                                                                                                     >  5:10374/1‑149 (MQ=255)
                                aTCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAAc                                                                                                                     >  8:116460/1‑149 (MQ=255)
                                aTCCACCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCAAGTAATGCAGTGGGTCGTAGGCGTACTCCTCGTCACCCAGGAACAAGCTTTTGCTGCCCTGGGCAAGGATGAACAGGGCAGCCTGCGCCAgtgtcg                                                                                                                              >  4:327911/1‑138 (MQ=255)
                                   cAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTa                                                                                                                  >  1:303722/1‑149 (MQ=255)
                                      ccTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACAt                                                                                                               >  2:127948/1‑149 (MQ=255)
                                        tCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTAc                                                                                                                 <  5:162997/145‑1 (MQ=255)
                                        tCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTAc                                                                                                                 >  6:162997/1‑145 (MQ=255)
                                          aCGGCGCTGAGTGGCAAGGTCACCGAAATGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTc                                                                                                           <  2:214045/149‑1 (MQ=255)
                                             gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                        <  5:102279/149‑1 (MQ=255)
                                             gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                        <  6:28229/149‑1 (MQ=255)
                                             gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                        <  6:10374/149‑1 (MQ=255)
                                                        cAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAAc                                                                                                   <  3:138585/143‑1 (MQ=255)
                                                        cAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAAc                                                                                                   >  4:138585/1‑143 (MQ=255)
                                                         aaGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAg                                                                                            >  2:141164/1‑149 (MQ=255)
                                                                         ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATG                                                                            <  7:216165‑M1/149‑1 (MQ=255)
                                                                             gTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCC                                                                         >  8:91943‑M1/1‑148 (MQ=255)
                                                                             gTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCC                                                                         >  4:124411‑M1/1‑148 (MQ=255)
                                                                                aTGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATg                                                                     <  1:141164/149‑1 (MQ=255)
                                                                                    aGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                        <  1:35764/110‑1 (MQ=255)
                                                                                    aGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                                        >  2:35764/1‑110 (MQ=255)
                                                                                          tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCGTTTTGCTGCCCTGGGCAAGGATGCACAGGGAAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCATCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                                             >  1:237029/1‑148 (MQ=255)
                                                                                            gTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAg                                                          <  1:210263/148‑1 (MQ=255)
                                                                                              aGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACa                                                                                       >  1:251547/1‑117 (MQ=255)
                                                                                              aGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACa                                                                                       <  2:251547/117‑1 (MQ=255)
                                                                                              aGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGc                                                       >  4:199451/1‑149 (MQ=255)
                                                                                                gCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTc                                                     >  4:246516/1‑149 (MQ=255)
                                                                                                  gTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAg                                                          >  1:180730/1‑142 (MQ=255)
                                                                                                  gTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAg                                                          <  2:180730/142‑1 (MQ=255)
                                                                                                   tACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGcc                                                  >  4:63490/1‑149 (MQ=255)
                                                                                                     cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGccc                                                 <  3:124411/148‑1 (MQ=255)
                                                                                                        cTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCAt                                             <  3:63490/149‑1 (MQ=255)
                                                                                                              aCCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGt                                       <  4:319858/149‑1 (MQ=255)
                                                                                                                ccAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCTATGCTGGACCAAAGGGCTCGCCCGCATGGCGGtt                                      >  1:206253/1‑148 (MQ=255)
                                                                                                                  aGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGTAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCATCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGAt                                   <  2:237029/149‑1 (MQ=255)
                                                                                                                        aGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTgg                                                                 <  4:337397/113‑1 (MQ=255)
                                                                                                                        aGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTgg                                                                 >  3:337397/1‑113 (MQ=255)
                                                                                                                         ggCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCtcg                            >  2:303722‑M2/1‑146 (MQ=255)
                                                                                                                         ggCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAAAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCtcg                            <  1:100982‑M2/149‑4 (MQ=255)
                                                                                                                          gCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCtcgg                           <  8:210280‑M2/149‑5 (MQ=255)
                                                                                                                               ttGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGAgctgatcaaccgccatgcgggcgagccctttggtccagcatcggcca                        <  7:116460‑M1/148‑48 (MQ=255)
                                                                                                                               agcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAg                       >  7:167126/4‑149 (MQ=255)
                                                                                                                                   tgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTc                    >  6:57989/1‑148 (MQ=255)
                                                                                                                                    gcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGAGGTTGATCAGCTCGTCCAGCTc                    >  5:368548/1‑147 (MQ=255)
                                                                                                                                     cccTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGTGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCtcg                            >  3:116009‑M2/1‑135 (MQ=255)
                                                                                                                                     cccTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCAGCAGGTACAGGTCATCGATGGCGGATGCAGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCCGCtcg                            <  4:116009‑M2/138‑4 (MQ=255)
                                                                                                                                           gcatccttGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGCTCAGCTCGGCCAGCTCCAGGCGCTg           >  6:178351/8‑149 (MQ=255)
                                                                                                                                            cAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTg           >  3:230233/1‑148 (MQ=255)
                                                                                                                                            cAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGc          <  7:91943/149‑1 (MQ=255)
                                                                                                                                             aaGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCa         <  5:54836/149‑1 (MQ=255)
                                                                                                                                                  tGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTcc    <  3:199451/149‑1 (MQ=255)
                                                                                                                                                    cacaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTcc    <  6:274313/147‑1 (MQ=255)
                                                                                                                                                    cacaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAg  >  3:213040/1‑149 (MQ=255)
                                                                                                                                                    cacaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAg  >  3:217878/1‑149 (MQ=255)
                                                                                                                                                    cacaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGAGGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTcc    >  5:274313/1‑147 (MQ=255)
                                                                                                                                                     acaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGAgctgatcaaccgccatgcgggcgagccctttggtccagcatcgg                           <  4:52641‑M1/123‑45 (MQ=255)
                                                                                                                                                     acaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGAgctgatcaaccgccatgcgggcgagccctttggtccagcatcgg                           >  3:52641‑M1/1‑79 (MQ=255)
                                                                                                                                                     acaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAgg                >  6:316881/1‑134 (MQ=255)
                                                                                                                                                     acaGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAgg                <  5:316881/134‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTGCAGTTCCAG  >  NC_002947/1591568‑1591863

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