Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,172,360 G→A 7.7% G216G (GGC→GGT fadD‑I ← long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,172,3600GA 7.7% 300.2 / 10.4 104G216G (GGC→GGTfadD‑Ilong‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (54/42);  new base A (3/5);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 117 total aligned reads displayed.

ACGCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACGCT  >  NC_002947/5172212‑5172509
                                                                                                                                                    |                                                                                                                                                      
aCGCATAGATGTGGTAGATCGGCAGCGGGGTTATGAGGATTTCCCCGGCTTTTTCCTTCGTGATATCTAGCATTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                                       <  1:1897147/150‑1 (MQ=19)
aCGCATAGATGTGGTACAGCCGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                                       <  2:2238919/150‑1 (MQ=255)
  gCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgccg                                                                                                                                                     <  1:1207321/150‑1 (MQ=255)
    aTAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCggg                                                                                                                                                   >  1:722328/1‑149 (MQ=255)
    aTAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCggg                                                                                                                                                   >  1:1821303/1‑149 (MQ=255)
    aTAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGgt                                                                                                                                                   >  1:2804120/1‑150 (MQ=255)
       gATGTTTTACATCGGCAGCGGGGTGATGAGGATTTCGCACCATTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGCTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtg                                                                                                                                                <  2:1244738/150‑2 (MQ=37)
        aTGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtgg                                                                                                                                               <  1:1231046/150‑3 (MQ=255)
        aTGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtgg                                                                                                                                               <  1:164169/150‑3 (MQ=255)
         tgtgGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                                       >  2:2489955/1‑141 (MQ=255)
         tgtgGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGACTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGCGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                                       <  1:2489955/141‑1 (MQ=255)
           tgGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgc                                                                                                                                            >  1:2672174/1‑150 (MQ=255)
             gTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgcag                                                                                                                                          >  1:1314030/1‑150 (MQ=255)
             gTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgcag                                                                                                                                          >  1:1556856/1‑150 (MQ=255)
               aCAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgcagca                                                                                                                                        >  1:739516/1‑150 (MQ=255)
                            ggTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGtt                                                                                                                           <  2:3002243/150‑1 (MQ=255)
                                atgaGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGcc                                                                                                                      >  2:1625788/1‑150 (MQ=255)
                                  gaGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGcc                                                                                                                                   >  1:63729/1‑136 (MQ=255)
                                  gaGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGcc                                                                                                                                   <  2:63729/136‑1 (MQ=255)
                                  gaGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTg                                                                                                                    >  2:911381/1‑150 (MQ=255)
                                   aGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGc                                                                                                                   <  1:2016618/150‑1 (MQ=255)
                                   aGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGc                                                                                                                   <  1:2104737/150‑1 (MQ=255)
                                    ggATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                                       <  1:438540/114‑1 (MQ=255)
                                        ttCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCgggct                                                                                                              <  2:1584426/150‑3 (MQ=255)
                                         tCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTg                                                                                                             <  1:1874049/150‑1 (MQ=255)
                                           gCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGc                                                                                                           <  1:253370/150‑1 (MQ=255)
                                                 ttCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgcag                                                                                                                                          <  2:1388966/114‑1 (MQ=255)
                                                 ttCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGt                                                                                                     <  1:1777337/150‑1 (MQ=255)
                                                    gTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCAc                                                                                                  >  1:1933853/1‑150 (MQ=255)
                                                    gTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCAc                                                                                                  <  1:911381/150‑1 (MQ=255)
                                                    gTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCAc                                                                                                  >  2:817377/1‑150 (MQ=255)
                                                        aGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCTGCGGGACAGCATGGCGCCTTTGGTCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGc                                                                                              >  2:3049825/1‑150 (MQ=255)
                                                         ggTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCt                                                                                             <  2:2993125/150‑1 (MQ=255)
                                                          gTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTg                                                                                            >  2:1661711/1‑150 (MQ=255)
                                                          gTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTg                                                                                            >  1:812565/1‑150 (MQ=255)
                                                          gTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTg                                                                                            <  2:1933853/150‑1 (MQ=255)
                                                           ttGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTgg                                                                                           <  2:2672174/150‑1 (MQ=255)
                                                                gCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCtt                                                                                      >  1:2249022/1‑150 (MQ=255)
                                                                 cccATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTg                                                                                     >  2:2764438/1‑150 (MQ=255)
                                                                 cccATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCATTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTg                                                                                     <  2:958889/150‑1 (MQ=255)
                                                                        gTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGc                                                                                                           >  1:1179721/1‑121 (MQ=255)
                                                                        gTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGc                                                                                                           <  2:1179721/121‑1 (MQ=255)
                                                                            aCGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACCCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAg                                                                          >  1:255487/1‑150 (MQ=255)
                                                                             cGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAgg                                                                         <  1:2986561/150‑1 (MQ=255)
                                                                             cGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAgg                                                                         <  2:1314030/150‑1 (MQ=255)
                                                                             cGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAgg                                                                         <  1:817377/150‑1 (MQ=255)
                                                                               gCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGc                                                                       >  1:503766/1‑150 (MQ=255)
                                                                                   tGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGGGTTGCCGGTGGCTTCGGTCACGTGCTGGCCCTTGCCCAGCGCCAGGGCGTCa                                                                   >  2:859548/1‑150 (MQ=37)
                                                                                   tGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCa                                                                   >  1:1592548/1‑150 (MQ=255)
                                                                                     cagcaTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCAtt                                                                 >  1:1547804/1‑150 (MQ=255)
                                                                                      agcaTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTg                                                                >  2:49346/1‑150 (MQ=255)
                                                                                       gcaTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGa                                                               <  1:49346/150‑1 (MQ=255)
                                                                                       gcaTGGTGGCGACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGTTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGTCTTCGGTCACGGGCGGGCCCTGGCCCAGCGCCAGGGCGTCATTGa                                                               <  2:255487/150‑1 (MQ=255)
                                                                                         aTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                             <  1:893194/150‑1 (MQ=255)
                                                                                         aTGTTGGCCACCAGGTTGCGGTGGGCCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                             <  1:223248/150‑1 (MQ=255)
                                                                                          tGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAAc                                                            <  2:1556856/150‑1 (MQ=255)
                                                                                            ttGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCAt                                                                  >  2:501277/1‑142 (MQ=255)
                                                                                            ttGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCAt                                                                  <  1:501277/142‑1 (MQ=255)
                                                                                             tGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCa                                                         <  2:2018130/150‑1 (MQ=255)
                                                                                              ggCCACCAGGTTGCGGTGGGTCAGCATGGGGCCTGTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGCTTGGCTTCGGTAAAGCGCTTGGTCTTGTCCAGCGCCAGGGCGTCATTGAAACGCAc                                                        >  2:3023025/1‑150 (MQ=37)
                                                                                              ggCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCAc                                                        <  1:1514399/150‑1 (MQ=255)
                                                                                               gCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACg                                                       <  1:1661711/150‑1 (MQ=255)
                                                                                                 caccaGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGc                                                     >  2:1912243/1‑150 (MQ=255)
                                                                                                    caGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTg                                                  >  1:283450/1‑150 (MQ=255)
                                                                                                     aGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGc                                                 <  2:643587/150‑1 (MQ=255)
                                                                                                          gCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAg                                            <  1:1912243/150‑1 (MQ=255)
                                                                                                                 gTCAGCATGGGGCCTTTGGACACTCCGGTCGTACCGCCGGTGTACCGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTTATTGGAACGCACGGCCTGCGGGAGGTTGTAg                                     <  2:2039809/150‑1 (MQ=255)
                                                                                                                 gTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTg                                                                                     >  1:643392/1‑102 (MQ=255)
                                                                                                                 gTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTg                                                                                     <  2:643392/102‑1 (MQ=255)
                                                                                                                  tCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAgg                                    >  1:1537254/1‑150 (MQ=255)
                                                                                                                    aGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGcc                                  >  1:2331747/1‑150 (MQ=255)
                                                                                                                         ggCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCAc                             <  2:2804120/150‑1 (MQ=255)
                                                                                                                         ggCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCATTGGCCTGCGGGTTGGCTTCGGTCAAGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCAc                             >  1:2371621/1‑150 (MQ=255)
                                                                                                                            gcCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAt                          >  2:1224907/1‑150 (MQ=255)
                                                                                                                             ccTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtc                         >  2:97059/1‑150 (MQ=255)
                                                                                                                                  ggCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATcttctt                    <  2:1821303/150‑1 (MQ=255)
                                                                                                                                  cagtacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATcttctt                    >  1:543884/5‑150 (MQ=255)
                                                                                                                                   gCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCTGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTCCCGCAGGTTGTAGGCCGGCACCATCttcttc                   <  1:859548/150‑1 (MQ=255)
                                                                                                                                   agtacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAg                                            <  1:1068692/122‑1 (MQ=255)
                                                                                                                                   agtacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAg                                            >  2:1068692/4‑125 (MQ=255)
                                                                                                                                    gtacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAgc                                                                               >  2:1742805/3‑89 (MQ=39)
                                                                                                                                    gtacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAgc                                                                               <  1:1742805/87‑1 (MQ=39)
                                                                                                                                    ccACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCa                  >  1:2863556/1‑150 (MQ=255)
                                                                                                                                     cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                             <  1:2287117/106‑1 (MQ=255)
                                                                                                                                     cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCAc                 >  1:1017269/1‑150 (MQ=255)
                                                                                                                                     cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTAGTAGGCCGGCACCATCTTCTTCAc                 >  1:419437/1‑150 (MQ=255)
                                                                                                                                      acacCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACg                >  2:3035/1‑150 (MQ=255)
                                                                                                                                      acacCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACg                >  1:776595/1‑150 (MQ=255)
                                                                                                                                      acacCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACg                >  1:1630867/1‑150 (MQ=255)
                                                                                                                                         ccGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCa                              <  1:531444/133‑1 (MQ=255)
                                                                                                                                         ccGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCa                              >  2:531444/1‑133 (MQ=255)
                                                                                                                                            gtggtACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtcttc                      >  1:2460933/1‑138 (MQ=255)
                                                                                                                                            gtggtACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtcttc                      <  2:2460933/138‑1 (MQ=255)
                                                                                                                                                tACCGCCGGTGTACTGAAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCa      <  1:1875522/150‑1 (MQ=255)
                                                                                                                                                  ccgccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACg    <  2:1017269/150‑1 (MQ=255)
                                                                                                                                                    gccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                             >  1:1611411/1‑91 (MQ=255)
                                                                                                                                                    gccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                             <  2:1611411/91‑1 (MQ=255)
                                                                                                                                                    gccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtcttc                      >  2:3044743/1‑130 (MQ=255)
                                                                                                                                                    gccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtcttc                      <  1:3044743/130‑1 (MQ=255)
                                                                                                                                                    gccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACGCt  <  1:1748276/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACGCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACGCT  >  NC_002947/5172212‑5172509

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