Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,568,129 C→G 7.3% T229S (ACT→AGT)  PP_3150 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,568,1290CG 7.3% 252.7 / 6.5 79T229S (ACT→AGT) PP_3150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base G (2/4);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACCAGCAGAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCACCAAGCTG  >  NC_002947/3567982‑3568277
                                                                                                                                                    |                                                                                                                                                    
aCCAGCAGAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCAGCTGGCCTTGACt                                                                                                                                                     >  5:278837/1‑149 (MQ=255)
   agcagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTg                                                                                                                                                    >  5:21974/1‑147 (MQ=255)
   agcagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGaa                                                                                                                                                  <  6:73917/149‑1 (MQ=255)
    gcagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAg                                                                                                                                                 >  4:16478/1‑149 (MQ=255)
     cagAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAgg                                                                                                                                                >  6:237330/1‑149 (MQ=255)
       gAACCGCCAGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCAGTGGCTGTGCAAGCTTTCGCCCGCACCCGCAAGACCTTCCTGAATCACCCCTACGCCCGGCTGCTCTAGCCACAGCGAACCATCCTGATGGCCTTGACTGAAGGc                                                                                                                                               >  8:102533/1‑148 (MQ=255)
        aaCCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTCAGTCAAGGCCa                                                                                                                                             <  8:52063/149‑1 (MQ=255)
           cGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGtt                                                                                                                                          <  7:47072/149‑1 (MQ=255)
            gCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTCAGTCAAGGCCAgcag                                                                                                                                         <  8:25477/149‑4 (MQ=255)
             ccTGCTGCGCGTATTCGTCAGTGGCATGGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGAACCGGCAAGACCTTCCTGATTCTCCAGTTCGCCCGGCTGCTGGAGCCACAGGGCTGCTTGCTGCTGGCCTTGAGTGAAGGCCAGTTgc                                                                                                                                        <  7:102533/149‑1 (MQ=255)
                 ctgcGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTgcgc                                                                                                                                      <  2:349437/147‑1 (MQ=255)
                   gcgcGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCCTGCTGCTGGCCTTGACTGAAGGCCAGTTgcgcgcgc                                                                                                                                  >  5:73991/1‑149 (MQ=255)
                         aTTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAg                                                                                                                             <  1:32643/149‑1 (MQ=255)
                          ttCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGgc                                                                                                                           <  2:158697/149‑1 (MQ=255)
                          ttCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGgc                                                                                                                           >  1:356010/1‑149 (MQ=255)
                               cAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCggg                                                                                                                      >  6:81012/1‑149 (MQ=255)
                                aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCggg                                                                                                                      <  2:202282/148‑1 (MQ=255)
                                aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCggg                                                                                                                      <  6:53440/148‑1 (MQ=255)
                                aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCgg                                                                                                                       <  6:21974/147‑1 (MQ=255)
                                aGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGt                                                                                                                     >  8:141096/1‑149 (MQ=255)
                                  tGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTc                                                                                                                    >  8:156489/1‑148 (MQ=255)
                                    gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGaa                                                                                                                                                  >  2:251558/1‑116 (MQ=255)
                                    gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGaa                                                                                                                                                  <  1:251558/116‑1 (MQ=255)
                                    gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAg                                                                                                                 >  8:115235/1‑149 (MQ=255)
                                    gCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAg                                                                                                                 <  6:278837/149‑1 (MQ=255)
                                       tCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGAc                                                                                                              <  5:237330/149‑1 (MQ=255)
                                             aGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgcc                                                                                                        <  7:141096/149‑1 (MQ=255)
                                              gTCTGTGGCTGTGCAAGGTTTCGCCGGCACTGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgccg                                                                                                       >  3:318933/1‑149 (MQ=255)
                                               tCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgccgc                                                                                                      <  2:170998/149‑1 (MQ=255)
                                               tCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTgccgc                                                                                                      <  7:115235/149‑1 (MQ=255)
                                                        gtgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTgcgc                                                                                                                                      >  7:8485/1‑108 (MQ=255)
                                                        gtgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTgcgc                                                                                                                                      <  8:8485/108‑1 (MQ=255)
                                                        gtgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgc                                                                                             >  3:111794/1‑149 (MQ=255)
                                                         tgCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgcg                                                                                            >  3:21367/1‑149 (MQ=255)
                                                          gCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCgcgc                                                                                           >  5:53755/1‑149 (MQ=255)
                                                           cAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGGCGCTGCCGCCTATACCGCGCt                                                                                          >  2:199742/1‑149 (MQ=255)
                                                            aaGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCt                                                                                          >  1:156852/1‑148 (MQ=255)
                                                            aaGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTg                                                                                         >  6:238977/1‑149 (MQ=255)
                                                                tttCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCgg                                                                                                                       <  4:252740/115‑1 (MQ=255)
                                                                tttCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCgg                                                                                                                       >  3:252740/1‑115 (MQ=255)
                                                                tttCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCt                                                                                     >  2:43725/1‑149 (MQ=255)
                                                                   cGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTc                                                                                   <  1:43725/148‑1 (MQ=255)
                                                                    gCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCg                                                                                  <  6:231233/148‑1 (MQ=255)
                                                                      cggcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgcg                                                                                            >  8:281178/1‑136 (MQ=255)
                                                                      cggcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACcgcg                                                                                            <  7:281178/136‑1 (MQ=255)
                                                                       ggcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCa                                                                              >  8:314244/1‑149 (MQ=255)
                                                                        gcaccggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCCTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCaa                                                                             <  6:73991/149‑1 (MQ=255)
                                                                             ggcaAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGc                                                                        <  2:156852/149‑1 (MQ=255)
                                                                               caAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTg                                                                      <  2:356010/149‑1 (MQ=255)
                                                                                  gACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACg                                                                   <  4:318933/149‑1 (MQ=255)
                                                                                      ttCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAgct                                                               >  1:257700/1‑149 (MQ=255)
                                                                                          tGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGc                                                                                <  4:270380/128‑1 (MQ=255)
                                                                                          tGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGc                                                                                >  3:270380/1‑128 (MQ=255)
                                                                                              tCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCtata                                                                                                 <  8:249107/107‑1 (MQ=255)
                                                                                              tCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCtata                                                                                                 >  7:249107/1‑107 (MQ=255)
                                                                                                aCCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGt                                                     <  7:314244/149‑1 (MQ=255)
                                                                                                 ccAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGt                                                     >  3:95846/1‑148 (MQ=255)
                                                                                                  cAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  6:53755/149‑1 (MQ=255)
                                                                                                  cAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  3:248136/149‑1 (MQ=255)
                                                                                                  cAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  7:36927/149‑1 (MQ=255)
                                                                                                   aGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGAc                                                  <  5:81012/149‑1 (MQ=255)
                                                                                                        gCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   >  7:125777/1‑143 (MQ=255)
                                                                                                        gCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGa                                                   <  8:125777/143‑1 (MQ=255)
                                                                                                              ctgctgGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCaa                                       >  3:199953/1‑149 (MQ=255)
                                                                                                                          caGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGAtt                           <  7:156489/149‑1 (MQ=255)
                                                                                                                                aCCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGa                     <  3:16478/149‑1 (MQ=255)
                                                                                                                                  cTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATc                   <  2:257700/149‑1 (MQ=255)
                                                                                                                                   gcgcaactgGCCTTCAGTCAAGGCCAGCTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATcc                  >  8:256011/7‑149 (MQ=255)
                                                                                                                                       tgctgGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTAc              >  7:228587/1‑149 (MQ=255)
                                                                                                                                       aactgGCCTTCAGTCAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGATGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTAc              <  1:93259/147‑1 (MQ=255)
                                                                                                                                       aactgGCCTTCAGTCAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTAc              >  2:93259/3‑149 (MQ=255)
                                                                                                                                         ctgGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGGCGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGTAACGGCTGGCGATTGCGCGATCCGTAc              <  1:199742/147‑1 (MQ=255)
                                                                                                                                         ctgGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCg            <  4:21367/149‑1 (MQ=255)
                                                                                                                                         ctgGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCg            >  2:320362/1‑149 (MQ=255)
                                                                                                                                          tgGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCg            >  8:265346/1‑148 (MQ=255)
                                                                                                                                            gCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTa               <  8:185189/143‑1 (MQ=255)
                                                                                                                                            gCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTa               >  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)
                                                                                                                                                    |                                                                                                                                                    
ACCAGCAGAACCGCCTGCTGCGCGAATTCCTCAGTGGCATCGACGAGTCTGTGGCTGTGCAAGGTTTCGCCGGCACCGGCAAGACCTTCCTGATTCACCAGTTCGCCCGGCTGCTGGAGCCACAGCGCACCTTGCTGCTGGCCTTGACTGAAGGCCAGTTGCGCGCGCTGCAGGCGCGGGTCAAGGACGCTGCCGCCTATACCGCGCTGACCTTCGGCCAACTGGCTGACGAGCTGCTCAACCGTGACCTCACCAGCAACGGCTGGCGATTGCGCGATCCGTACCGCACCAAGCTG  >  NC_002947/3567982‑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: