Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,027,746 T→A 7.9% E228D (GAA→GAT)  PP_3552 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,027,7460TA 7.9% 178.1 / 5.5 63E228D (GAA→GAT) PP_3552two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base A (2/3);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

ATCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGT  >  NC_002947/4027599‑4027894
                                                                                                                                                   |                                                                                                                                                        
aTCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCtt                                                                                                                                                         >  7:252553/1‑149 (MQ=255)
         cAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATgcgc                                                                                                                                                 >  5:364182/1‑148 (MQ=255)
          aGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCATCGAAGCGCag                                                                                                                                               <  7:350392/149‑2 (MQ=255)
           ggCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAg                                                                                                                                              <  6:364182/149‑1 (MQ=255)
            gCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAgg                                                                                                                                             <  6:180767/149‑1 (MQ=255)
               cTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGa                                                                                                                                          <  5:329684/149‑1 (MQ=255)
                tGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGAt                                                                                                                                         >  2:53036/1‑149 (MQ=255)
                     gTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTcg‑‑‑‑c                                                                                                                                >  2:300753/1‑149 (MQ=255)
                     gTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGACGCGCAAGGTGATGTcg‑‑‑‑c                                                                                                                                >  2:185181/1‑149 (MQ=255)
                      tAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCATCGAAGCGCa                                                                                                                                                >  2:34162/1‑136 (MQ=255)
                      tAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGAGCGCCAGGCTTTCCTGCGCATCGAAGCGCa                                                                                                                                                <  1:34162/136‑1 (MQ=255)
                        aCGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑Cgg                                                                                                                              <  6:104485/148‑1 (MQ=255)
                              aCTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCACGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACg                                                                                                                        >  6:356079/1‑148 (MQ=255)
                                    gTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGacca                                                                                                                 >  6:365070/1‑149 (MQ=255)
                                        gCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTc                                                                                                             >  1:116503/1‑149 (MQ=255)
                                          aTAGTTGAGCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCAGAGGCCAGCTCGCCGAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGGGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACGACCTcca                                                                                                           <  4:257935/149‑1 (MQ=255)
                                          aTAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTcca                                                                                                           >  7:260878/1‑149 (MQ=255)
                                           tAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTccac                                                                                                          <  5:365070/149‑1 (MQ=255)
                                            aGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTccac                                                                                                          >  7:35353/1‑148 (MQ=255)
                                                   caccGCGGCCAGTGGCTGATTGTTCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACcaccgcac                                                                                                  >  1:251689/1‑149 (MQ=255)
                                                   caccGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTc                                                                                                             >  6:243221/1‑138 (MQ=255)
                                                   caccGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTc                                                                                                             <  5:243221/138‑1 (MQ=255)
                                                         ggCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGgtgt                                                                                            <  3:299053/149‑1 (MQ=255)
                                                           ccAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAg                                                                                          <  2:220734/149‑1 (MQ=255)
                                                            cAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAgg                                                                                         >  8:59874/1‑149 (MQ=255)
                                                                   tgattgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAgg                                                                                                                                             <  1:254904/94‑1 (MQ=255)
                                                                   tgattgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAgg                                                                                                                                             >  2:254904/1‑94 (MQ=255)
                                                                        gatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACAACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCg                                                                             >  8:100419/1‑149 (MQ=255)
                                                                          tctcGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACAAGCTGACAACCTCCACCACCGCACCGGTCTAGGTATCACGGATCGcc                                                                           >  5:48632/1‑149 (MQ=255)
                                                                                 ggCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTg                                                                     <  8:35353/148‑1 (MQ=255)
                                                                                  gCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                                   <  5:356079/149‑1 (MQ=255)
                                                                                  gCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                                   <  1:76732/149‑1 (MQ=255)
                                                                                       gCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCTCCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGAttt                                                              >  4:243733/1‑149 (MQ=255)
                                                                                        cccGAGGCGAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTc                                                             <  4:144012/149‑1 (MQ=255)
                                                                                        cccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTc                                                             <  2:116503/149‑1 (MQ=255)
                                                                                         ccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTc                                                             >  2:250063/1‑148 (MQ=255)
                                                                                           gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCTCCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGaa                                                          <  3:243733/149‑1 (MQ=255)
                                                                                           gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGaa                                                          <  1:53036/149‑1 (MQ=255)
                                                                                              gCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAAcc                                                        >  3:239824/1‑148 (MQ=255)
                                                                                               ccAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCCTCCTGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGaa                                                          <  1:300753/149‑1 (MQ=255)
                                                                                                 aGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACg                                                    <  2:251689/149‑1 (MQ=255)
                                                                                                  gCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                   <  1:250063/149‑1 (MQ=255)
                                                                                                  gCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                   <  1:136178/149‑1 (MQ=255)
                                                                                                    tCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                   >  1:41734/1‑147 (MQ=255)
                                                                                                     cgcgCAGGGTAGAGAGTCGGGCTCTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCAGCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTa                                                 <  6:48632/148‑1 (MQ=255)
                                                                                                      gCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGACGATTTCGAACCAACGGTa                                                 <  5:50784/147‑1 (MQ=255)
                                                                                                            ggTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctcct                                          <  8:260878/148‑1 (MQ=255)
                                                                                                                gCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGACGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGc                                     <  1:185181/149‑1 (MQ=255)
                                                                                                                   aGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTg                                  <  3:166311/149‑1 (MQ=255)
                                                                                                                    ggCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTgg                                 >  5:330854/1‑149 (MQ=255)
                                                                                                                                gcgcCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCa                      <  4:239824/148‑1 (MQ=255)
                                                                                                                                  gcCCGGCTTTCCTGCGCTTCGATGCGCCCGGTGATGTCG‑‑‑‑CGGGACACGCTGCCCACCTCCACCACCGCACCGGTGTCGATCTCACTGATCGCCCTAAAATAAAAAAAAAAAAAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAgg                    >  6:143968/1‑148 (MQ=255)
                                                                                                                                  gcCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                                   <  4:253941/101‑1 (MQ=255)
                                                                                                                                    cAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGGAGCCTCCTGTCGCCTGGCGTACCCGGCAGGGCa                 >  8:333026/1‑149 (MQ=255)
                                                                                                                                    cAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGAAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCa                 >  5:270178/1‑149 (MQ=255)
                                                                                                                                      ggCTTTCCTGCGCTTCGATGCGCAAGGTGATGGCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc               <  8:252553/149‑1 (MQ=255)
                                                                                                                                       gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGc                                >  2:138555/1‑131 (MQ=255)
                                                                                                                                       gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGc                                <  1:138555/131‑1 (MQ=255)
                                                                                                                                           ccttGCGCATCGAAGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg          >  4:17920/4‑149 (MQ=255)
                                                                                                                                             cTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc               >  5:45548/1‑142 (MQ=255)
                                                                                                                                             cTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc               <  6:45548/142‑1 (MQ=255)
                                                                                                                                                cgcATCGAAGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAg      <  3:17920/148‑1 (MQ=255)
                                                                                                                                                  cTTCGATGCGCAAGGAGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                   >  5:232904/1‑101 (MQ=255)
                                                                                                                                                  cTTCGATGCGCAAGGAGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                   <  6:232904/101‑1 (MQ=255)
                                                                                                                                                   ttCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGt  >  4:84695/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
ATCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGT  >  NC_002947/4027599‑4027894

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