Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,568,131 G→C 6.4% E230Q (GAA→CAA)  PP_3150 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,568,1310GC 6.4% 251.4 / 6.0 78E230Q (GAA→CAA) PP_3150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base C (2/3);  total (35/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 = 6.95e-01

AGCAGAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCACCAAGCTG  >  NC_002947/3567985‑3568277
                                                                                                                                                   |                                                                                                                                                  
agcagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTg                                                                                                                                                    >  5:21974/1‑147 (MQ=255)
agcagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGaa                                                                                                                                                  <  6:73917/149‑1 (MQ=255)
 gcagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAg                                                                                                                                                 >  4:16478/1‑149 (MQ=255)
  cagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAgg                                                                                                                                                >  6:237330/1‑149 (MQ=255)
    gAACCGCCAGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCAGTGGCTGTGCAAGCTTTCGCCCGCACCCGCAAGACCTTCCTGAATCACCCCTACGCCCGGCTGCTCTAGCCACAGCGAACCATCCTGATGG‑CCTTGACTGAAGGc                                                                                                                                               >  8:102533/1‑148 (MQ=255)
     aaCCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTCAGTCAAGGCCa                                                                                                                                             <  8:52063/149‑1 (MQ=255)
        cGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGtt                                                                                                                                          <  7:47072/149‑1 (MQ=255)
         gCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTCAGTCAAGGCCAgcag                                                                                                                                         <  8:25477/149‑4 (MQ=255)
          ccTGCTGCGCGTATTCGTCAGTGGCATGGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGAACCGGCAAGACCTTCCTGATTCTCCAGTTCGCCCGGCTGCTGGAGCCACAGGGCTGCTTGCTGCTGG‑CCTTGAGTGAAGGCCAGTTgc                                                                                                                                        <  7:102533/149‑1 (MQ=255)
              ctgcGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTgcgc                                                                                                                                      <  2:349437/147‑1 (MQ=255)
                gcgcGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCCTGCTGCTGG‑CCTTGACTGAAGGCCAGTTgcgcgcgc                                                                                                                                  >  5:73991/1‑149 (MQ=255)
                      aTTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAg                                                                                                                             <  1:32643/149‑1 (MQ=255)
                       ttCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGgc                                                                                                                           >  1:356010/1‑149 (MQ=255)
                       ttCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGgc                                                                                                                           <  2:158697/149‑1 (MQ=255)
                            cAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCggg                                                                                                                      >  6:81012/1‑149 (MQ=255)
                             aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCggg                                                                                                                      <  6:53440/148‑1 (MQ=255)
                             aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCggg                                                                                                                      <  2:202282/148‑1 (MQ=255)
                             aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCgg                                                                                                                       <  6:21974/147‑1 (MQ=255)
                             aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGt                                                                                                                     >  8:141096/1‑149 (MQ=255)
                               tGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTc                                                                                                                    >  8:156489/1‑148 (MQ=255)
                                 gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGaa                                                                                                                                                  <  1:251558/116‑1 (MQ=255)
                                 gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGaa                                                                                                                                                  >  2:251558/1‑116 (MQ=255)
                                 gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAg                                                                                                                 <  6:278837/149‑1 (MQ=255)
                                 gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAg                                                                                                                 >  8:115235/1‑149 (MQ=255)
                                    tCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGAc                                                                                                              <  5:237330/149‑1 (MQ=255)
                                          aGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgcc                                                                                                        <  7:141096/149‑1 (MQ=255)
                                           gTCTGTGGCTGTGCAAGGTTTCGCCGGCACTGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgccg                                                                                                       >  3:318933/1‑149 (MQ=255)
                                            tCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgccgc                                                                                                      <  2:170998/149‑1 (MQ=255)
                                            tCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgccgc                                                                                                      <  7:115235/149‑1 (MQ=255)
                                                     gtgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTgcgc                                                                                                                                      >  7:8485/1‑108 (MQ=255)
                                                     gtgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTgcgc                                                                                                                                      <  8:8485/108‑1 (MQ=255)
                                                     gtgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgc                                                                                             >  3:111794/1‑149 (MQ=255)
                                                      tgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgcg                                                                                            >  3:21367/1‑149 (MQ=255)
                                                       gCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCgcgc                                                                                           >  5:53755/1‑149 (MQ=255)
                                                        cAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGGCGCTGCCGCCTATACCGCGCt                                                                                          >  2:199742/1‑149 (MQ=255)
                                                         aaGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCt                                                                                          >  1:156852/1‑148 (MQ=255)
                                                         aaGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTg                                                                                         >  6:238977/1‑149 (MQ=255)
                                                             tttCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCgg                                                                                                                       <  4:252740/115‑1 (MQ=255)
                                                             tttCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCgg                                                                                                                       >  3:252740/1‑115 (MQ=255)
                                                             tttCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCt                                                                                     >  2:43725/1‑149 (MQ=255)
                                                                cGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTc                                                                                   <  1:43725/148‑1 (MQ=255)
                                                                 gCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCg                                                                                  <  6:231233/148‑1 (MQ=255)
                                                                   cggcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgcg                                                                                            <  7:281178/136‑1 (MQ=255)
                                                                   cggcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgcg                                                                                            >  8:281178/1‑136 (MQ=255)
                                                                    ggcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCa                                                                              >  8:314244/1‑149 (MQ=255)
                                                                     gcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCCTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCaa                                                                             <  6:73991/149‑1 (MQ=255)
                                                                          ggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGc                                                                        <  2:156852/149‑1 (MQ=255)
                                                                            caAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTg                                                                      <  2:356010/149‑1 (MQ=255)
                                                                               gACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACg                                                                   <  4:318933/149‑1 (MQ=255)
                                                                                   ttCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAgct                                                               >  1:257700/1‑149 (MQ=255)
                                                                                       tGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGc                                                                                >  3:270380/1‑128 (MQ=255)
                                                                                       tGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGc                                                                                <  4:270380/128‑1 (MQ=255)
                                                                                           tCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCtata                                                                                                 >  7:249107/1‑107 (MQ=255)
                                                                                           tCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCtata                                                                                                 <  8:249107/107‑1 (MQ=255)
                                                                                             aCCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGt                                                     <  7:314244/149‑1 (MQ=255)
                                                                                              ccAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGt                                                     >  3:95846/1‑148 (MQ=255)
                                                                                               cAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  6:53755/149‑1 (MQ=255)
                                                                                               cAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  7:36927/149‑1 (MQ=255)
                                                                                               cAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  3:248136/149‑1 (MQ=255)
                                                                                                aGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGAc                                                  <  5:81012/149‑1 (MQ=255)
                                                                                                     gCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   >  7:125777/1‑143 (MQ=255)
                                                                                                     gCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  8:125777/143‑1 (MQ=255)
                                                                                                           ctgctgGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCaa                                       >  3:199953/1‑149 (MQ=255)
                                                                                                                       caGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGAtt                           <  7:156489/149‑1 (MQ=255)
                                                                                                                             aCCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGa                     <  3:16478/149‑1 (MQ=255)
                                                                                                                               cTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATc                   <  2:257700/149‑1 (MQ=255)
                                                                                                                                gcgcaactgG‑CCTTCAGTCAAGGCCAGCTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATcc                  >  8:256011/7‑149 (MQ=255)
                                                                                                                                    tgctgG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTAc              >  7:228587/1‑149 (MQ=255)
                                                                                                                                    aactgG‑CCTTCAGTCAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGATGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTAc              <  1:93259/147‑1 (MQ=255)
                                                                                                                                    aactgG‑CCTTCAGTCAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTAc              >  2:93259/3‑149 (MQ=255)
                                                                                                                                      ctgG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGGCGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGTAACGGCTGGCGATTGCGCGATCCGTAc              <  1:199742/147‑1 (MQ=255)
                                                                                                                                      ctgG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCg            <  4:21367/149‑1 (MQ=255)
                                                                                                                                      ctgG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCg            >  2:320362/1‑149 (MQ=255)
                                                                                                                                       tgG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCg            >  8:265346/1‑148 (MQ=255)
                                                                                                                                         g‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTa               <  8:185189/143‑1 (MQ=255)
                                                                                                                                         g‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTa               >  7:185189/1‑143 (MQ=255)
                                                                                                                                            cTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCAcc        >  6:90899/1‑148 (MQ=255)
                                                                                                                                               gACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCACCaa      <  5:90899/147‑1 (MQ=255)
                                                                                                                                                 cTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCACCAAGCTg  <  4:111794/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCAGAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGG‑CCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCACCAAGCTG  >  NC_002947/3567985‑3568277

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