Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,076,099 A→C 5.9% S155A (TCC→GCC)  PP_5328 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,076,0990AC 5.9% 214.4 / 4.0 68S155A (TCC→GCC) PP_5328phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base C (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ACGCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACACG  >  NC_002947/6075951‑6076246
                                                                                                                                                      |                                                                                                                                                   
acgCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGa                                                                                                                                                     <  2:206359/149‑1 (MQ=255)
 cgCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGAc                                                                                                                                                    <  4:311557/149‑1 (MQ=255)
  gCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACa                                                                                                                                                   <  7:291071/149‑1 (MQ=255)
   cTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGCGAAAACCAGCGCCTGGCCATTGGACa                                                                                                                                                   <  1:3938/149‑1 (MQ=255)
        gTGGTCCAGCGCACGCCCCTGCTAGTCGAGCACGAACGATTGATCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGAAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACATGccaa                                                                                                                                             >  8:83322/1‑147 (MQ=255)
              cAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacc                                                                                                                                        <  1:3993/148‑1 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                      <  8:216496/142‑1 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                      >  7:216496/1‑142 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgcg                                                                                                                               >  8:147950/1‑149 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgcg                                                                                                                               <  8:60108/149‑1 (MQ=255)
                        cccTGCTCGTCGAGCACGAACGATTGACTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgcgc                                                                                                                              >  1:4286/1‑149 (MQ=255)
                        cccTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCt                                                                                                                             >  2:280501/1‑149 (MQ=255)
                              tcgtcgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCggg                                                                                                                       >  4:285892/1‑149 (MQ=255)
                                gtcgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCt                                                                                                                     <  7:147950/149‑1 (MQ=255)
                                     gCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGt                                                                                                                  >  8:174443/1‑147 (MQ=255)
                                        cgaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTa                                                                                                             <  7:174443/149‑1 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGcgc                                                                                                                                              >  6:259203/1‑113 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGcgc                                                                                                                                              <  5:259203/115‑3 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTAc                                                                                                            >  2:11531/1‑149 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGATGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACAACCGGCGCGCTGCCCGGCTGGTCAGTAc                                                                                                            >  7:331220/1‑149 (MQ=255)
                                                      gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGt                                                                                                  <  2:289702/147‑1 (MQ=255)
                                                               gTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGgtgt                                                                                       >  4:286512/1‑149 (MQ=255)
                                                                   tCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGATCATATTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGCCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACTTGACCTGGGTGTCTGc                                                                                   >  3:58559/1‑149 (MQ=255)
                                                                    cGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                                    >  3:345519/1‑131 (MQ=255)
                                                                    cGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                                    <  4:345519/131‑1 (MQ=255)
                                                                     gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGa                                                                                                         >  7:217027/1‑125 (MQ=255)
                                                                     gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGa                                                                                                         <  8:217027/125‑1 (MQ=255)
                                                                        gCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCa                                                                              >  8:167149/1‑149 (MQ=255)
                                                                            gggCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACggg                                                                          <  5:37785/149‑1 (MQ=255)
                                                                            gggCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACggg                                                                          <  1:280501/149‑1 (MQ=255)
                                                                              gCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCCGGGTGTCTGCCGGCACGGGTa                                                                        <  1:11531/149‑1 (MQ=255)
                                                                                tGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGt                                                                      >  8:59484/1‑149 (MQ=255)
                                                                                  aTGGTTTTCTTGTTGTCCGTGTAGGTGATCTTGTAGGGGTGCTGGAAAATCAGCGACGGGCCATTGGACAGGCCCAGCACCACCAGCGCGCGGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGtt                                                                     <  8:331220/148‑1 (MQ=255)
                                                                                   tGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                                    >  2:144708/1‑148 (MQ=255)
                                                                                     gTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGt                                                                                                  <  5:47708/116‑1 (MQ=255)
                                                                                     gTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGt                                                                                                  >  6:47708/1‑116 (MQ=255)
                                                                                       tttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGacac                                                               >  8:301309/1‑149 (MQ=255)
                                                                                             gttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCtt                                                           >  6:246790/1‑147 (MQ=255)
                                                                                             gttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGa                                                         >  3:95989/1‑149 (MQ=255)
                                                                                              ttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                           <  7:337489/130‑1 (MQ=255)
                                                                                              ttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                           >  8:337489/1‑130 (MQ=255)
                                                                                                gtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGtt                                                      >  7:43139/1‑149 (MQ=255)
                                                                                                 tCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTc                                                     <  7:59484/149‑1 (MQ=255)
                                                                                                         ggTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCg                                              <  7:301309/148‑1 (MQ=255)
                                                                                                          gTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGca                                                                                                                                          >  7:109537/1‑55 (MQ=255)
                                                                                                          gTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGca                                                                                                                                          <  8:109537/55‑1 (MQ=255)
                                                                                                                  gTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCCGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTGTGCCGGGACGGGTAGTTCGACACGCGTGAGTTCGTTACCGGTCTTGGTag                                    >  4:41903/1‑148 (MQ=255)
                                                                                                                  gTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTg                                        <  2:4286/145‑1 (MQ=255)
                                                                                                                    aGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGATGGTCAGTACTGATGGACGTGACCTGGGTGTATGCAGGCACTGGTAGTTCGAAACGCTTGAGTTCGTTAACGGTCGTGATATc                                   >  8:331437/1‑148 (MQ=255)
                                                                                                                     ggTGTGCTGGAAAAGCAGCGCCTGGCCATTGGACAGGCCCAGCACCAGCAGCGCGCTGCCCGGGTGGTCAGTACTTATGGACGTGACCTGTGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCg                                  <  4:58559/148‑1 (MQ=255)
                                                                                                                     ggTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCg                                  <  1:144708/148‑1 (MQ=255)
                                                                                                                      gtgtGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaa                                <  8:43139/149‑1 (MQ=255)
                                                                                                                            tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACaa                          >  4:301755/1‑149 (MQ=255)
                                                                                                                             ggAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                                    >  3:192221/1‑106 (MQ=255)
                                                                                                                             ggAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                                    <  4:192221/106‑1 (MQ=255)
                                                                                                                               gAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATg                        >  8:202713/1‑149 (MQ=255)
                                                                                                                                aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAcc                                               <  3:59081/125‑1 (MQ=255)
                                                                                                                                aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAcc                                               >  4:59081/1‑125 (MQ=255)
                                                                                                                                aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGc                       <  5:246790/149‑1 (MQ=255)
                                                                                                                                  aaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCTGTCTTGGTATCGAAGAACAATGCCt                     <  3:286512/149‑1 (MQ=255)
                                                                                                                                    ccAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaag                               >  3:105879/1‑137 (MQ=255)
                                                                                                                                    ccAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaag                               <  4:105879/137‑1 (MQ=255)
                                                                                                                                     cAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGc                  >  5:301412/1‑149 (MQ=255)
                                                                                                                                     cAGCGCCTGGACATTTGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGc                  <  7:331437/149‑1 (MQ=255)
                                                                                                                                       gggcCTGTCCAATGGCCAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTa                                     >  4:315014/3‑128 (MQ=255)
                                                                                                                                       gggcCTGTCCAATGGCCAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTa                                     <  3:315014/126‑1 (MQ=255)
                                                                                                                                           cTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGt            <  7:167149/149‑1 (MQ=255)
                                                                                                                                                    ggACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAaca    <  4:95989/148‑1 (MQ=255)
                                                                                                                                                     gACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACACg  >  6:136580/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ACGCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACACG  >  NC_002947/6075951‑6076246

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: