Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,248,217 G→A 7.8% intergenic (+30/+38) argG → / ← PP_1089 argininosuccinate synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,248,2170GA 7.8% 222.9 / 7.2 77intergenic (+30/+38)argG/PP_1089argininosuccinate synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base A (3/3);  total (40/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 = 6.04e-01

CCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGG  >  NC_002947/1248075‑1248359
                                                                                                                                              |                                                                                                                                              
ccAGCTTTGAAGAAGATGGTGGTGCCTACTACCCTGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGGCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCggg                                                                                                                                          <  4:298424/149‑1 (MQ=255)
ccACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCggg                                                                                                                                          <  8:128721/149‑1 (MQ=255)
  aCCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGt                                                                                                                                        >  8:167448/1‑149 (MQ=255)
   ccTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGtt                                                                                                                                       >  2:54570/1‑149 (MQ=255)
       tGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGttttt                                                                                                                                    <  6:168476/148‑1 (MQ=255)
       tGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTAGCGTGTGCGAAAACCCCGGACCGAGCCGGGGttttc                                                                                                                                    >  3:357496/1‑147 (MQ=255)
       tGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTAGCGTGTGCGAAAACCCCGGACCGAGCCGGGGttttc                                                                                                                                    <  4:357496/148‑2 (MQ=255)
         aagaagATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTa                                                                                                                                 >  1:297648/1‑149 (MQ=255)
           gaagaTGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATg                                                                                                                               <  6:109194/149‑1 (MQ=255)
                     gtgCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCa                                                                                                                       >  6:331461/1‑147 (MQ=255)
                        ccTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGc                                                                                                                  <  2:297648/149‑1 (MQ=255)
                        ccTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCAGTCCGGGGTTTTTTTATGCCCGCGCAAAGGc                                                                                                                  <  2:210351/149‑1 (MQ=255)
                           ggttgtaGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                               >  7:317265/7‑149 (MQ=255)
                           aCAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                               >  4:180399/1‑149 (MQ=255)
                            cAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCGATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCgg                                                                                                                <  7:310516/147‑1 (MQ=255)
                            cAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                               <  6:12705/148‑1 (MQ=255)
                            cAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGAGTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                               <  4:41819/148‑1 (MQ=255)
                             aaCCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTc                                                                                                             >  3:118242/1‑149 (MQ=255)
                               ccAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGt                                                                                                           >  1:9400/1‑149 (MQ=255)
                                cAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGtt                                                                                                          <  2:180935/149‑1 (MQ=255)
                                cAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGt                                                                                                           >  3:331730/1‑148 (MQ=255)
                                       aTGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                               >  2:26283/1‑137 (MQ=255)
                                       aTGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                               <  1:26283/137‑1 (MQ=255)
                                       aTGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTc                                                                                                   >  6:80092/1‑149 (MQ=255)
                                       aTGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTc                                                                                                   <  6:27651/149‑1 (MQ=255)
                                        tgctgcTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCAGCGCAAAGGCGGGTCGTTCAACttag                                                                                                  >  7:52016/1‑147 (MQ=255)
                                          ctgctgGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGt                                                                                                <  1:22880/149‑1 (MQ=255)
                                           tgctgGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTg                                                                                               >  1:295726/1‑149 (MQ=255)
                                                gCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCt                                                                                            <  1:54570/147‑1 (MQ=255)
                                                gCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTca                                                                                          >  5:273992/1‑149 (MQ=255)
                                                 cTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTcac                                                                                         <  5:331461/149‑1 (MQ=255)
                                                 cTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTcac                                                                                         <  3:180399/149‑1 (MQ=255)
                                                  ttCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTcaca                                                                                        >  5:353435/1‑149 (MQ=255)
                                                          gCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACatcatc                                                                                   <  3:304269/146‑1 (MQ=255)
                                                          gCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACatcatc                                                                                   >  4:304269/1‑146 (MQ=255)
                                                           cTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGccc                                                                               >  8:22700/1‑149 (MQ=255)
                                                            tCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCt                                                                              >  7:199233/1‑149 (MQ=255)
                                                              aaTGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTc                                                                             >  8:235094/1‑148 (MQ=255)
                                                                  cACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGt                                                                         <  7:180372/148‑1 (MQ=255)
                                                                   aCTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgc                                                                       >  7:167670/1‑149 (MQ=255)
                                                                   aCTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgc                                                                       <  5:317298/149‑1 (MQ=255)
                                                                   aCTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgc                                                                       >  8:21334/1‑149 (MQ=255)
                                                                   aCTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgc                                                                       >  6:317298/1‑149 (MQ=255)
                                                                       cGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgcgc                                                                     <  1:229986/147‑1 (MQ=255)
                                                                          aagcccATCGCCGCCAACAACCGCCGTTCGCTGCTCTGACTGCTGGCGTGGGCGAACACCCCCGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCCAGTGCGCTGATc                                                                <  8:317265/144‑1 (MQ=255)
                                                                            gcgcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCg                                                               <  4:118242/148‑1 (MQ=255)
                                                                             cgcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGa                                                              >  4:19836/1‑148 (MQ=255)
                                                                              gcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTAt                                                                                                                                <  7:293200/81‑1 (MQ=255)
                                                                              gcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTAt                                                                                                                                >  8:293200/1‑81 (MQ=255)
                                                                              gcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGaaa                                                            <  2:327339/149‑1 (MQ=255)
                                                                              gcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGaaa                                                            <  8:54864/149‑1 (MQ=255)
                                                                               cATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGCGTGCGCGGAGCGAAAg                                                           >  8:280014/1‑149 (MQ=255)
                                                                                   gccgccAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTg                                                        <  8:171520/148‑1 (MQ=255)
                                                                                      gccAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCTAACGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGca                                                    >  3:125004/1‑149 (MQ=255)
                                                                                          acaaGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACAtt                                                >  7:50013/1‑149 (MQ=255)
                                                                                           caaGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACAtt                                                >  6:347241/1‑148 (MQ=255)
                                                                                              gggCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACg                                            >  5:179985/1‑149 (MQ=255)
                                                                                                   gTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGt                                       <  2:9400/149‑1 (MQ=255)
                                                                                                         tgctCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGAt                                 >  6:98714/1‑149 (MQ=255)
                                                                                                         tgctCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGAt                                 >  7:288423/1‑149 (MQ=255)
                                                                                                         tgctCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGAt                                 >  6:312764/1‑149 (MQ=255)
                                                                                                               gATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGAGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaa                            >  6:313311/1‑148 (MQ=255)
                                                                                                                aTTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaa                            >  7:92359/1‑147 (MQ=255)
                                                                                                                 ttGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAg                         <  8:50013/149‑1 (MQ=255)
                                                                                                                   gCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAgc                        <  3:19836/148‑1 (MQ=255)
                                                                                                                       gCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGATTCCTGTCCAGATTACAAAAGCGcc                     <  7:235094/147‑1 (MQ=255)
                                                                                                                       gCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAg                   <  8:288423/149‑1 (MQ=255)
                                                                                                                             aaaaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGaaa                                                            <  5:143552/99‑1 (MQ=38)
                                                                                                                             aaaaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGaaa                                                            >  6:143552/4‑102 (MQ=38)
                                                                                                                              aaaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGAcgcgc             >  6:247348/3‑148 (MQ=255)
                                                                                                                              aaaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGAcgcgc             <  5:247348/146‑1 (MQ=255)
                                                                                                                               gAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACgcgcgc           <  7:22700/149‑1 (MQ=255)
                                                                                                                                  aaCCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTgg        >  1:39909/1‑149 (MQ=255)
                                                                                                                                     cccGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTc     >  4:215003/1‑149 (MQ=255)
                                                                                                                                        ggCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCggg  <  6:273992/149‑1 (MQ=255)
                                                                                                                                           tCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTcc                                     >  7:93488/1‑111 (MQ=255)
                                                                                                                                           tCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTcc                                     <  8:93488/111‑1 (MQ=255)
                                                                                                                                             ggTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaa                            >  5:53263/1‑118 (MQ=255)
                                                                                                                                             ggTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaa                            <  6:53263/118‑1 (MQ=255)
                                                                                                                                             ggTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAg                         <  6:84503/121‑1 (MQ=255)
                                                                                                                                             ggTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAg                         >  5:84503/1‑121 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGG  >  NC_002947/1248075‑1248359

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: