Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,635,440 C→A 8.1% intergenic (+30/+87) PP_2304 → / ← PP_2305 peptidyl‑prolyl cis‑trans isomerase D/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,635,4400CA 8.1% 265.8 / 11.7 86intergenic (+30/+87)PP_2304/PP_2305peptidyl‑prolyl cis‑trans isomerase D/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/35);  new base A (2/5);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CGCCACCGACGAAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGA  >  NC_002947/2635293‑2635588
                                                                                                                                                   |                                                                                                                                                    
cgccACCGACGAAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGATGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCt                                                                                                                                                     >  2:203339/1‑149 (MQ=255)
      cgacgaAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCg                                                                                                                                               >  5:154104/1‑149 (MQ=255)
       gacgaAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCACGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCgg                                                                                                                                              <  1:85148/149‑1 (MQ=255)
       gacgaAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCgg                                                                                                                                              >  8:48389/1‑149 (MQ=255)
        acgaAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGc                                                                                                                                             >  7:101531/1‑149 (MQ=255)
             gagaAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGtt                                                                                                                                        >  3:23466/1‑149 (MQ=255)
              agagGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCACGAGGCGGCCTGttt                                                                                                                                       <  3:332823/149‑1 (MQ=255)
              agaAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGttt                                                                                                                                       >  5:164876/1‑149 (MQ=255)
              agaAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCACGAGGCGGCCTGttt                                                                                                                                       <  1:235560/149‑1 (MQ=255)
               gaAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCACGAGGCGGCCTGTTTc                                                                                                                                      <  1:181726/149‑1 (MQ=255)
                aaGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCg                                                                                                                                     >  5:351285/1‑149 (MQ=255)
                  gcagcaGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCACGAGGCGGCCTGTTTcttc                                                                                                                                   <  5:248577/149‑4 (MQ=255)
                    agcagATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGtttt                                                                                                                                 >  5:25360/1‑149 (MQ=255)
                     gcagATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTc                                                                                                                                <  3:119237/149‑1 (MQ=255)
                       agATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGcc                                                                                                                             <  2:141386/149‑1 (MQ=255)
                        gATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGc                                                                                                                              >  7:322348/1‑148 (MQ=255)
                              ccgcTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCcgc                                                                                                                       <  4:23466/149‑1 (MQ=255)
                              ccgcTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCACGAGGCGGCCTGTTTctt                                                                                                                                    <  4:233150/136‑3 (MQ=255)
                              ccgcTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCACGAGGCGGCCTGTTTctt                                                                                                                                    >  3:233150/1‑134 (MQ=255)
                                    ccTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCa                                                                                                                  <  6:321516/148‑1 (MQ=255)
                                    ccTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTACGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGcc                                                                                                                         >  5:341011/1‑141 (MQ=255)
                                    ccTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTACGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGcc                                                                                                                         <  6:341011/141‑1 (MQ=255)
                                     cTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACg                                                                                                                <  6:154104/149‑1 (MQ=255)
                                     cTCGCGTCACGTGCTGGCCAGCAGGACTTAGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACg                                                                                                                >  1:221612/1‑149 (MQ=255)
                                           tCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCAtt                                                                                                          >  1:52912/1‑149 (MQ=255)
                                              cGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTc                                                                                                       <  8:101531/149‑1 (MQ=255)
                                                  cTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGaa                                                                                                    <  6:351285/148‑1 (MQ=255)
                                                     gCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGt                                                                                                 >  7:100785/1‑148 (MQ=255)
                                                      ccAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTca                                                                                               >  2:59585/1‑149 (MQ=255)
                                                          cagGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACGCGTTACTAAGGTTTTAGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCATGTTTCTTTTTCGAAAGCCGCGAGCACGTTCATTTTCGATGGTCACAGt                                                                                           >  6:72172/1‑149 (MQ=255)
                                                              aCTTTGCCGCTTACCGCAAGTAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGCGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAg                                                                                                   >  1:45923/1‑137 (MQ=255)
                                                              aCTTTGCCGCTTACCGCAAGTAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGCGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAg                                                                                                   <  2:45923/137‑1 (MQ=255)
                                                               cTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAg                                                                                                   <  1:278151/136‑1 (MQ=255)
                                                               cTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATCAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAg                                                                                                   >  2:278151/1‑136 (MQ=255)
                                                                tttGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTca                                                                                               <  5:254456/139‑1 (MQ=255)
                                                                tttGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTca                                                                                               >  6:254456/1‑139 (MQ=255)
                                                                          aCCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCTTTGTCGAAGGTCACAGTGCCTTCAGGGGTgg                                                                             <  1:203339/147‑1 (MQ=255)
                                                                          aCCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGgcg                                                                           >  7:22086/1‑149 (MQ=255)
                                                                          aCCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGgc                                                                            <  4:42368/148‑1 (MQ=255)
                                                                          aCCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGgc                                                                            <  8:322348/148‑1 (MQ=255)
                                                                             gCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTa                                                                        <  8:100785/149‑1 (MQ=255)
                                                                                agcagcTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAAcc                                                                     <  2:52912/149‑1 (MQ=255)
                                                                                  cagcTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAg                                                                   >  4:176062/1‑149 (MQ=255)
                                                                                   agcTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAgg                                                                  >  1:263070/1‑149 (MQ=255)
                                                                                    gcTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGt                                                                 <  7:48389/149‑1 (MQ=255)
                                                                                      tGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGc                                                               >  8:82502/1‑149 (MQ=255)
                                                                                       ggAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTaa                                                                       <  2:295124/140‑1 (MQ=255)
                                                                                       ggAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAAcagatagg                                                               >  1:295124/1‑141 (MQ=255)
                                                                                          aGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTAcc                                                           <  6:25360/149‑1 (MQ=255)
                                                                                              aaTGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGCGTGGCGCTAACCAGGTGCTACCGCtt                                                       >  4:296804/1‑149 (MQ=255)
                                                                                               aTGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTc                                                      <  2:188267/149‑1 (MQ=255)
                                                                                                   aGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCAccc                                                  >  5:283190/1‑149 (MQ=255)
                                                                                                       aTCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCtt                                               >  6:263581/1‑148 (MQ=255)
                                                                                                            ccGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAAc                                         >  7:128304/1‑149 (MQ=255)
                                                                                                              gTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTg                                       >  7:35627/1‑149 (MQ=255)
                                                                                                                 aCTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTc                                     <  3:176062/148‑1 (MQ=255)
                                                                                                                 aCTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTc                                     >  4:94158/1‑148 (MQ=255)
                                                                                                                  cTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTc                                              <  5:385730/138‑1 (MQ=255)
                                                                                                                  cTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTc                                              >  6:385730/1‑138 (MQ=255)
                                                                                                                  cTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTc                                     >  1:87634/1‑147 (MQ=255)
                                                                                                                   tAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCa                                                                    >  7:106279/1‑115 (MQ=255)
                                                                                                                   tAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCa                                                                    <  8:106279/115‑1 (MQ=255)
                                                                                                                   tAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTAc                                                            >  3:320010/1‑123 (MQ=255)
                                                                                                                   tAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTAc                                                            <  4:320010/123‑1 (MQ=255)
                                                                                                                      ggTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGgc                                                                            >  6:371759/1‑104 (MQ=255)
                                                                                                                      ggTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGgc                                                                            <  5:366379/104‑1 (MQ=255)
                                                                                                                      ggTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGgc                                                                            <  5:371759/104‑1 (MQ=255)
                                                                                                                      ggTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTAACAGTGCCTTCAGGGGTGgc                                                                            >  6:366379/1‑104 (MQ=255)
                                                                                                                        tgtTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTAt                             >  5:137269/1‑149 (MQ=255)
                                                                                                                             gTTGGATGAAGATACAGGCCTCCTGGTGCGGGCTGTTTCGTTTTCGCCAGCCGCGAGGACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGt                        <  2:263070/149‑1 (MQ=255)
                                                                                                                             gTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGt                        <  6:164876/149‑1 (MQ=255)
                                                                                                                              ttGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTg                       >  1:140079/1‑149 (MQ=255)
                                                                                                                                gCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATcgc                          >  5:172001/1‑144 (MQ=255)
                                                                                                                                gCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATcgc                          <  6:172001/144‑1 (MQ=255)
                                                                                                                                gCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGc                      >  5:316684/1‑148 (MQ=255)
                                                                                                                                 gaaaacgaaACAGGCCGCACGAGGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCaa                    >  3:87706/7‑149 (MQ=39)
                                                                                                                                  aTGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATcgc                          >  7:336514/1‑142 (MQ=255)
                                                                                                                                  aTGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATcgc                          <  8:336514/142‑1 (MQ=255)
                                                                                                                                     aagaaACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATcgc                          >  7:276680/1‑139 (MQ=255)
                                                                                                                                     aagaaACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATcgc                          <  8:276680/139‑1 (MQ=255)
                                                                                                                                     aagaaACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGt                        <  5:324450/141‑1 (MQ=255)
                                                                                                                                     aagaaACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGt                        >  6:324450/1‑141 (MQ=255)
                                                                                                                                     aagaaACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAAtctc                >  2:249332/1‑149 (MQ=255)
                                                                                                                                      agaaACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAAtctct               <  3:94158/149‑1 (MQ=255)
                                                                                                                                         aaCAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGc            <  2:87634/149‑1 (MQ=255)
                                                                                                                                          aCAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAggg                                                                                 >  6:297463/1‑79 (MQ=255)
                                                                                                                                               ccgccTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAgg      <  3:296804/149‑1 (MQ=255)
                                                                                                                                                cgccTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAggg     >  8:211791/1‑149 (MQ=255)
                                                                                                                                                   cTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGa  >  6:302619/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCACCGACGAAGAGAAGCAGCAGATTCGCCGCTACCTCGCGTCACGTGCTGGCCAGCAGGACTTTGCCGCTTACCGCAAGCAGCTGGAAGCCAATGCAGACATCACCCGTTACTAAGGTGTTCGTTGCATGAAGAAACAGGCCGCCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGA  >  NC_002947/2635293‑2635588

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: