Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,489,778 G→C 6.2% P538R (CCA→CGA)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,489,7780GC 6.2% 260.7 / 5.6 81P538R (CCA→CGA) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base G (35/40);  new base C (2/3);  total (38/43)
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.83e-01

TCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGC  >  NC_002947/3489631‑3489924
                                                                                                                                                   |                                                                                                                                                    
tCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTgg                                                                                                                                                     >  7:323743/1‑149 (MQ=255)
tCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTgg                                                                                                                                                     <  1:376837/149‑1 (MQ=255)
     gCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTgg                                                                                                                                                <  8:233497/148‑1 (MQ=255)
      cccAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTc                                                                                                                                              >  5:238411/1‑149 (MQ=255)
       ccAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCg                                                                                                                                             <  5:253038/149‑1 (MQ=255)
         aGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCAGGGTATCGATTTCGTCACGCAGCCTCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGcc                                                                                                                                           <  2:164962/149‑1 (MQ=255)
         aGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGcc                                                                                                                                           <  5:347038/149‑1 (MQ=255)
           cACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCtt                                                                                                                                         <  5:222548/149‑1 (MQ=255)
             cgcgCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGc                                                                                                                                       <  2:69766/149‑1 (MQ=255)
              gcgcTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCa                                                                                                                                      >  3:278476/1‑149 (MQ=255)
               cgcTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCAc                                                                                                                                     <  6:23718/149‑1 (MQ=255)
                gcTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGcc                                                                                                                                           >  5:328992/1‑142 (MQ=255)
                gcTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGcc                                                                                                                                           <  6:328992/142‑1 (MQ=255)
                      cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGcc                                                                                                                             <  8:200438/149‑1 (MQ=255)
                      cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGcc                                                                                                                             <  8:221727/149‑1 (MQ=255)
                      cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGT‑CAGGCACGGATGCCCATGGCGATGGCCCCCCAGTACCCCAGCGCCTGGACCAGTTG‑CTGGTAGCCTTTCACCCC‑CGcc                                                                                                                             >  8:114821/1‑148 (MQ=255)
                       aaGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACgcg‑cgc                                                                                                                              >  2:240720/1‑147 (MQ=255)
                       aaGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCCACTGGTCCACGCG‑CTGGTCga                                                                                                                                            <  3:138480/134‑2 (MQ=255)
                       aaGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCCACTGGTCCACGCG‑CTGGTCga                                                                                                                                            >  4:138480/1‑133 (MQ=255)
                           aGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTc                                                                                                                        >  3:146487/1‑149 (MQ=255)
                           aGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTc                                                                                                                        >  2:305476/1‑149 (MQ=255)
                            gggCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCg                                                                                                                       >  2:175966/1‑149 (MQ=255)
                             ggCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCgg                                                                                                                      >  2:56109/1‑149 (MQ=255)
                             ggCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCgg                                                                                                                      >  5:384722/1‑149 (MQ=255)
                               ccGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCa                                                                                                                    <  8:323743/149‑1 (MQ=255)
                                      gcgcTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCCACTGGTCCACGCG‑CTGGTCgaca                                                                                                                                          >  8:290793/1‑118 (MQ=38)
                                      gcgcTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCCACTGGTCCACGCG‑CTGGTCgaca                                                                                                                                          <  7:290793/121‑4 (MQ=38)
                                        gcTGCACGGTATCGATTTCGTCACGCAGCACCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCACTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTa                                                                                                           >  2:103732/1‑149 (MQ=255)
                                            cACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTg                                                                                                        <  5:314793/148‑1 (MQ=255)
                                            cACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTg                                                                                                        <  5:408178/148‑1 (MQ=255)
                                            cACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGAGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTg                                                                                                        <  1:305476/148‑1 (MQ=255)
                                                   tgcaTTTCGTCACGCCGCAGCTTGCCACGGGGTATGCCGGT‑CAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTgg                                                                                               <  7:114821/146‑1 (MQ=255)
                                                           gTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCgg                                                                                        >  7:236255/1‑149 (MQ=255)
                                                            tCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCCTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCgg                                                                                        <  1:343328/149‑1 (MQ=255)
                                                             cACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCa                                                                                      >  6:124566/1‑149 (MQ=255)
                                                                gcagcaTCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCg                                                                                   <  1:240720/149‑1 (MQ=255)
                                                                  agcaTCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGtt                                                                                 <  4:146487/149‑1 (MQ=255)
                                                                  agcaTCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGtt                                                                                 >  1:189195/1‑149 (MQ=255)
                                                                           cccAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTc                                                                        >  3:423449/1‑149 (MQ=255)
                                                                            ccAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTc                                                                        >  4:183336/1‑148 (MQ=255)
                                                                               gCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAgc                                                                    >  3:99833/1‑149 (MQ=255)
                                                                                cGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGGACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGACGCCGCCAACGATGCCGACGCCCAgca                                                                   <  5:10249/149‑2 (MQ=255)
                                                                                   gATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCt                                                                >  8:365493/1‑149 (MQ=255)
                                                                                    aTGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCtt                                                               <  4:278476/149‑1 (MQ=255)
                                                                                       ccGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAg                                                            <  6:384722/149‑1 (MQ=255)
                                                                                       ccGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAg                                                            <  7:352284/149‑1 (MQ=255)
                                                                                       ccGGTCCAGCCACTGATGACCTGGGCGATGGCGCCACTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAg                                                            <  1:103732/149‑1 (MQ=255)
                                                                                             cAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGtt                                                       <  2:74583/148‑1 (MQ=255)
                                                                                               gCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCGCCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCAt                                                                                     <  1:89840/117‑1 (MQ=255)
                                                                                               gCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCAt                                                                                     >  2:89840/1‑116 (MQ=255)
                                                                                               gCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGa                                                    <  8:236255/149‑1 (MQ=255)
                                                                                                ccACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGAc                                                   <  2:189195/149‑1 (MQ=255)
                                                                                                 cACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGAcc                                                  <  3:183336/149‑1 (MQ=255)
                                                                                                 cACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGAc                                                   >  8:302703/1‑148 (MQ=255)
                                                                                                    tgatgaCCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGt                                               >  4:193594/1‑149 (MQ=255)
                                                                                                      atgaCCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTc                                             >  4:82854/1‑149 (MQ=255)
                                                                                                        gaCCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGAc                                                   <  5:16862/141‑1 (MQ=255)
                                                                                                        gaCCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGAc                                                   >  6:16862/1‑141 (MQ=255)
                                                                                                           cTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGc                                         >  4:116016/1‑148 (MQ=255)
                                                                                                                cGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCctgc                                   >  8:70416/1‑149 (MQ=255)
                                                                                                                      cgccgcTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCAc                             <  1:175966/149‑1 (MQ=255)
                                                                                                                         cgcTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTgg                          >  2:61205/1‑149 (MQ=255)
                                                                                                                              tCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCAtt                     <  6:238411/149‑1 (MQ=255)
                                                                                                                               cGACCAGCGCGTGTACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTg                    <  6:36948/149‑1 (MQ=255)
                                                                                                                                 aCCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAg                  <  6:186738/149‑1 (MQ=255)
                                                                                                                                 aCCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAg                  <  6:193273/149‑1 (MQ=255)
                                                                                                                                 aCCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGACGCCGCCAACGATGCCGACGccc                                                                       >  2:140554/1‑96 (MQ=37)
                                                                                                                                 aCCAGCGCGTGGACCAGTCG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGACGCCGCCAACGATGCCGACGccc                                                                       <  1:140554/96‑1 (MQ=25)
                                                                                                                                    aGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctg               <  4:423449/149‑1 (MQ=255)
                                                                                                                                     gcgcGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctgc              >  5:408547/1‑149 (MQ=255)
                                                                                                                                     gcgcGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctgc              >  6:358720/1‑149 (MQ=255)
                                                                                                                                     gcgcGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctg               >  2:75642/1‑148 (MQ=255)
                                                                                                                                     gcgcGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctg               <  1:75642/148‑1 (MQ=255)
                                                                                                                                      cgcgTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctgct             <  7:70416/149‑1 (MQ=255)
                                                                                                                                      cgcgTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctgct             >  2:401851/1‑149 (MQ=255)
                                                                                                                                             aCCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctg      <  3:193594/149‑1 (MQ=255)
                                                                                                                                              ccAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCCGCGCCTTGAGCTGTTCGCCCAGTTCCAGCTCCTGCTGGCACGGGGCATTGAGCTGCTGTTCctgc     >  3:315944/1‑149 (MQ=255)
                                                                                                                                               cAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctgct    >  1:402818/1‑149 (MQ=255)
                                                                                                                                               cAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTGCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctgcg    >  7:283914/1‑148 (MQ=255)
                                                                                                                                                aGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctgctg   <  1:56109/149‑1 (MQ=255)
                                                                                                                                                 gTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctgctgc  >  2:226906/1‑149 (MQ=255)
                                                                                                                                                 gTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctgctgc  <  1:61205/149‑1 (MQ=255)
                                                                                                                                                 gTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTActgctgc  >  3:25085/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGG‑CTGGTCGCCTTGCACGCG‑CGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGC  >  NC_002947/3489631‑3489924

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