Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,871,259 C→G 6.9% P164A (CCC→GCC) ‡ xdhC → xanthine dehydrogenase accessory factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,871,2590CG 6.9% 231.0 / 6.5 72P164A (CCC→GCC) ‡xdhCxanthine dehydrogenase accessory factor
Reads supporting (aligned to +/- strand):  ref base C (37/30);  new base G (2/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

AGGGCTCTGGTACCCTTGCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGAGCACCGCCTGCGCGAAC  >  NC_002947/4871115‑4871406
                                                                                                                                                |                                                                                                                                                   
aGGGCTCTGGTACCCTTGCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCACGGTCGGCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCgg                                                                                                                                                 <  3:235564/149‑1 (MQ=255)
       tGGTACCCTTGCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTAc                                                                                                                                          >  1:362898/1‑149 (MQ=255)
         gTACCCTTGCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGGGCGGCAggtca                                                                                                                                          <  4:290595/147‑5 (MQ=255)
            cccTTGCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCa                                                                                                                                      <  7:320065/148‑1 (MQ=255)
            cccTTACTCGCCGCGCTCCCCTGCGGGGTGCGCTGGATCGACTCTCGTGAACTGCAATTCCCCGAACTCATCCCCAATGGGGTGACCAATGTCGTCAGTGATGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGGAGCTACTGCa                                                                                                                                      <  4:113506/148‑1 (MQ=255)
                 gCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGt                                                                                                                                  <  5:3446/147‑1 (MQ=255)
                   tCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATg                                                                                                                              >  3:282546/1‑149 (MQ=255)
                     gccgcGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGGGCGGCAggtca                                                                                                                                          >  3:116122/1‑131 (MQ=255)
                     gccgcGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGGGCGGCAggtca                                                                                                                                          <  4:116122/135‑5 (MQ=255)
                     gccgcGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTc                                                                                                                                 <  4:318119/144‑1 (MQ=255)
                     gccgcGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGa                                                                                                                             <  7:84757/148‑1 (MQ=255)
                     gccgcGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCAAGat                                                                                                                            >  3:318119/1‑148 (MQ=255)
                            tCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAac                                                                                                                     >  7:19525/1‑149 (MQ=255)
                            tCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAac                                                                                                                     >  3:160939/1‑149 (MQ=255)
                              cccTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGATCCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCCTGACCCACAgcca                                                                                                                   >  8:196719/1‑149 (MQ=255)
                                cTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAaccacc                                                                                                                 <  6:69906/149‑1 (MQ=255)
                                       gTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGc                                                                                                                                                >  7:108156/1‑111 (MQ=255)
                                       gTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGc                                                                                                                                                <  8:108156/111‑1 (MQ=255)
                                        tGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGAc                                                                                                         <  4:282546/149‑1 (MQ=255)
                                        tGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGAc                                                                                                         <  8:19525/149‑1 (MQ=255)
                                             tGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTgg                                                                                                     >  5:356806/1‑148 (MQ=255)
                                              ggATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGaa                                                                                                   >  3:11877/1‑149 (MQ=255)
                                                    aCTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGGGCGGCAggtc                                                                                                                                           >  5:391137/1‑100 (MQ=255)
                                                    aCTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGGGCGGCAggtc                                                                                                                                           <  6:391137/103‑4 (MQ=255)
                                                     cTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAactgactg                                                                                            >  2:22824/1‑149 (MQ=255)
                                                        gcgTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTgccg                                                                                         >  2:285550/1‑149 (MQ=255)
                                                        gcgTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTgcc                                                                                          <  2:346845/148‑1 (MQ=255)
                                                          gTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTgccgc                                                                                        >  4:255819/1‑148 (MQ=255)
                                                               cAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATtct                                                                                  >  5:321330/1‑149 (MQ=255)
                                                                  gAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAactga                                                                                               >  3:49617/1‑133 (MQ=255)
                                                                  gAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAactga                                                                                               <  4:49617/133‑1 (MQ=255)
                                                                  gAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCaa                                                                               <  4:11877/149‑1 (MQ=255)
                                                                       ccccGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAgcgc                                                                           >  3:326419/1‑148 (MQ=255)
                                                                        cccGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCaa                                                                         >  4:20424/1‑149 (MQ=255)
                                                                            aaCTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGAt                                                                     <  4:160939/149‑1 (MQ=255)
                                                                               tCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTAAAGCGCAACGAttt                                                                   >  4:202165/1‑148 (MQ=255)
                                                                                 aTCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAactgact                                                                                             >  5:135753/1‑120 (MQ=255)
                                                                                 aTCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAactgact                                                                                             <  6:135753/120‑1 (MQ=255)
                                                                                 aTCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTc                                                                  <  5:79233/147‑1 (MQ=255)
                                                                                 aTCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCa                                                                 <  2:362898/148‑1 (MQ=255)
                                                                                    cccAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTg                                                             >  8:106119/1‑149 (MQ=255)
                                                                                     ccAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTgg                                                            >  4:128877/1‑149 (MQ=255)
                                                                                     ccAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTgg                                                            >  8:125961/1‑149 (MQ=255)
                                                                                       aaTGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGc                                                                                                              >  3:6239/1‑97 (MQ=255)
                                                                                       aaTGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGc                                                                                                              <  4:6239/97‑1 (MQ=255)
                                                                                       aaTGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGtt                                                          >  6:239313/1‑149 (MQ=255)
                                                                                           gtgTGACCAAGGTCGTCATGGAGTAACCGTTCTACGAAGTGTCTGACCTGCCGCCCGGCAGCTACTGCATCGGCATGACCCACAACTACCAGCTGTTCCTGGAACGGACTTCCGGCATTCTCAAGCGCAACGATTTCACCTGGTTCGGc                                                      <  2:263966/147‑1 (MQ=255)
                                                                                            ggTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGc                                                      <  6:321330/148‑1 (MQ=255)
                                                                                               gACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGa                                                  >  5:190250/1‑149 (MQ=255)
                                                                                                  cAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCg                                               >  6:80882/1‑149 (MQ=255)
                                                                                                   aaGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATc                                                <  8:266121/147‑1 (MQ=255)
                                                                                                   aaGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATc                                                >  7:266121/1‑147 (MQ=255)
                                                                                                      gtcgtcAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTc                                           >  2:129950/1‑149 (MQ=255)
                                                                                                      gtcgtcAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTc                                           >  5:20585/1‑149 (MQ=255)
                                                                                                       tcgtcAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCg                                          <  1:285550/149‑1 (MQ=255)
                                                                                                         gtcAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGaa                                        <  7:106119/149‑1 (MQ=255)
                                                                                                           cAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGa                                      <  4:326419/149‑1 (MQ=255)
                                                                                                            aGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGa                                      <  1:22824/148‑1 (MQ=255)
                                                                                                            aGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGAc                                     >  2:239516/1‑149 (MQ=255)
                                                                                                              tGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGAcc                                    >  1:202070/1‑148 (MQ=255)
                                                                                                                  gAACCGGTCGACGAAGTCGCTGACCTGCAGCCCGGCAGCTACTGCATCGTCATGACCAACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCCCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACg                               >  2:78080/1‑149 (MQ=255)
                                                                                                                   aaCCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACgg                              >  3:215614/1‑149 (MQ=255)
                                                                                                                    aCCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACggg                             <  7:125961/149‑1 (MQ=255)
                                                                                                                        gTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCaa                         >  1:325762/1‑149 (MQ=255)
                                                                                                                           gacgaAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGtt                      >  4:47101/1‑149 (MQ=255)
                                                                                                                           gacgaAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGtt                      >  2:272474/1‑149 (MQ=255)
                                                                                                                              gaAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGa                   <  4:215614/149‑1 (MQ=255)
                                                                                                                              gaAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGa                   <  6:356806/149‑1 (MQ=255)
                                                                                                                                aGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCAcc                                                               >  7:237028/1‑103 (MQ=255)
                                                                                                                                aGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCAcc                                                               <  8:237028/103‑1 (MQ=255)
                                                                                                                                 gTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGAGCa                >  2:209783/1‑149 (MQ=255)
                                                                                                                                         ccTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGAGCACCGCCTgc        <  3:255819/149‑1 (MQ=255)
                                                                                                                                            gccgccCGGCAGTTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGaa                                        >  8:311726/1‑114 (MQ=255)
                                                                                                                                            gccgccCGGCAGTTACCGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGaa                                        <  7:311726/114‑1 (MQ=255)
                                                                                                                                             ccgccCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGAGCACCGCCTgcgc      <  3:128877/147‑1 (MQ=255)
                                                                                                                                               gccCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGAGCACCGCCTGCGCGAAc  >  1:255189/1‑149 (MQ=255)
                                                                                                                                                cccGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGAACAAACGGGTCAAGTTCGAGCACCGCCTGCGCGAAc  >  4:126199/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGGGCTCTGGTACCCTTGCTCGCCGCGCTCCCCTGCCGGGTGCGCTGGATCGACTCGCGTGAACAGGAATTCCCCGAACTCATCCCCAATGGGGTGACCAAGGTCGTCAGTGAGGAACCGGTCGACGAAGTCGCTGACCTGCCGCCCGGCAGCTACTGCATCGTCATGACCCACAACCACCAGCTCGACCTGGAACTGACTGCCGCCATTCTCAAGCGCAACGATTTCACCTGGTTCGGCCTGATCGGCTCGAAGACCAAACGGGTCAAGTTCGAGCACCGCCTGCGCGAAC  >  NC_002947/4871115‑4871406

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: