Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,623 T→A 5.3% K183I (AAA→ATA)  fkl ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947243,6230TA 5.3% 227.2 / 2.1 76K183I (AAA→ATA) fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGCCGCCCGCT  >  NC_002947/243476‑243769
                                                                                                                                                    |                                                                                                                                                   
gAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCAAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCAttt                                                                                                                                                    >  6:99013/1‑149 (MQ=255)
   gagGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGgc                                                                                                                                                 <  7:301330/149‑1 (MQ=255)
   gagGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCAAATGgc                                                                                                                                                 >  8:297674/1‑149 (MQ=255)
       gtgcagtgcaGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGCAGGAATTACC‑AGGCGCCATTTGGCGCCg                                                                                                                                             <  3:277007/149‑1 (MQ=255)
        tgcagtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCAAATGGCGcctg                                                                                                                                            <  2:218386/149‑3 (MQ=255)
        tgcagtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCAAATGGCGcctg                                                                                                                                            <  7:297674/149‑3 (MQ=255)
          cagtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGgt                                                                                                                                           <  2:166731/148‑1 (MQ=255)
             tgcagGCCTGAGGTTGGCGTCAGTCGCGAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAAATACC‑AGGCGCCATTTGGCGCCGGTGTGc                                                                                                                                       <  7:172253/149‑1 (MQ=255)
             tgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGtgtg                                                                                                                                        >  5:385890/1‑148 (MQ=255)
             tgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGtgtg                                                                                                                                        >  3:188603/1‑148 (MQ=255)
                  gCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCAtttt                                                                                                                                  <  6:385890/149‑1 (MQ=255)
                           tGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCg                                                                                                                                             <  4:50329/129‑1 (MQ=255)
                           tGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCg                                                                                                                                             >  3:50329/1‑129 (MQ=255)
                             gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCg                                                                                                                        >  2:155786/1‑148 (MQ=255)
                              cGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCg                                                                                                                                <  7:319684/139‑1 (MQ=255)
                              cGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCg                                                                                                                                >  8:319684/1‑139 (MQ=255)
                                    tCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGc                                                                                                                  >  6:241701/1‑147 (MQ=255)
                                    tCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCa                                                                                                                <  6:124264/149‑1 (MQ=255)
                                      gCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGc                                                                                                              <  3:184789/149‑1 (MQ=255)
                                       ccAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTTGAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGcc                                                                                                             <  8:313358/149‑1 (MQ=255)
                                           ccGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTc                                                                                                          >  1:127114/1‑148 (MQ=255)
                                             gCCAGCAGTTCGATCTCAAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGAt                                                                                                       <  5:99013/149‑1 (MQ=255)
                                                   agTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGa                                                                                                 <  4:188603/149‑1 (MQ=255)
                                                    gTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGTGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGaa                                                                                                >  5:5005/1‑149 (MQ=255)
                                                     ttCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAAt                                                                                               >  3:404744/1‑149 (MQ=255)
                                                           ctcGTATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGTAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAgg                                                                                          <  4:338698/148‑1 (MQ=255)
                                                            tcGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGc                                                                                                                  <  7:156350/123‑1 (MQ=255)
                                                            tcGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGc                                                                                                                  >  8:156350/1‑123 (MQ=255)
                                                              gAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAgg                                                                                          <  4:404744/145‑1 (MQ=255)
                                                              gAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGtt                                                                                        <  2:415836/147‑1 (MQ=255)
                                                              gAATACCAAGGGGGTGTAGGGTGCAATCAGCTCTCCCGCCCCCTCGGCGCAATAGGATAGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTACGATAACCGAATCCAGGTTg                                                                                       >  6:13981/1‑148 (MQ=255)
                                                                  aCCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCa                                                                                  <  3:312157/149‑1 (MQ=255)
                                                                   ccAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCa                                                                                  <  1:54476/148‑1 (MQ=255)
                                                                                gggTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggt                                                                    >  1:57681/1‑149 (MQ=255)
                                                                                  gTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttc                                                                  >  2:265499/1‑149 (MQ=255)
                                                                                   tGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttct                                                                 >  6:76851/1‑149 (MQ=255)
                                                                                      aaTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGTAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                                <  1:162474/147‑1 (MQ=255)
                                                                                      aaTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctggt                                                              >  1:202270/1‑149 (MQ=255)
                                                                                      aaTCAGGTCTCCCGCACCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                                >  2:162474/1‑147 (MQ=255)
                                                                                       aTCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACT‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTc                                                             >  5:329593/1‑149 (MQ=255)
                                                                                              ctcCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCGAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAatat                                                       >  1:90583/1‑149 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACT‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGa                                                    <  6:329593/148‑1 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGCACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAg                                                   >  6:37912/1‑149 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATAtt                                                      <  2:202270/146‑1 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTg                                                     <  5:241701/147‑1 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAg                                                   <  1:265499/149‑1 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAg                                                   <  6:400587/149‑1 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAg                                                   <  2:90583/149‑1 (MQ=255)
                                                                                                 ccGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGAGAACCGAATCCAGGTTGATCGATTGGGGGGTCTGGGTCTGGTCGAATATTGAg                                                   >  7:297955/1‑149 (MQ=255)
                                                                                                    cccccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCa                                                                                            <  1:118253/105‑1 (MQ=255)
                                                                                                    cccccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCa                                                                                            >  2:118253/1‑105 (MQ=255)
                                                                                                    cccccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                                <  1:153605/133‑1 (MQ=255)
                                                                                                    cccccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGtctggttctg                                                                >  2:153605/1‑133 (MQ=255)
                                                                                                       ccTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑Tcc                                            >  3:421269/1‑149 (MQ=255)
                                                                                                         tCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCgg                                          >  7:41490/1‑149 (MQ=255)
                                                                                                         tCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCgg                                          >  2:2815/1‑149 (MQ=255)
                                                                                                                cATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTg                                   <  1:155786/149‑1 (MQ=255)
                                                                                                                 aTAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAACCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGc                                  >  4:183647/1‑149 (MQ=255)
                                                                                                                  tAGGCTTGCGCTGAAGGAATTACC‑AGGCGGCATTTGGCGCCGGTGTGAATTTTCTGCAAGGCGACCTGGCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGcc                                 <  5:13981/149‑1 (MQ=255)
                                                                                                                  tAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGTGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGcc                                 <  6:5005/149‑1 (MQ=255)
                                                                                                                    ggCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGccc                                <  4:421269/148‑1 (MQ=255)
                                                                                                                     gCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCca                               >  8:284878/1‑148 (MQ=255)
                                                                                                                       ttGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCcaca                             >  7:84648/1‑148 (MQ=255)
                                                                                                                         gcgcTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAg                          >  4:383165/1‑149 (MQ=255)
                                                                                                                             tGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAgcgc                       >  1:332780/1‑148 (MQ=255)
                                                                                                                                 ggAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCgcacttg                  <  8:84648/149‑1 (MQ=255)
                                                                                                                                   aaTTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCgcacttgca                >  6:362370/1‑149 (MQ=255)
                                                                                                                                   aaTTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCgcacttgca                >  3:12824/1‑149 (MQ=255)
                                                                                                                                    aTTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCgcacttgcac               >  3:384976/1‑149 (MQ=255)
                                                                                                                                       aCC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACttt            >  1:366275/1‑149 (MQ=255)
                                                                                                                                        cc‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAACCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTg           <  3:183647/149‑1 (MQ=255)
                                                                                                                                         c‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTgc          >  5:113347/1‑149 (MQ=255)
                                                                                                                                           cggCGCCAAATGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTgcc         >  7:319756/2‑149 (MQ=255)
                                                                                                                                           aGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCGCTCCGGCAATCTGCCCACATAGCGCACTTGCACTTTg           <  1:71880/148‑1 (MQ=255)
                                                                                                                                           aGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTg           >  2:71880/1‑147 (MQ=255)
                                                                                                                                               gcCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGccgccc     <  4:384976/149‑1 (MQ=255)
                                                                                                                                                  aTTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGTACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGCCGCCCGCt  >  1:356960/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACC‑AGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCAGCCTTCGATAACCGAATCCAGGTTGAACCATTGGGGGGTCTGGTTCTGGTCGAATATTGAGCCG‑TCCGGCAATCTGCCCACATAGCGCACTTGCACTTTGCCGCCCGCT  >  NC_002947/243476‑243769

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: ‑