Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,923,896 C→G 5.7% E306D (GAG→GAC peaB ← quinohemoprotein amine dehydrogenase modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,923,8960CG 5.7% 233.2 / 3.2 70E306D (GAG→GACpeaBquinohemoprotein amine dehydrogenase modification protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base G (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

CACACGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGCGCCCCCAC  >  NC_002947/3923754‑3924042
                                                                                                                                              |                                                                                                                                                  
cacaCGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCg                                                                                                                                              >  2:27902/1‑149 (MQ=255)
 acacGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGc                                                                                                                                             >  2:280040/1‑149 (MQ=255)
  cacGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGAGCTGGTGCAGGTCTTCGcc                                                                                                                                            <  4:250698/149‑1 (MQ=255)
   acGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCa                                                                                                                                           <  4:15659/149‑1 (MQ=255)
   acGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCa                                                                                                                                           <  1:97284/149‑1 (MQ=255)
     ggCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGt                                                                                                                                         <  7:56812/149‑1 (MQ=255)
      gCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTg                                                                                                                                        <  7:11868/149‑1 (MQ=255)
            acttgttGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCGAGTGCCGCTTCCAGTTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTg                                                                                                                                  <  7:214031/148‑1 (MQ=255)
                     ggCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCATTGCCGCTTCCACGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGc                                                                                                                                                    <  1:245612/122‑1 (MQ=255)
                     ggCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGc                                                                                                                                                    >  2:245612/1‑122 (MQ=255)
                      gCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGcc                                                                                                                         <  3:33690/148‑1 (MQ=255)
                        cTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAgg                                                                                                                                     >  8:115087/1‑134 (MQ=255)
                        cTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATg                                                                                                                      >  4:267597/1‑149 (MQ=255)
                        cTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATg                                                                                                                      >  2:123711/1‑149 (MQ=255)
                            gTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCGTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGACGc                                                                                                                  >  3:146879/1‑149 (MQ=255)
                                 tGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCagag                                                                                                             <  1:27902/149‑1 (MQ=255)
                                         aTCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCgg                                                                                                     <  3:142978/149‑1 (MQ=255)
                                          tCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCggg                                                                                                    <  5:300077/149‑1 (MQ=255)
                                            aaCTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCgggg                                                                                                   <  3:197592/148‑1 (MQ=255)
                                            aaCTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGc                                                                                                  <  4:120205/149‑1 (MQ=255)
                                                 gTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCagag                                                                                                             <  4:72052/133‑1 (MQ=255)
                                                 gTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCagag                                                                                                             >  3:72052/1‑133 (MQ=255)
                                                 gTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGc                                                                                             >  1:53453/1‑149 (MQ=255)
                                                 gTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGAGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGc                                                                                             >  3:133958/1‑149 (MQ=255)
                                                       gTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCt                                                                                       >  5:84296/1‑149 (MQ=255)
                                                          ggAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCGGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCgg                                                                                    >  7:218809/1‑149 (MQ=255)
                                                            agaAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGc                                                                                   >  6:175410/1‑148 (MQ=255)
                                                               aGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGCCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGaaa                                                                               <  1:240824/149‑1 (MQ=255)
                                                                gCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAAc                                                                              >  5:331216/1‑149 (MQ=255)
                                                                          cGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGtt                                                                    <  5:228043/149‑1 (MQ=255)
                                                                          cGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGtt                                                                    <  8:91141/149‑1 (MQ=255)
                                                                          cGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGt                                                                     <  1:4643/148‑1 (MQ=255)
                                                                           ggTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGtt                                                                    >  2:209611/1‑148 (MQ=255)
                                                                                  cAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAggt                                                            >  5:305159/1‑149 (MQ=255)
                                                                                   aGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAggtg                                                           >  4:84845/1‑149 (MQ=255)
                                                                                     tGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCCCCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTCACCGGGGCAAAGCCGCCCTCGGCGAACCCCATGTCGTTGAACAggtggt                                                         >  7:151688/1‑149 (MQ=255)
                                                                                       ccGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTc                                                       <  3:84845/149‑1 (MQ=255)
                                                                                       ccGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTc                                                       >  1:7230/1‑149 (MQ=255)
                                                                                        cGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTcc                                                      <  6:305159/149‑1 (MQ=255)
                                                                                        cGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTcc                                                      <  1:209611/149‑1 (MQ=255)
                                                                                         gCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTc                                                       >  7:145736/1‑147 (MQ=255)
                                                                                            tCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGa                                                   <  4:146879/148‑1 (MQ=255)
                                                                                             ccAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGAt                                                  <  4:133958/148‑1 (MQ=255)
                                                                                              cAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATgg                                                <  4:6156/149‑1 (MQ=255)
                                                                                              cAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATg                                                 <  6:331216/148‑1 (MQ=255)
                                                                                               aGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAgg                                                             <  5:159525/135‑1 (MQ=255)
                                                                                               aGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAgg                                                             >  6:159525/1‑135 (MQ=255)
                                                                                                ggtcccgcCGGCCGAGGGCCTTCATCGTGGCGACGCCCTGCACCAGCTCTTCGCCAGTGAGGGTGCAGGCGGCCATGTCGCCAGAGCTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTCGCACAGGTGGTTCCAGATgggt                                              <  8:151688/144‑3 (MQ=255)
                                                                                                 gTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGtt                                              >  3:8627/1‑148 (MQ=255)
                                                                                                    gcgcCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCaa                                          <  6:84296/149‑1 (MQ=255)
                                                                                                          gccgAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCgg                                    <  8:145736/149‑1 (MQ=255)
                                                                                                               gTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATg                               <  2:53453/149‑1 (MQ=255)
                                                                                                                   cTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCAc                           >  7:224731/1‑149 (MQ=255)
                                                                                                                        tGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTc                      <  4:258471/149‑1 (MQ=255)
                                                                                                                           tGGCGAAGAGCTGGTGCAGGTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCAc                           >  1:213500/1‑141 (MQ=255)
                                                                                                                           tGGCGAAGAGCTGGTGCAGGTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCAc                           <  2:213500/141‑1 (MQ=255)
                                                                                                                               gAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGa               >  5:158959/1‑149 (MQ=255)
                                                                                                                               gAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGa               >  2:284793/1‑149 (MQ=255)
                                                                                                                               gAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGa               >  2:22127/1‑149 (MQ=255)
                                                                                                                                aaGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGAc              >  5:18157/1‑149 (MQ=255)
                                                                                                                                     cTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGttt                                             >  5:69845/1‑113 (MQ=255)
                                                                                                                                     cTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGttt                                             <  6:69845/113‑1 (MQ=255)
                                                                                                                                     cTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCcgcg         >  5:128307/1‑149 (MQ=255)
                                                                                                                                     cTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCcgcg         >  2:107965/1‑149 (MQ=255)
                                                                                                                                        gtgcAGGTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGCGccc      >  3:10817/4‑149 (MQ=255)
                                                                                                                                        cACCAGCTCTTCGCCGGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCgcgc        <  8:218809/147‑1 (MQ=255)
                                                                                                                                        cACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGCGcc       <  1:123711/148‑1 (MQ=255)
                                                                                                                                           cAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGCGCCCCCa   >  7:81707/1‑149 (MQ=255)
                                                                                                                                            aGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGaa                                                                                <  7:254647/71‑1 (MQ=255)
                                                                                                                                            aGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGaa                                                                                >  8:254647/1‑71 (MQ=255)
                                                                                                                                            aGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGCGCCCCCAc  <  2:7230/149‑1 (MQ=255)
                                                                                                                                            aGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGAGCCCCCAc  >  4:157530/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CACACGGCAGGGCCTTCTTGTGGCCTTCGTGGATGTCGGTGATCAACTGGTGCAGGTTGGAGAAGCCGATGTTGCGGTGCTCCAGTGCCGCTTCCAGGTAGCGCCGGCCGAGTGCCTTCATGTTGGCGAAGACCTGCACCAGCTCTTCGCCAGTGAGGTTGAAGTCGGCCATGTCGCCAGAGGTAACCGGGGCAAAGCCGACCTCGGCGAAACCCATTTCGTTGAACAGGTGGTTCCAGATGGTTTCAACGTCGGTAATGCCACGGGTCAGCGTGACCCGCGCCCCCAC  >  NC_002947/3923754‑3924042

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: