Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,640,011 G→T 6.9% L314M (CTG→ATG) ‡ 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,0110GT 6.9% 185.6 / 3.9 60L314M (CTG→ATG) ‡cmoBtRNA 5‑methoxyuridine(34) synthase CmoB
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base T (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 = 3.03e-01

GGCTGTTGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACGCAGCG  >  NC_002947/1639864‑1640152
                                                                                                                                                     |                                                                                                                                             
ccgttctGCTGAGCTCGAGCGCAAGGCTCAAGTGAAGGCGAGTGAGGTTGCGTATGATTGAGGTATCGACACCTGCTCCGCGAC‑GGCTCGTCG‑CCGCAGTTGTATTGAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCcagcag                                                                                                                                               <  6:305947/142‑1 (MQ=255)
 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)
        ccGAGCGCGAGCGCAAGGCGCAAGAGAAGGAGCGTGAGGTGGACAAGGATTGAGGTATCGACACATGCGCCGCGAC‑GGCTCCACG‑CCGCAGGTGTCTTAATGCGCACCTGCGGGCGCCTTTTTTATTTGCTCCCCAGCAggagatctg                                                                                                                                       >  5:305947/1‑141 (MQ=255)
         cgagcgcgagcgCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCgcgg                                                                                                                                     >  6:171348/1‑149 (MQ=255)
         cgagcgcgagcgCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCgcgg                                                                                                                                     >  2:191718/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)
                                                                                                                                                     |                                                                                                                                             
GGCTGTTGCCGAGCGCGAGCGCAAGGCGCAAGAGAAGGCGCGTGAGGTGGCCAAGGATTGAGGTATCGACACCTGCGCCGCGAC‑GGCTCCTCG‑CCGCAGGTGTATTCAGGCGCCCCTGCGGGCGCCTTTTTTATTTGCGCGCCAGCAGGGTGGCGCGGCGTGGGGCAGGCAGGCCTTCGACGGTCTTGCTGTGGTCGTTCGGGTCGAGGAAGTCGCTGAGTGACTGGTAGCGCATCCAATCGGTGCTGCGCTGCTCCTCGACGCTGGTCACGCTCACGTCGACGCAGCG  >  NC_002947/1639864‑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: ‑