Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 495,966 A→T 6.2% L23Q (CTG→CAG)  PP_0408 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947495,9660AT 6.2% 244.6 / 5.0 81L23Q (CTG→CAG) PP_0408hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/46);  new base T (2/3);  total (32/49)
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.94e-01

CGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTT  >  NC_002947/495820‑496109
                                                                                                                                                  |                                                                                                                                                  
cGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAgg                                                                                                                                                 <  3:252903/149‑1 (MQ=255)
cGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAgg                                                                                                                                                 >  7:118842/1‑149 (MQ=255)
 ggCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAggg                                                                                                                                                <  5:107120/149‑1 (MQ=255)
    cTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGccgc                                                                                                                                            >  5:98906/1‑149 (MQ=255)
     tGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGcc                                                                                                                                           >  5:28389/1‑149 (MQ=255)
     tGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGcc                                                                                                                                           >  3:133988/1‑149 (MQ=255)
       cAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGcc                                                                                                                                           <  8:118842/147‑1 (MQ=255)
        aGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAg                                                                                                                                        <  5:253584/149‑1 (MQ=255)
        aGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAg                                                                                                                                        >  3:179023/1‑149 (MQ=255)
         gggTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGt                                                                                                                                       <  5:225688/149‑1 (MQ=255)
         gggTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGt                                                                                                                                       <  3:156991/149‑1 (MQ=255)
          ggTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTg                                                                                                                                      <  2:186836/149‑1 (MQ=255)
          ggTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGCCGatg                                                                                                                                       <  8:31283/148‑3 (MQ=255)
           gTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGCCGatgcc                                                                                                                                     <  6:98906/149‑5 (MQ=255)
            tCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTgg                                                                                                                                     >  5:53296/1‑148 (MQ=255)
                ctgctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGaa                                                                                                                                   <  3:352264/146‑1 (MQ=255)
                ctgctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGaa                                                                                                                                   <  1:1015/146‑1 (MQ=255)
                ctgctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGaa                                                                                                                                   <  6:248522/146‑1 (MQ=255)
                      aGCCTTGGGCAGCAGCCGTTCCCGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTTGAAGGACCGa                                                                                                                            <  1:326905/147‑1 (MQ=255)
                      aGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGa                                                                                                                            <  1:332125/147‑1 (MQ=255)
                      aGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAgg                                                                                                                          >  5:139076/1‑149 (MQ=255)
                      aGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAg                                                                                                                           <  6:205843/148‑1 (MQ=255)
                      aGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAg                                                                                                                           <  5:64293/148‑1 (MQ=255)
                       gCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGc                                                                                                                         <  8:128244/149‑1 (MQ=255)
                       gCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGc                                                                                                                         <  6:28389/149‑1 (MQ=255)
                            gggCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAg                                                                                                                           <  2:282466/142‑1 (MQ=255)
                            gggCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAg                                                                                                                           >  1:282466/1‑142 (MQ=255)
                            gggCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAgg                                                                                                                     <  5:257656/148‑1 (MQ=255)
                             ggCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGgc                                                                                                                    <  4:133988/148‑1 (MQ=255)
                                        ttCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGc                                                                                                          >  6:127933/1‑147 (MQ=255)
                                        ttCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGc                                                                                                          >  6:122985/1‑147 (MQ=255)
                                           cAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACAtt                                                                                                            >  6:33465/1‑142 (MQ=255)
                                           cAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACAtt                                                                                                            <  5:33465/142‑1 (MQ=255)
                                             gcgcCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGcgtc                                                                                                  <  4:100498/149‑1 (MQ=255)
                                               gcCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCTGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGCCCGAGGCCAGCCAGCACATTGCCCGcgtcgt                                                                                                >  4:230684/1‑149 (MQ=255)
                                                   gTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGa                                                                                            <  5:342381/149‑1 (MQ=255)
                                                    tGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGTCGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGa                                                                                            <  1:271240/149‑1 (MQ=255)
                                                     gCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGAATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGGCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGACAGCGTGGTGCGAt                                                                                           >  7:54050/1‑148 (MQ=255)
                                                      cATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATa                                                                                          >  3:353323/1‑148 (MQ=255)
                                                        tcctccTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATaaa                                                                                        <  3:257490/148‑1 (MQ=255)
                                                        tcctccTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATa                                                                                          >  4:257490/1‑146 (MQ=255)
                                                          ctcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgtg                                                                                     <  6:139076/149‑1 (MQ=255)
                                                             ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGtgtg                                                                                   <  3:46471/148‑1 (MQ=255)
                                                             ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGtgtg                                                                                   <  7:335037/148‑1 (MQ=255)
                                                                gCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgt                                                                                      <  6:371878/142‑1 (MQ=255)
                                                                gCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgt                                                                                      >  5:371878/1‑142 (MQ=255)
                                                                   gCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGa                                                                            <  2:320485/149‑1 (MQ=255)
                                                                         cAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAg                                                                                                                           <  4:182076/97‑1 (MQ=255)
                                                                         cAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAg                                                                                                                           >  3:182076/1‑97 (MQ=255)
                                                                          aGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAcc                                                                     >  8:48534/1‑149 (MQ=255)
                                                                           gCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAccc                                                                    <  4:179023/149‑1 (MQ=255)
                                                                                 gTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGc                                                               >  4:339982/1‑148 (MQ=255)
                                                                                  tGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCTGCATCGGCGGCCATGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGTACATGAGTCAACCCAAGGCCa                                                             <  3:230684/149‑1 (MQ=255)
                                                                                  tGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCa                                                             >  7:84657/1‑149 (MQ=255)
                                                                                     gCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGaa                                                          <  4:353323/149‑1 (MQ=255)
                                                                                             ccTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCg                                                  >  6:274440/1‑149 (MQ=255)
                                                                                                 ggCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCa                                               <  6:53296/148‑1 (MQ=255)
                                                                                                   cAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGa                                                                            <  5:19661/117‑1 (MQ=255)
                                                                                                   cAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGa                                                                            >  6:19661/1‑117 (MQ=255)
                                                                                                      ccGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTaa                                         <  3:339982/149‑1 (MQ=255)
                                                                                                        ggTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATa                                       <  5:127933/149‑1 (MQ=255)
                                                                                                        ggTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATa                                       <  8:36987/149‑1 (MQ=255)
                                                                                                         gTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTTGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAg                                      <  1:103840/149‑1 (MQ=255)
                                                                                                            ccATGGCATTTGCAGGTGGTTTTTCGGCGTCGGCGGCCAGGGCGGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTa                                    >  1:251669/1‑149 (MQ=255)
                                                                                                              aTGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACAGGCAGCCTAATAGTAttt                                 >  6:68215/1‑149 (MQ=255)
                                                                                                                 ccATTTGCAGGTGTTTTATCGGCATCGGCGGCCAGGGCGCCGAGTGTAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGc                               <  8:54050/147‑1 (MQ=255)
                                                                                                                  cATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGccc                             <  4:213742/149‑1 (MQ=255)
                                                                                                                   aTTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCg                            >  1:334934/1‑149 (MQ=255)
                                                                                                                        cAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAATGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGAc                        <  6:341655/148‑1 (MQ=255)
                                                                                                                        cAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACg                       <  2:111641/149‑1 (MQ=255)
                                                                                                                          ggtggtTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGATTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGc                     <  8:84657/149‑1 (MQ=255)
                                                                                                                            tggtTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGTCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGc                     >  1:289537/1‑148 (MQ=255)
                                                                                                                             ggtTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAt                  >  3:100913/1‑149 (MQ=255)
                                                                                                                             ggtTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAt                  >  1:185883/1‑149 (MQ=255)
                                                                                                                                     ccactCGGCGCCCTGGCCGCCGATGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAcc                                                                     >  4:31414/5‑89 (MQ=37)
                                                                                                                                     ccactCGGCGCCCTGGCCGCCGATGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAcc                                                                     <  3:31414/85‑1 (MQ=37)
                                                                                                                                      gCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATgg         >  7:38947/1‑149 (MQ=255)
                                                                                                                                        aTCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc       >  2:236028/1‑149 (MQ=255)
                                                                                                                                        aTCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc       >  7:206047/1‑149 (MQ=255)
                                                                                                                                           ggcggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTc    <  7:48534/149‑1 (MQ=255)
                                                                                                                                             cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac  <  2:334934/149‑3 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTT  >  NC_002947/495820‑496109

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: