Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,088,661 T→G 6.2% intergenic (+157/‑97) pmbA → / → PP_0943 regulatory protease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,088,6610TG 6.2% 205.3 / 4.2 65intergenic (+157/‑97)pmbA/PP_0943regulatory protease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGC  >  NC_002947/1088516‑1088802
                                                                                                                                                  |                                                                                                                                             
cgcCGGCCCGTTCGCGGGTAAACGCGCTCCTACAGGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTc                                                                                                                                              <  3:80513/147‑1 (MQ=255)
 gcCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATGCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTcgc                                                                                                                                            >  1:99916/1‑149 (MQ=255)
   cGGCCCTTTCGCGTGTAAACCCGCTCCTACACGGTGGTGGACGTACTTCGGGGTTGGACACTGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGcc                                                                                                                                         <  5:59059/149‑1 (MQ=255)
      cccTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCggg                                                                                                                                      >  5:393427/1‑149 (MQ=255)
         tttCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCGCGAGGGCCGGGtct                                                                                                                                   >  1:322548/1‑147 (MQ=255)
         tttCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGttt                                                                                                                                   >  6:228276/1‑149 (MQ=255)
            cgcgGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTa                                                                                                                                >  4:152576/1‑149 (MQ=255)
            cgcgGGTAAACCCGCTCCTACACGGTGGTGGACGTAACTCGGGGATTGACACATGCAGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCGCGAGGGCCGGGtcttc                                                                                                                                 <  5:11803/148‑2 (MQ=37)
             gcgGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTAt                                                                                                                               >  8:39560/1‑149 (MQ=255)
              cgGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCGCGAGGGCCGGGtcttcgga                                                                                                                              <  2:322548/149‑5 (MQ=37)
                    aaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGGTGGACACAGGCCGGTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTtcatc                                                                                                                         <  5:195927/148‑1 (MQ=255)
                       ccGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCt                                                                                                                     >  6:275939/1‑149 (MQ=255)
                                 acGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCTGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGt                                                                                                                                     <  8:327/123‑1 (MQ=255)
                                 acGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCTGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGt                                                                                                                                     >  7:327/1‑123 (MQ=255)
                                       gtggACGTACATCGGTTTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGAtt                                                                                                     <  1:158629/149‑1 (MQ=255)
                                       gtggACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGAt                                                                                                      <  5:228276/148‑1 (MQ=255)
                                          gACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCt                                                                                                   <  6:146043/148‑1 (MQ=255)
                                           aCGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTAt                                                                                                 <  5:301953/149‑1 (MQ=255)
                                             gTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTAttt                                                                                               <  6:188703/149‑1 (MQ=255)
                                             gTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTAtt                                                                                                >  3:32314/1‑148 (MQ=255)
                                                cATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATc                                                                                            <  3:152576/149‑1 (MQ=255)
                                                  tCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGAtt                                                                                                     <  3:75550/138‑1 (MQ=255)
                                                  tCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGAtt                                                                                                     >  4:75550/1‑138 (MQ=255)
                                                  tCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCAc                                                                                          >  2:253959/1‑149 (MQ=255)
                                                      gggTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACtt                                                                                        >  3:200547/1‑148 (MQ=255)
                                                           ggACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaataa                                                                                   <  4:279452/148‑1 (MQ=255)
                                                           ggACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTaataat                                                                                  >  1:69783/1‑149 (MQ=255)
                                                            gACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTaataat                                                                                  >  1:214024/1‑148 (MQ=255)
                                                               acaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAat                                                                              <  4:27766/149‑1 (MQ=255)
                                                                  ggCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAatat                                                                            >  5:402946/1‑148 (MQ=255)
                                                                   gCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCa                                                                          >  8:325586/1‑149 (MQ=255)
                                                                     cGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCAtt                                                                        <  6:393427/149‑1 (MQ=255)
                                                                         gTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATc                                                                     <  2:99916/148‑1 (MQ=255)
                                                                         gTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATc                                                                     <  7:39560/148‑1 (MQ=255)
                                                                            ggAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAg                                                                  <  7:325586/148‑1 (MQ=255)
                                                                             gAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGt                                                                 >  5:109/1‑148 (MQ=255)
                                                                                   gTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATg                                                                                 <  2:294702/126‑1 (MQ=255)
                                                                                   gTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAAACCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATg                                                                                 >  1:294702/1‑126 (MQ=255)
                                                                                       aCCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCga                                                       <  6:402946/148‑1 (MQ=255)
                                                                                            cgAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGt                                                 <  1:253959/149‑1 (MQ=255)
                                                                                                   gCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGa                                           >  8:211626/1‑148 (MQ=255)
                                                                                                     cgtaaCATCACTACACCTTTCCCTGCAAATCCGAAGACCCTGCCCTCGCTCCTGTTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAt                                        <  5:260/146‑1 (MQ=255)
                                                                                                     cGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAt                                        >  7:342353/1‑149 (MQ=255)
                                                                                                        aaCAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCg                                     >  3:136175/1‑149 (MQ=255)
                                                                                                        aaCAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCg                                     <  6:109/149‑1 (MQ=255)
                                                                                                               cTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcg                              >  5:168475/1‑149 (MQ=255)
                                                                                                                   acCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGaa                          <  5:275939/149‑1 (MQ=255)
                                                                                                                   acCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGaa                          >  4:20029/1‑149 (MQ=255)
                                                                                                                     cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcg                              <  6:125137/143‑1 (MQ=255)
                                                                                                                     cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcg                              >  5:125137/1‑143 (MQ=255)
                                                                                                                        tCCCTGAAAATACGAAGAACCGGCCCTCGCGGCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGc                      <  8:342353/148‑1 (MQ=255)
                                                                                                                          ccTGCAAATCCGAAGACCCGGCCCTCGCGCCGTGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGccc                    <  2:69783/148‑1 (MQ=255)
                                                                                                                           cTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAGTAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGc                      <  4:150773/145‑1 (MQ=255)
                                                                                                                           cTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAGTAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGc                      >  3:150773/1‑145 (MQ=255)
                                                                                                                           cTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTa                  >  5:93755/1‑149 (MQ=255)
                                                                                                                             gCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTAcc                >  6:272563/1‑149 (MQ=255)
                                                                                                                                 aTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTgg             >  5:55390/1‑148 (MQ=255)
                                                                                                                                   ataaaaaCCCGGCGCGAGGGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGaa           >  7:81493/7‑148 (MQ=255)
                                                                                                                                     gAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAAc        >  1:18119/1‑149 (MQ=255)
                                                                                                                                     gAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATACGTCCTCCGCGAACTGCCCTACCTGGAAAAc        >  6:52695/1‑149 (MQ=255)
                                                                                                                                      aaGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGCCGATGATGTTCATGAAATCCGTCCTCAGCGAACTGCCCTACATGGAAAACt       >  3:145736/1‑149 (MQ=255)
                                                                                                                                        gACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTgg     <  6:93755/149‑1 (MQ=255)
                                                                                                                                           ccGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGgcgc  <  5:272563/149‑1 (MQ=255)
                                                                                                                                              gCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGaa                                                                               >  2:107172/1‑69 (MQ=255)
                                                                                                                                              gCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGaa                                                                               <  1:107172/69‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGC  >  NC_002947/1088516‑1088802

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: