Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,640,013 G→C 6.7% T313S (ACC→AGC)  cmoB ← tRNA 5‑methoxyuridine(34) synthase CmoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,640,0130GC 6.7% 190.4 / 3.9 60T313S (ACC→AGC) cmoBtRNA 5‑methoxyuridine(34) synthase CmoB
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

GCTGTTGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACGCAGCG  >  NC_002947/1639865‑1640152
                                                                                                                                                      |                                                                                                                                           
gCTGTTGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCAATCAggg                                                                                                                                             >  6:190796/1‑149 (MQ=255)
 cTGTTGCCGAGCGCGAGCGCCAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGCGt                                                                                                                                            >  3:265742/1‑149 (MQ=255)
  tGTTGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTg                                                                                                                                           >  3:276462/1‑149 (MQ=255)
    ttGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTTAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCGCCCTCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTg                                                                                                                                           <  1:84136/148‑1 (MQ=255)
    ttGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGgc                                                                                                                                         >  2:293869/1‑149 (MQ=255)
     tGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGgcg                                                                                                                                        >  7:286972/1‑149 (MQ=255)
       ccGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCACCCTGCTGgcgc                                                                                                                                       <  3:84781/148‑1 (MQ=255)
        cgagcgcgagcgCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCgcgg                                                                                                                                     >  2:191718/1‑149 (MQ=255)
        cgagcgcgagcgCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCgcgg                                                                                                                                     >  6:171348/1‑149 (MQ=255)
             gcgagcgcAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTgg                                                                                                                                <  6:116025/149‑1 (MQ=255)
                agcgcAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGggc                                                                                                                             <  1:191718/149‑1 (MQ=255)
                      aggcgcaagAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGgcaggcag                                                                                                                       >  1:290290/1‑149 (MQ=255)
                      aggcgcaagAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGgcaggcag                                                                                                                       <  7:165782/149‑1 (MQ=255)
                              gagaAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTGTTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGa                                                                                                               <  4:54533/149‑1 (MQ=255)
                                      gcgTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCCGCGTGGGGCAGGCAGGCCTTCGACGGTCTTg                                                                                                       >  3:73390/1‑149 (MQ=255)
                                        gTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGACGGGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGc                                                                                                      >  1:58817/1‑149 (MQ=255)
                                        gTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCt                                                                                                     <  1:293869/149‑1 (MQ=255)
                                         tGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCtg                                                                                                    <  2:290290/149‑1 (MQ=255)
                                          gAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCtgt                                                                                                   >  2:265572/1‑149 (MQ=255)
                                            ggTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTgg                                                                                                 <  1:265572/149‑1 (MQ=255)
                                                   aaGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCg                                                                                          >  1:178306/1‑149 (MQ=255)
                                                           aGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTg                                                                                                       >  7:123230/1‑128 (MQ=255)
                                                           aGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTg                                                                                                       <  8:123230/128‑1 (MQ=255)
                                                             gTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAgg                                                                                  >  8:68368/1‑147 (MQ=255)
                                                                  gACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGc                                                                           >  4:103589/1‑149 (MQ=255)
                                                                       cTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGt                                                                      <  4:265742/149‑1 (MQ=255)
                                                                       cTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTGTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGt                                                                      >  7:123419/1‑149 (MQ=255)
                                                                         gcgcCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGa                                                                    >  2:132045/1‑149 (MQ=255)
                                                                            ccgcGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCgg                                                                                         <  3:191386/125‑1 (MQ=255)
                                                                            ccgcGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCgg                                                                                         >  4:191386/1‑125 (MQ=255)
                                                                               cgAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTa                                                              <  2:58817/149‑1 (MQ=255)
                                                                                       tcctcG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTg                                                                 >  8:25709/1‑139 (MQ=255)
                                                                                       tcctcG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTg                                                                 <  7:25709/139‑1 (MQ=255)
                                                                                         ctcG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGTGGAAGTCGCTGAGTGACTGGTAGCGCATcc                                                      <  1:26330/148‑1 (MQ=255)
                                                                                         ctcG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATcc                                                      >  2:26330/1‑148 (MQ=255)
                                                                                            g‑ccgcAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAAt                                                   <  5:171348/148‑1 (MQ=255)
                                                                                                gcAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGt                                               >  2:143442/1‑149 (MQ=255)
                                                                                                   ggTGTATTCAGGCGCGCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGGCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGtgct                                            <  3:103589/149‑1 (MQ=255)
                                                                                                     tgtATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGtgctgc                                          >  1:166683/1‑149 (MQ=255)
                                                                                                         ttCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCgg                                                >  2:250941/1‑139 (MQ=255)
                                                                                                         ttCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCACGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCgg                                                <  1:250941/139‑1 (MQ=255)
                                                                                                                gcCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGctcctc                                >  4:153979/1‑148 (MQ=255)
                                                                                                                gcCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGctcctc                                >  1:89005/1‑148 (MQ=255)
                                                                                                                       cGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTg                         <  4:276462/148‑1 (MQ=255)
                                                                                                                       cGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTg                         <  7:68368/148‑1 (MQ=255)
                                                                                                                          gcgcCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCa                     <  2:178306/149‑1 (MQ=255)
                                                                                                                            gcCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTg                         <  2:277656/143‑1 (MQ=255)
                                                                                                                            gcCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTg                         >  1:277656/1‑143 (MQ=255)
                                                                                                                            gcCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACg                   <  3:153979/149‑1 (MQ=255)
                                                                                                                              cTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGtgctg                                           >  4:93161/1‑123 (MQ=255)
                                                                                                                              cTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGtgctg                                           <  3:93161/123‑1 (MQ=255)
                                                                                                                              cTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCt                 <  1:143442/149‑1 (MQ=255)
                                                                                                                                  tttATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCAc              <  5:190796/148‑1 (MQ=255)
                                                                                                                                      tttGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGa         <  1:132045/149‑1 (MQ=255)
                                                                                                                                      acgccgcgcCACCCTGCTGGCGC‑GCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTg                         >  6:143871/5‑132 (MQ=255)
                                                                                                                                      acgccgcgcCACCCTGCTGGCGC‑GCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTg                         <  5:143871/128‑1 (MQ=255)
                                                                                                                                        gccgcgcCACCCTGCTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACg       >  3:284739/3‑149 (MQ=255)
                                                                                                                                         gcgcgcCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACg       >  1:308359/1‑148 (MQ=255)
                                                                                                                                          cgcgcCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGGCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACg       >  6:270750/1‑147 (MQ=255)
                                                                                                                                           gcgcCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACGCAg    >  1:62998/1‑149 (MQ=255)
                                                                                                                                             gcCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACGCAgcg  <  5:270750/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
GCTGTTGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACGCAGCG  >  NC_002947/1639865‑1640152

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