Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,761,211 T→C 6.2% I88T (ATC→ACC)  gpmI → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,761,2110TC 6.2% 187.8 / 2.6 64I88T (ATC→ACC) gpmI2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (2/2);  total (33/31)
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.75e-01

ACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGC  >  NC_002947/5761069‑5761352
                                                                                                                                               |                                                                                                                                             
aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGa                                                                                                                                         <  2:265614/149‑1 (MQ=255)
aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGa                                                                                                                                         >  7:23020/1‑149 (MQ=255)
aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGATGGCCTTGGTCACCcgcgt                                                                                                                                         <  8:81127/149‑2 (MQ=255)
aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGATGGCCTTGGTCACCcgcgt                                                                                                                                         <  7:130614/149‑2 (MQ=255)
  ccAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACg                                                                                                                                       >  1:338948/1‑149 (MQ=255)
  ccAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACg                                                                                                                                       >  8:154178/1‑149 (MQ=255)
   cAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACgg                                                                                                                                      >  1:31126/1‑149 (MQ=255)
      ccGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCTACGGCGa                                                                                                                                   <  7:154178/149‑1 (MQ=255)
             gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            <  4:153795/149‑1 (MQ=255)
             gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            <  6:194379/149‑1 (MQ=255)
             gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            <  5:333513/149‑1 (MQ=255)
              ccTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            >  4:147296/1‑148 (MQ=255)
              ccTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            <  3:147296/148‑1 (MQ=255)
                  aTTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGaaa                                                                                                                        <  1:296389/148‑1 (MQ=255)
                   tttCCGGCTCGGGCATGGATGTCGGCCTGACGGAAGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAAc                                                                                                                      >  6:332084/1‑149 (MQ=255)
                         gCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTa                                                                                                                >  3:169490/1‑149 (MQ=255)
                             gggCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGa                                                                                                            <  8:23020/149‑1 (MQ=255)
                               gCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                          >  7:246526/1‑149 (MQ=255)
                               gCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                          >  2:77570/1‑149 (MQ=255)
                               gCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                          >  6:125385/1‑149 (MQ=255)
                                   ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                           >  8:215085/1‑144 (MQ=255)
                                   ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                           <  7:215085/144‑1 (MQ=255)
                                   ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCgg                                                                                                        <  3:203877/147‑1 (MQ=255)
                                   ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCgg                                                                                                        >  4:203877/1‑147 (MQ=255)
                                        tCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgt                                                                                                  <  5:166059/148‑1 (MQ=255)
                                          ggCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACAGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGa                                                                                               >  3:238509/1‑149 (MQ=255)
                                             cTGGCGGACGGGCAGATGGGCAACTCCGTAGTCGTTCATATGATAGTGGGCGCAGGTCGCGTTGTGTATGAGGACTTCACTCTGGTGTCCATGGCCATGCGTGACTGCGTGTTCTTGTATAAACCTGTAGTGTCCGGTGCTGTTGatta                                                                                            <  4:238509/149‑4 (MQ=255)
                                                ccGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAgg                                                                                          <  2:338948/148‑1 (MQ=255)
                                                ccGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAg                                                                                           <  2:131184/147‑1 (MQ=255)
                                                ccGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAg                                                                                           >  1:131184/1‑147 (MQ=255)
                                                 cGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCa                                                                                        <  3:120010/149‑1 (MQ=255)
                                                     cGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCg                                                                                    <  4:169490/149‑1 (MQ=255)
                                                        gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                          >  6:322479/1‑124 (MQ=255)
                                                        gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                          <  5:322479/124‑1 (MQ=255)
                                                        gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgt                                                                                                  >  5:117637/1‑132 (MQ=255)
                                                        gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgt                                                                                                  <  6:117637/132‑1 (MQ=255)
                                                         cAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGc                                                                                <  2:31126/149‑1 (MQ=255)
                                                          aGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCt                                                                               >  2:188969/1‑149 (MQ=255)
                                                          aGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGAGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCt                                                                               >  8:210021/1‑149 (MQ=255)
                                                             tGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACaa                                                                                            <  3:198910/133‑1 (MQ=255)
                                                             tGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACaa                                                                                            >  4:198910/1‑133 (MQ=255)
                                                              gggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            <  6:14475/100‑1 (MQ=255)
                                                              gggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                            >  5:14475/1‑100 (MQ=255)
                                                               ggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCaa                                                                          <  5:125385/149‑1 (MQ=255)
                                                               ggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCaa                                                                          <  5:332084/149‑1 (MQ=255)
                                                                gCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAg                                                                         >  5:282555/1‑149 (MQ=255)
                                                                                      gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTACACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGActg                                                    >  1:56489/1‑149 (MQ=255)
                                                                                      gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAgg                                                                                          >  6:265739/1‑110 (MQ=255)
                                                                                      gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAgg                                                                                          <  5:265739/110‑1 (MQ=255)
                                                                                               cgcAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtg                                          >  8:164326/1‑149 (MQ=255)
                                                                                                  aGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                       <  1:77570/149‑1 (MQ=255)
                                                                                                             tgtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGc                            >  4:229458/1‑149 (MQ=255)
                                                                                                               tATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                       <  8:286290/136‑1 (MQ=255)
                                                                                                               tATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                       >  7:286290/1‑136 (MQ=255)
                                                                                                                aTCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                         >  5:48969/1‑149 (MQ=255)
                                                                                                                 tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCCGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                         <  8:246526/148‑1 (MQ=255)
                                                                                                                        ttCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCa                 >  3:189472/1‑149 (MQ=255)
                                                                                                                          cgtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCAc                >  7:19803/4‑149 (MQ=255)
                                                                                                                           gtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACt               >  1:292940/3‑149 (MQ=255)
                                                                                                                                 ggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg         <  7:164326/149‑1 (MQ=255)
                                                                                                                                      ccAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATg    >  8:268669/1‑149 (MQ=255)
                                                                                                                                       cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg         <  2:131244/143‑1 (MQ=255)
                                                                                                                                       cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg         >  1:131244/1‑143 (MQ=255)
                                                                                                                                        aaGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGc  >  7:8303/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                             
ACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGGC  >  NC_002947/5761069‑5761352

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: