Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,195,990 T→C 5.7% A175A (GCA→GCG ybgJ ← allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,195,9900TC 5.7% 256.8 / 3.1 88A175A (GCA→GCGybgJallophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (42/40);  new base C (2/3);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCG  >  NC_002947/5195843‑5196138
                                                                                                                                                    |                                                                                                                                                       
aTCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGGGCTCCAGGTCACCGTGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACAGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTg                                                                                                                                                        <  7:163832/149‑1 (MQ=255)
 tCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGc                                                                                                                                                       >  8:215259/1‑149 (MQ=255)
   accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGcc                                                                                                                                                      <  8:370998/148‑1 (MQ=255)
      accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCGGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGt                                                                                                                                                  <  2:287339/149‑1 (MQ=255)
       cccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtg                                                                                                                                                 <  3:252692/149‑2 (MQ=255)
        ccATGGCCAGGAATTCGCTGCTCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGTGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAg                                                                                                                                                <  8:2143/149‑1 (MQ=255)
        ccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAg                                                                                                                                                >  4:288027/1‑149 (MQ=255)
          aTGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGtctt                                                                                                                                               >  2:162036/1‑145 (MQ=255)
             gCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgc                                                                                                                                           <  7:298228/149‑1 (MQ=255)
             gCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgc                                                                                                                                           >  2:224771/1‑149 (MQ=255)
               cAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                         >  7:118602/1‑149 (MQ=255)
                 ggAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGCTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCACGAACAGGTTCCAGCCACCCGCCGATACAGCCCGGTAGGCCGCCGt                                                                                                                                        >  4:120317/1‑148 (MQ=255)
                   aaTTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTg                                                                                                                                     <  6:30370/149‑1 (MQ=255)
                           cagcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCTAACAACCGCTCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                             <  7:155303/147‑1 (MQ=255)
                               gAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACCCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATAAGGTTCCAGCCACCCGGAGATACTGACGGGTAGCCCGCCGTCTGGCGCTAGGCaa                                                                                                                          >  8:232093/1‑148 (MQ=255)
                                   ccGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATAccc                                                                                                                    >  1:99673/1‑149 (MQ=255)
                                     ggCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCAc                                                                                                                  >  2:117314/1‑149 (MQ=255)
                                        aCGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCt                                                                                                               >  5:36038/1‑149 (MQ=255)
                                         cGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTg                                                                                                              <  8:118602/149‑1 (MQ=255)
                                          gAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGc                                                                                                             >  2:102672/1‑149 (MQ=255)
                                          gAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGc                                                                                                             >  7:107301/1‑149 (MQ=255)
                                           aaGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGt                                                                                                                                                  >  8:136114/1‑112 (MQ=255)
                                           aaGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGt                                                                                                                                                  <  7:136114/112‑1 (MQ=255)
                                            aGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGccc                                                                                                           <  8:366385/149‑1 (MQ=255)
                                             gcgcACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCccg                                                                                                          >  5:161267/1‑149 (MQ=255)
                                                cACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCccgccg                                                                                                       >  2:41134/1‑149 (MQ=255)
                                                 aCCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCcgccgc                                                                                                      >  2:164449/1‑149 (MQ=255)
                                                   ccGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCa                                                                                                   <  1:224771/149‑1 (MQ=255)
                                                   ccGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCa                                                                                                   <  8:107301/149‑1 (MQ=255)
                                                   ccGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACAGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGCGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCa                                                                                                   >  7:168394/1‑149 (MQ=255)
                                                      gTCATCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAcgc                                                                                                <  3:239811/149‑1 (MQ=255)
                                                         aCCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCt                                                                                             >  5:147326/1‑149 (MQ=255)
                                                                 gCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                       <  1:41134/148‑1 (MQ=255)
                                                                 gCAAGAGGCTGTAACCCTCACGCTCACCGTCGAACAACCGCACCGGGGTGCGCCCGCTCAGGGTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCGCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                      <  4:256133/148‑1 (MQ=255)
                                                                  cAAGTGGCTGTAACCCTCAGGCTCACGGTCTAACTACCGCACCGGGGTGCGCCCGATCAGGTTCGAGCCACCCGGCGATACTGCCGGGTATGCCGCGGTCTGGCGCTCGGCAATACCCACGCTGCCCGGCGGCACGTGCTTGCGTggg                                                                                     <  6:36038/148‑1 (MQ=255)
                                                                  cAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                                    <  7:230625/149‑1 (MQ=255)
                                                                   aaGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                                    >  5:23727/1‑148 (MQ=255)
                                                                     gagGCTGTAACCCTCACGCTCACGGTCGAACACCCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCCCGCGCTTGCGTGGGGTAc                                                                                 >  4:25396/1‑149 (MQ=255)
                                                                       ggCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACt                                                                                <  2:125861/148‑1 (MQ=255)
                                                                             aaCCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGa                                                                         >  6:244827/1‑149 (MQ=255)
                                                                              aCCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAg                                                                        >  3:231118/1‑149 (MQ=255)
                                                                               cccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAgg                                                                       >  4:234519/1‑149 (MQ=255)
                                                                               cccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAgg                                                                       <  3:6880/149‑1 (MQ=255)
                                                                               cccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAgg                                                                       >  7:105270/1‑149 (MQ=255)
                                                                                ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAggg                                                                      >  3:204284/1‑149 (MQ=255)
                                                                                 ctcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                     >  5:137996/1‑149 (MQ=255)
                                                                                    acgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTg                                                                   <  4:231118/148‑1 (MQ=255)
                                                                                       ctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCa                                                               <  3:288027/149‑1 (MQ=255)
                                                                                          acgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGgc                                                            <  5:244827/149‑1 (MQ=255)
                                                                                          acgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGgc                                                            >  2:62910/1‑149 (MQ=255)
                                                                                            ggTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGgcgc                                                          >  8:86914/1‑149 (MQ=255)
                                                                                                gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgt                                                      <  6:23727/149‑1 (MQ=255)
                                                                                                 aacaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                     >  8:308823/1‑149 (MQ=255)
                                                                                                 aacaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                     >  3:90196/1‑149 (MQ=255)
                                                                                                  acaacCGCCCCGGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                     <  2:165572/149‑1 (MQ=255)
                                                                                                  acaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                     >  5:276204/1‑148 (MQ=255)
                                                                                                    aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTgg                                                                  >  8:149997/1‑133 (MQ=255)
                                                                                                    aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTgg                                                                  <  7:149997/133‑1 (MQ=255)
                                                                                                    aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGAc                                                  <  2:99673/149‑1 (MQ=255)
                                                                                                      ccGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCa                                               <  5:65800/149‑1 (MQ=255)
                                                                                                         cACCGGGGTGCGCCCGATCGGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGc                                            >  3:122972/1‑149 (MQ=255)
                                                                                                         cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCTCCGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAgg                                             <  1:102672/149‑1 (MQ=255)
                                                                                                          aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGcc                                           >  8:165929/1‑149 (MQ=255)
                                                                                                              ggggTGCGCCCGATCAGGGTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCAt                                        <  5:383157/149‑1 (MQ=255)
                                                                                                               gggTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCAt                                        >  4:348003/1‑148 (MQ=255)
                                                                                                                     gcCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCACTACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGcc                                  >  5:273459/1‑148 (MQ=255)
                                                                                                                         gATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAc                             >  6:307153/1‑149 (MQ=255)
                                                                                                                          aTCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAcc                            >  8:304318/1‑149 (MQ=255)
                                                                                                                           tCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAccc                           <  4:204284/149‑1 (MQ=255)
                                                                                                                              ggTTCCAGCCACCCGGCGATACTGTCGGGTAGGCCGGCGTCTGGCGGTCGGCAATACGCTCGCTGCCCGCCGCGACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGGCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        <  8:168394/149‑1 (MQ=255)
                                                                                                                              ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        <  1:62910/149‑1 (MQ=255)
                                                                                                                              ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        >  8:161203/1‑149 (MQ=255)
                                                                                                                              ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        <  7:304318/149‑1 (MQ=255)
                                                                                                                              ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        <  7:86914/149‑1 (MQ=255)
                                                                                                                              ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        <  5:55733/149‑1 (MQ=255)
                                                                                                                               gTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                       >  7:316207/1‑149 (MQ=255)
                                                                                                                               gTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                       <  7:161203/149‑1 (MQ=255)
                                                                                                                                ttCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGc                      <  2:198175/149‑1 (MQ=255)
                                                                                                                                  ccAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCg                     <  3:234519/148‑1 (MQ=255)
                                                                                                                                   cAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGaa                   <  7:165929/149‑1 (MQ=255)
                                                                                                                                    aGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                      >  3:292208/1‑82 (MQ=255)
                                                                                                                                    aGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                      <  4:292208/82‑1 (MQ=255)
                                                                                                                                    aGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAg                  >  4:53867/1‑149 (MQ=255)
                                                                                                                                       taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAc                                                                                                  >  6:47916/2‑67 (MQ=37)
                                                                                                                                       taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAc                                                                                                  <  5:47916/66‑1 (MQ=37)
                                                                                                                                        aCCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc                >  8:295268/1‑147 (MQ=255)
                                                                                                                                        aCCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGAGCTGGCCAGGCGCTCGTCGACCAAGGCCCATGAAGCCGAAACCCGGTGCGAAGccc               <  1:53039/149‑1 (MQ=255)
                                                                                                                                          ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg            >  5:125814/1‑149 (MQ=255)
                                                                                                                                                    tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCg  >  8:32711/1‑149 (MQ=255)
                                                                                                                                                    tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCg  >  2:299770/1‑149 (MQ=255)
                                                                                                                                                    tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCg  >  5:362392/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
ATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCG  >  NC_002947/5195843‑5196138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: