Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,525,770 C→A 9.1% intergenic (+60/‑63) PP_4858 → / → PP_4859 hypothetical protein/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,525,7700CA 9.1% 192.3 / 7.6 65intergenic (+60/‑63)PP_4858/PP_4859hypothetical protein/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base A (2/4);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

CGGCTCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGT  >  NC_002947/5525626‑5525913
                                                                                                                                                |                                                                                                                                               
cggcTGCTATCGGCTTATCGAACTATCTTGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCAtt                                                                                                                                              <  7:274409/148‑1 (MQ=255)
   cTCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGAcc                                                                                                                                          <  4:68536/149‑1 (MQ=255)
      cGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCgg                                                                                                                                        <  6:51173/148‑1 (MQ=255)
       gATGGTCTGATCGAACTAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                      <  7:248282/149‑1 (MQ=255)
       gATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                      <  5:18150/149‑1 (MQ=255)
       gATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                      <  8:53747/149‑1 (MQ=255)
         tCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCtt                                                                                                                                    <  5:242113/149‑1 (MQ=255)
               gATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGtttt                                                                                                                              >  5:269097/1‑149 (MQ=255)
                aTCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                             >  4:361715/1‑149 (MQ=255)
                aTCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                             >  2:60891/1‑149 (MQ=255)
                aTCGAACAAGCTGGCCGTGTAGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTTTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttt                                                                                                                               <  8:5416/147‑1 (MQ=255)
                    aaCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTAt                                                                                                                         >  2:151757/1‑149 (MQ=255)
                       aaGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATc                                                                                                                        >  4:355442/1‑147 (MQ=255)
                         gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTaga                                                                                                                    <  1:60891/149‑1 (MQ=255)
                         gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTa                                                                                                                      >  6:183269/1‑147 (MQ=255)
                         gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTa                                                                                                                      <  4:266687/147‑1 (MQ=255)
                           tGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCa                                                                                                                                                <  2:193569/119‑1 (MQ=255)
                           tGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCa                                                                                                                                                >  1:193569/1‑119 (MQ=255)
                              ccGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGcc                                                                                                               <  2:261792/149‑1 (MQ=255)
                                      cTGAAGGTGGCTCGGATCGCCTGGTTTAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGa                                                                                                       <  3:361715/149‑1 (MQ=255)
                                      cTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGa                                                                                                       >  8:297738/1‑149 (MQ=255)
                                       tGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGAc                                                                                                      >  8:274274/1‑149 (MQ=255)
                                            gTGGCTCGGATCGCCTGGTTGAACTCCGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGAACGGGCTTTGTTTTTTTATATAGAGTGCCCCGCCTGACCGTTa                                                                                                 >  6:285224/1‑149 (MQ=255)
                                            gTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTa                                                                                                 <  2:302760/149‑1 (MQ=255)
                                             tGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTAc                                                                                                <  1:151757/149‑1 (MQ=255)
                                               gCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACa                                                                                               <  5:129338/148‑1 (MQ=255)
                                                     aTCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGccc                                                                                        <  7:274274/149‑1 (MQ=255)
                                                         ccTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATATAGAGTGCCAAGCCTGACCGTTACAGAAGCCCGca                                                                                     >  7:140103/1‑148 (MQ=255)
                                                                gAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGc                                                                                      >  8:249959/1‑140 (MQ=255)
                                                                gAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGc                                                                                      <  7:249959/140‑1 (MQ=255)
                                                                 aaCTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCAtgctg                                                                             <  3:355442/148‑1 (MQ=255)
                                                                   cTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCAtgctgc                                                                            >  2:225653/1‑147 (MQ=255)
                                                                   cTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCAtgctgc                                                                            <  1:225653/147‑1 (MQ=255)
                                                                      aGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgcgc                                                                       >  8:41330/1‑149 (MQ=255)
                                                                      aGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTT‑TCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgcgcg                                                                      >  6:19518/1‑149 (MQ=255)
                                                                       gCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgcgcg                                                                      >  3:236937/1‑149 (MQ=255)
                                                                             tgAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGaa                                                                 <  6:269097/148‑1 (MQ=255)
                                                                             tgAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGaa                                                                 >  2:132193/1‑148 (MQ=255)
                                                                                  tGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAg                                                              >  5:265863/1‑146 (MQ=255)
                                                                                       gCCATGGCCGAAAAGATCAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTAc                                                       <  3:245098/148‑1 (MQ=255)
                                                                                       gCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACa                                                      <  4:236937/149‑1 (MQ=255)
                                                                                                 aaaaGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAg                                            <  7:297738/149‑1 (MQ=255)
                                                                                                 aaaaGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAg                                            >  5:155871/1‑149 (MQ=255)
                                                                                                    aGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGcc                                         <  6:265863/149‑1 (MQ=255)
                                                                                                      aacaacTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCg                                                                                                    >  3:293895/1‑88 (MQ=255)
                                                                                                      aacaacTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCg                                                                                                    <  4:293895/88‑1 (MQ=255)
                                                                                                       acaacTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTg                                      <  1:132193/149‑1 (MQ=255)
                                                                                                          acTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTT‑TCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCg                                  <  5:19518/149‑1 (MQ=255)
                                                                                                                   acaACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgcc                          >  4:277337/1‑149 (MQ=255)
                                                                                                                   acaACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgcc                          >  5:99306/1‑149 (MQ=255)
                                                                                                                    caACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGACGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCgcct                         >  5:176446/1‑149 (MQ=255)
                                                                                                                         cTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACACGCCAGGCCCTGCTCGATGACGcctgcctg                    >  1:137666/1‑149 (MQ=255)
                                                                                                                             cTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGAt                                <  2:272603/133‑1 (MQ=255)
                                                                                                                             cTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGAt                                >  1:272603/1‑133 (MQ=255)
                                                                                                                                cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCACGATGCCGCCTGCCTGCTCATg              >  4:233002/1‑148 (MQ=255)
                                                                                                                                 ccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga            <  6:155871/149‑1 (MQ=255)
                                                                                                                                  cTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATg              >  7:113454/1‑146 (MQ=255)
                                                                                                                                  cTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATg              >  8:320151/1‑146 (MQ=255)
                                                                                                                                  aacaaagCCCGGTCAATGGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATg              >  7:146491/7‑146 (MQ=255)
                                                                                                                                     cAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGgaga          <  5:183269/147‑1 (MQ=255)
                                                                                                                                     cAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGt        <  3:277337/149‑1 (MQ=255)
                                                                                                                                       agCCCGGTCAATGGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTgg      >  7:43968/2‑149 (MQ=255)
                                                                                                                                        gCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTgg      >  4:220548/1‑148 (MQ=255)
                                                                                                                                        gCCCGGTCAATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTgg      <  3:220548/148‑1 (MQ=255)
                                                                                                                                           cGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGt  <  2:137666/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CGGCTCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGT  >  NC_002947/5525626‑5525913

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: ‑