Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,450,925 C→T 15.0% R37W (CGG→TGG)  PP_1270 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,450,9250CT15.0% 203.9 / 23.0 86R37W (CGG→TGG) PP_1270LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base T (6/7);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CTACCATCGCGTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGC  >  NC_002947/1450778‑1451068
                                                                                                                                                   |                                                                                                                                                   
cTACCATCGCGTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGcc                                                                                                                                                     <  3:288346/148‑1 (MQ=255)
    cATCGCGTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACa                                                                                                                                                >  2:295991/1‑149 (MQ=255)
      tCGCGTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACAtt                                                                                                                                             >  3:53313/1‑149 (MQ=255)
      tCGCGTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACAtt                                                                                                                                             >  7:317886/1‑149 (MQ=255)
         cgTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTggg                                                                                                                                          <  6:190418/149‑1 (MQ=255)
         cgTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTggg                                                                                                                                          <  5:411992/149‑1 (MQ=255)
           tttCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGc                                                                                                                                        >  5:166055/1‑149 (MQ=255)
           tttCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGc                                                                                                                                        <  7:244345/149‑1 (MQ=255)
                cGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGTGCGTATc                                                                                                                                   >  3:326487/1‑149 (MQ=255)
                 gCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTggg                                                                                                                                          <  4:315635/141‑1 (MQ=255)
                 gCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTggg                                                                                                                                          >  3:315635/1‑141 (MQ=255)
                           cGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCTGGACATTgg                                                                                                                                           >  8:237378/1‑130 (MQ=255)
                           cGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCTGGACATTgg                                                                                                                                           <  7:237378/130‑1 (MQ=255)
                            ggAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTgg                                                                                                                                           >  6:245270/1‑129 (MQ=255)
                            ggAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTgg                                                                                                                                           <  5:245270/129‑1 (MQ=255)
                            ggAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTg                                                                                                                       >  2:315274/1‑149 (MQ=255)
                            ggAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTg                                                                                                                       >  6:349897/1‑149 (MQ=255)
                             gAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTgg                                                                                                                      <  8:287494/149‑1 (MQ=255)
                             gAAAACGGACATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGTTGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTgg                                                                                                                      <  6:226201/149‑1 (MQ=255)
                                   ggCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGAcc                                                                                                                <  2:173822/149‑1 (MQ=255)
                                     ccATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCg                                                                                                              <  6:166055/149‑1 (MQ=255)
                                           aaCTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTg                                                                                                        <  8:183169/149‑1 (MQ=255)
                                             cTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGt                                                                                                      <  1:295991/149‑1 (MQ=255)
                                              tCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTg                                                                                                        <  6:37784/146‑1 (MQ=255)
                                              tCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTg                                                                                                        >  5:37784/1‑146 (MQ=255)
                                              tCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTg                                                                                                     >  5:178582/1‑149 (MQ=255)
                                              tCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTg                                                                                                     >  3:301580/1‑149 (MQ=255)
                                              tCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTAGGGTGACCCGATCCTGGTg                                                                                                     >  8:126348/1‑149 (MQ=255)
                                                 cGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTgcgc                                                                                                  >  7:205753/1‑149 (MQ=255)
                                                      aTGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCgccgg                                                                                             >  3:241263/1‑149 (MQ=255)
                                                        gAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCgccgg‑‑cc                                                                                         >  3:343566/1‑149 (MQ=255)
                                                          aCCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGt                                                                                       >  7:192952/1‑149 (MQ=255)
                                                           ccTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCGGG‑‑CCGTg                                                                                     >  1:140824/1‑149 (MQ=255)
                                                           ccTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTg                                                                                     >  5:262168/1‑149 (MQ=255)
                                                               cTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATCTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGc                                                                                   <  1:315274/148‑1 (MQ=255)
                                                               cTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGcc                                                                                  <  4:241263/148‑1 (MQ=255)
                                                                  gtcgCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTgg                                                                               <  2:246180/148‑1 (MQ=255)
                                                                  gtcgCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTg                                                                                <  6:178582/147‑1 (MQ=255)
                                                                       ccTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTa                                                                         <  2:9918/149‑1 (MQ=255)
                                                                       ccTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTa                                                                         >  3:128024/1‑149 (MQ=255)
                                                                       ccTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTa                                                                         <  8:317886/149‑1 (MQ=255)
                                                                         tCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACg                                                                       >  8:47127/1‑149 (MQ=255)
                                                                           gATGCCTTGCTCGACGAGGGGAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGGCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGcc                                                                     <  6:134255/149‑1 (MQ=255)
                                                                              gCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCg                                                                  >  7:382977/1‑149 (MQ=255)
                                                                              gCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGAGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCg                                                                  >  1:59433/1‑149 (MQ=255)
                                                                                cTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGCCATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTg                                                                >  3:5314/1‑149 (MQ=255)
                                                                                cTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGta                                                                >  4:26759/1‑148 (MQ=255)
                                                                                cTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTg                                                                >  5:275239/1‑149 (MQ=255)
                                                                                     gcgacgaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTg                                                           <  8:308085/148‑1 (MQ=255)
                                                                                       gacgaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCAACGCACGGGGCGTAGGc                                                         <  3:194004/149‑1 (MQ=255)
                                                                                         cgaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGccc                                                       >  2:355970/1‑149 (MQ=255)
                                                                                          gaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCt                                                      <  8:192952/149‑1 (MQ=255)
                                                                                              gCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTgcgc                                                  >  3:240320/1‑149 (MQ=255)
                                                                                              gCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTgcgc                                                  <  7:47127/149‑1 (MQ=255)
                                                                                               cAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTgcgc                                                  >  4:136356/1‑148 (MQ=255)
                                                                                                aGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGa                                                <  4:343566/149‑1 (MQ=255)
                                                                                                 gCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAg                                               <  4:240320/149‑1 (MQ=255)
                                                                                                    tggtggGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAg                                            <  6:172687/149‑1 (MQ=255)
                                                                                                      gtggGCGCAGCGCAGCGCTTGAACGTGTGCGCGGCGGGGATGAGCCGGAGTTTGGGGCGTATGCGCGATGCCTTGGGTGTCCCGTTCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGGCTACGCGCCGTGCGTTGGGCCTGCGCGAGCAGGt                                          <  2:59433/149‑1 (MQ=255)
                                                                                                      gtggGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGt                                          <  5:358323/149‑1 (MQ=255)
                                                                                                      gtggGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGt                                          >  6:358323/1‑149 (MQ=255)
                                                                                                        gggCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGc                                        <  1:355970/149‑1 (MQ=255)
                                                                                                           cgcagcgcagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATg                                     >  3:158274/1‑149 (MQ=255)
                                                                                                           cgcagcgcagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATg                                     <  4:128024/149‑1 (MQ=255)
                                                                                                            gcagcgcagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATgg                                    >  6:263855/1‑149 (MQ=255)
                                                                                                                 gcagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTCGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtg                                <  1:309461/149‑1 (MQ=255)
                                                                                                                 gcagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtgg                               >  8:12350/1‑149 (MQ=255)
                                                                                                                   agcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGCGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtggtg                             <  5:349897/149‑1 (MQ=255)
                                                                                                                     cgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGa                           <  7:126348/149‑1 (MQ=255)
                                                                                                                        aggAACCTGAGGCCGGCGGCCATGTGCCGGACATTGGGGCGTATCCGGGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCa                        <  4:5314/147‑1 (MQ=255)
                                                                                                                        aTGAACCTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAg                          <  1:120852/149‑1 (MQ=255)
                                                                                                                             ccTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtgg                               <  4:77856/137‑1 (MQ=255)
                                                                                                                             ccTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtgg                               >  3:77856/1‑137 (MQ=255)
                                                                                                                             ccTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCg                   <  6:262168/149‑1 (MQ=255)
                                                                                                                              cTGAGCCCGGCGGCCATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCgg                  >  3:263524/1‑149 (MQ=255)
                                                                                                                                gAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTg                                                                                <  6:3670/85‑1 (MQ=255)
                                                                                                                                gAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTg                                                                                >  5:3670/1‑85 (MQ=255)
                                                                                                                                 aGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTGGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCggg                 <  4:53313/147‑1 (MQ=255)
                                                                                                                                  gCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAg              >  5:18011/1‑149 (MQ=255)
                                                                                                                                      ggcggcCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCTCGAGCAGGTGCATGGCCTGTTGGAGCAGGCCGGGCAGgtgt          <  6:77956/149‑1 (MQ=255)
                                                                                                                                       gcggcCATGAGCCGGACATGGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGtt         <  4:301580/149‑1 (MQ=255)
                                                                                                                                        cggcCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTc        <  4:158274/149‑1 (MQ=255)
                                                                                                                                         ggcCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGGGGCCTGGTGCCTACGCCCCGGGCGTTGGCTCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTcc       <  8:382977/149‑1 (MQ=255)
                                                                                                                                            cATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCa                        <  2:172870/129‑1 (MQ=255)
                                                                                                                                            cATGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCa                        >  1:172870/1‑129 (MQ=255)
                                                                                                                                            aaTGAGCCGGACATTGGGGCGTATCCGCGCTGCCTTGGGGGACCCGATCCTGGTGCGCGCCGG‑‑CAGTGGACGGGTTCCTACGCCCCTTGAGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGGGGATCATGCCGGGCAGGTGTTCCGCa    <  4:326487/148‑1 (MQ=255)
                                                                                                                                              tGAGCTGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGc  >  2:250244/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTACCATCGCGTTTCCCGCAGCTGTAACGGAAAACGGCCATGCAACTCCCGGACATGAACCTGCTCGTCGCCCTCGATGCCTTGCTCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGG‑‑CCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGC  >  NC_002947/1450778‑1451068

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: