Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,871,421 A→G 6.8% Q232R (CAG→CGG)  cobC → adenosylcobinamide‑phosphate synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,871,4210AG 6.8% 169.4 / 2.4 59Q232R (CAG→CGG) cobCadenosylcobinamide‑phosphate synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base G (2/2);  total (28/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 = 9.15e-01

GGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGAGCAT  >  NC_002947/1871273‑1871566
                                                                                                                                                    |                                                                                                                                                 
gggCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCa                                                                                                                                                   <  7:25028/149‑1 (MQ=255)
   cTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGc                                                                                                                                                <  7:138500/149‑1 (MQ=255)
     gATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGc                                                                                                                                                >  3:128988/1‑147 (MQ=255)
      aTCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGcccg                                                                                                                                             >  8:306447/1‑147 (MQ=255)
      aTCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCa                                                                                                                                              <  3:220418/148‑1 (MQ=255)
            ctgcGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCTGGGCCAGcacccg                                                                                                                                       <  5:228152/149‑6 (MQ=255)
            ctgcGCTCGTTTGGCAAGATCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTc                                                                                                                                       <  3:269628/149‑1 (MQ=255)
               cgcgcgTTTGGCAAGTTCTTCCGCTGGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCt                                                                                                                                      <  8:150892/147‑1 (MQ=255)
                gcTCGTTCGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCGTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCtgat                                                                                                                                   <  4:48666/149‑1 (MQ=255)
                    gTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATgaga                                                                                                                               <  4:139570/149‑1 (MQ=255)
                        ggCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTg                                                                                                                            <  3:246124/148‑1 (MQ=255)
                          cAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTg                                                                                                                            >  7:112904/1‑146 (MQ=255)
                          cAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTg                                                                                                                            <  8:112904/146‑1 (MQ=255)
                           aaGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCa                                                                                                                        >  3:195803/1‑149 (MQ=255)
                              ttcttcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAGGACAGCTTGCAGCTGCGCCTGGCCGAGCTGCTCGGCCCCTTGCCCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgcc                                                                                                                     >  1:245747/1‑148 (MQ=255)
                                cttcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgcgcg                                                                                                                   >  2:73348/1‑149 (MQ=255)
                                   cGGCTTGGCCGGCGTGCGCCTGTGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgcgcg                                                                                                                   <  4:43124/146‑1 (MQ=255)
                                             ggCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTCGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACg                                                                                                      <  2:245747/149‑1 (MQ=255)
                                              gCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGc                                                                                                     >  3:293170/1‑149 (MQ=255)
                                                         ggCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGtgctgctgc                                                                                          >  5:313750/1‑149 (MQ=255)
                                                          gCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGtgctgctgc                                                                                          >  2:264704/1‑148 (MQ=255)
                                                           cTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGtgctgctgc                                                                                          >  6:197643/1‑147 (MQ=255)
                                                                    cgccgAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGGCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGc                                                                               >  2:252414/1‑149 (MQ=255)
                                                                     gccgAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCt                                                                              >  5:129476/1‑149 (MQ=255)
                                                                       cgAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCGGCGGCTTg                                                                            <  8:272588/149‑1 (MQ=255)
                                                                            acaGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCAc                                                                       <  6:313750/149‑1 (MQ=255)
                                                                                cTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACgtgt                                                                   >  3:24566/1‑149 (MQ=255)
                                                                                cTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACgtgt                                                                   <  4:24566/149‑1 (MQ=255)
                                                                                    ctgctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTgcg                                                               >  3:52699/1‑149 (MQ=255)
                                                                                    ctgctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTgcg                                                               <  1:264704/149‑1 (MQ=255)
                                                                                       ctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCa                                                             <  8:128896/148‑1 (MQ=255)
                                                                                                  cGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTgg                                                 >  7:230888/1‑149 (MQ=255)
                                                                                                         ctCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCCGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAg                                          >  5:118036/1‑149 (MQ=255)
                                                                                                         ctCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTa                                           <  1:73348/148‑1 (MQ=255)
                                                                                                          tCGGGCCGTTGACGGTCCAGGGGCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAgt                                         >  7:284010/1‑149 (MQ=255)
                                                                                                            ggCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAgtgg                                       >  2:26326/1‑149 (MQ=255)
                                                                                                            ggCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAgtg                                        <  8:54490/148‑1 (MQ=255)
                                                                                                            ggCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAgt                                         <  5:197643/147‑1 (MQ=255)
                                                                                                                  tGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTg                                 >  2:25229/1‑149 (MQ=255)
                                                                                                                      ccGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCt                                            >  4:42859/1‑134 (MQ=255)
                                                                                                                      ccGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGAc                             <  6:129476/149‑1 (MQ=255)
                                                                                                                        gTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCt                           >  5:203/1‑149 (MQ=255)
                                                                                                                        gTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCt                           <  4:195803/149‑1 (MQ=255)
                                                                                                                         tCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGgc                                                >  1:344911/1‑127 (MQ=255)
                                                                                                                         tCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGgc                                                <  2:344911/127‑1 (MQ=255)
                                                                                                                           aaCGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGt                         <  5:269267/148‑1 (MQ=255)
                                                                                                                            aCGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTc                       <  1:26326/149‑1 (MQ=255)
                                                                                                                              ggCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCa                     <  1:252414/149‑1 (MQ=255)
                                                                                                                              ggCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCa                     <  4:293170/149‑1 (MQ=255)
                                                                                                                                    aCCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGt               <  1:25229/149‑1 (MQ=255)
                                                                                                                                          aggCTGGCCTGGGCCAGCCTGTCTGATGAGACTGCCAAGCTCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAg         >  2:54027/3‑149 (MQ=255)
                                                                                                                                           tGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGc        <  8:230888/149‑1 (MQ=255)
                                                                                                                                            gCTGGCCTGGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGc                                <  7:175587/124‑1 (MQ=255)
                                                                                                                                            gCTGGCCTGGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGc                                >  8:175587/1‑124 (MQ=255)
                                                                                                                                              tGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGAg     >  6:49234/1‑149 (MQ=255)
                                                                                                                                              tGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCACAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGc        >  3:153843/1‑146 (MQ=255)
                                                                                                                                                gCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGAGCa   >  4:191508/1‑149 (MQ=255)
                                                                                                                                                gCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGAGCa   >  6:273434/1‑149 (MQ=255)
                                                                                                                                                gCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCACAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGa      <  4:153843/146‑1 (MQ=255)
                                                                                                                                                 cccAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGAGCAt  <  4:52699/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGCGCCTAGTGGTGGCTGCGACCTGTTCCAGTACGTGCGTAGCGAGCAT  >  NC_002947/1871273‑1871566

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