Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 863,903 T→A 6.0% intergenic (+17/+19) hemH → / ← uraA ferrochelatase/uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,9030TA 6.0% 415.8 / 14.4 134intergenic (+17/+19)hemH/uraAferrochelatase/uracil permease
Reads supporting (aligned to +/- strand):  ref base T (66/60);  new base A (3/5);  total (69/65)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 138 total aligned reads displayed.

GCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATC  >  NC_002947/863754‑864051
                                                                                                                                                     |                                                                                                                                                      
 cTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGAGGGTGAGTGTGTTGGCTGATGTGTGTGGAAAGGCCTGAATCGAAGGTGCCGCTTTg                                                                                                                                                       >  2:1063524/1‑150 (MQ=255)
  tGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGc                                                                                                                                                      <  1:475327/150‑1 (MQ=255)
   ggAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGCGGAGCCATGCGGCCAGGAGCTGGTGCTGGTGCCGTGTCGGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCg                                                                                                                                                     <  1:2214075/150‑1 (MQ=255)
    gAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGGGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCgg                                                                                                                                                    >  2:2093351/1‑150 (MQ=255)
      agagATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGcc                                                                                                                                                  <  2:1681464/150‑1 (MQ=255)
      agagATAGGCATGCGTGGCAGCAAGAGCGTGGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTTGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGcc                                                                                                                                                  <  1:2307706/150‑1 (MQ=255)
         gaTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCt                                                                                                                                               <  2:2919502/150‑1 (MQ=255)
          aTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCtt                                                                                                                                              <  1:2207442/150‑1 (MQ=255)
              gCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGa                                                                                                                                          >  2:2478603/1‑150 (MQ=255)
                 tGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCg                                                                                                                                       >  2:1546168/1‑150 (MQ=255)
                  gCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCgg                                                                                                                                      <  1:1546168/150‑1 (MQ=255)
                   cGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGGGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGt                                                                                                                                     <  1:1063524/150‑1 (MQ=255)
                    gTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTc                                                                                                                                    >  1:421556/1‑150 (MQ=255)
                       gCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGACGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATa                                                                                                                                 <  2:862247/150‑1 (MQ=255)
                         aGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTc                                                                                                                                             <  2:382891/136‑1 (MQ=255)
                          gCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTcga                                                                                                                                           <  2:2559727/137‑3 (MQ=255)
                          gCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTcga                                                                                                                                           >  1:2559727/1‑135 (MQ=255)
                           cGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATc                                                                                                                             <  1:946219/150‑1 (MQ=255)
                            gAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATcc                                                                                                                            >  1:3019224/1‑150 (MQ=255)
                                      tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgt                                                                                                                  <  1:2594608/150‑1 (MQ=255)
                                         aGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTc                                                                                                                                             <  1:2720773/120‑1 (MQ=255)
                                             gggcggCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCa                                                                                                           >  1:1637460/1‑150 (MQ=255)
                                               gcggcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGc                                                                                                         <  2:2964036/150‑1 (MQ=255)
                                                cggcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCt                                                                                                        <  1:2478603/150‑1 (MQ=255)
                                                 ggcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTg                                                                                                       <  1:1842321/150‑1 (MQ=255)
                                                   ccAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTc                                                                                                                                    <  1:2210240/119‑1 (MQ=255)
                                                   ccAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTc                                                                                                                                    >  2:2210240/1‑119 (MQ=255)
                                                       gAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTa                                                                                                 >  1:2294026/1‑150 (MQ=255)
                                                          ctggtgctggtgCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTc                                                                                                             <  2:417526/135‑1 (MQ=255)
                                                          ctggtgctggtgCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTc                                                                                                             >  1:417526/1‑135 (MQ=255)
                                                             gtgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACGGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGc                                                                                          >  2:534587/1‑150 (MQ=255)
                                                               gctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAg                                                                                        >  1:573374/1‑150 (MQ=255)
                                                                 tggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGaa                                                                                      >  2:2640398/1‑150 (MQ=255)
                                                                   gtgccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATc                                                                                    >  1:2639892/1‑150 (MQ=255)
                                                                     gccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg                                                                                  <  2:2599981/150‑1 (MQ=255)
                                                                      ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAgg                                                                                 <  2:1662904/150‑1 (MQ=255)
                                                                      ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAgg                                                                                 <  1:1819624/150‑1 (MQ=255)
                                                                        gtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCCTACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGtt                                                                               >  1:2107835/1‑150 (MQ=255)
                                                                           ccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGtt                                                                               >  1:1054166/1‑147 (MQ=255)
                                                                           ccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAg                                                                            >  1:2426580/1‑150 (MQ=255)
                                                                            cTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAgg                                                                           >  1:1947299/1‑150 (MQ=255)
                                                                              gAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACa                                                                                                                               >  2:2330519/1‑97 (MQ=255)
                                                                              gAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCa                                                                                                           >  1:85653/1‑117 (MQ=255)
                                                                              gAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCCAGTTAAGGGc                                                                         >  1:3030054/1‑150 (MQ=255)
                                                                               aaTGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCa                                                                        >  1:882544/1‑150 (MQ=255)
                                                                                accctCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCa                                                                                                           <  1:77299/111‑1 (MQ=255)
                                                                                accctCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCa                                                                                                           >  2:77299/5‑115 (MQ=255)
                                                                                  gATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCa                                                                             <  2:2454351/142‑1 (MQ=255)
                                                                                  gATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCa                                                                             >  1:2454351/1‑142 (MQ=255)
                                                                                  gATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATa                                                                     >  2:971011/1‑150 (MQ=255)
                                                                                    tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATcc                                                                                                                            <  2:51337/94‑1 (MQ=255)
                                                                                     cACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGcc                                                                  >  1:1657739/1‑150 (MQ=255)
                                                                                       cccGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCcac                                                                >  1:3048073/1‑150 (MQ=255)
                                                                                              tGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGgcg                                                         <  2:573374/150‑1 (MQ=255)
                                                                                                ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGca                                                       >  1:2032412/1‑150 (MQ=255)
                                                                                                ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCACCTTCCGACCGGTCATACATCCAGCTTCTTGTGCt                                                                                                               >  1:1575957/1‑95 (MQ=255)
                                                                                                ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCACCTTCCGACCGGTCATACATCCAGCTTCTTGTGCt                                                                                                               <  2:1575957/95‑1 (MQ=255)
                                                                                                 gTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTg                                                                                                       <  2:819523/102‑1 (MQ=255)
                                                                                                 gTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTg                                                                                                       >  1:819523/1‑102 (MQ=255)
                                                                                                 gTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGcac                                                      >  2:2677509/1‑150 (MQ=255)
                                                                                                 gTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGcac                                                      >  1:568956/1‑150 (MQ=255)
                                                                                                       gtgtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCCGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCg                                               <  2:1637460/150‑1 (MQ=255)
                                                                                                         gtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGAGACAATGGCGGGCAGGGCGac                                             >  2:1709489/1‑149 (MQ=255)
                                                                                                         gtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTACGCAGAAGAGCCACAATGGCGCACAGGGCGAt                                             <  1:2093351/150‑1 (MQ=255)
                                                                                                          ttGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATg                                            <  1:1282519/150‑1 (MQ=255)
                                                                                                           tGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCt                                                                                                        >  2:2693846/1‑91 (MQ=255)
                                                                                                             gCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGccc                                         <  2:1947299/150‑1 (MQ=255)
                                                                                                             gCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGccc                                         >  1:1214715/1‑150 (MQ=255)
                                                                                                             gCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGccc                                         <  1:547884/150‑1 (MQ=255)
                                                                                                                gATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTc                                      <  2:2294026/150‑1 (MQ=255)
                                                                                                                gATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTc                                      >  2:1537179/1‑150 (MQ=255)
                                                                                                                  tatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATc                                                                                                                             >  1:620923/1‑63 (MQ=255)
                                                                                                                  tatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATc                                                                                                                             <  2:620923/63‑1 (MQ=255)
                                                                                                                   atGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAgg                                   <  2:3048073/150‑1 (MQ=255)
                                                                                                                   atGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGATTCCAGCTGTCATTACCCGGCAGAATAAGGTTCAGGGCAATAGCCACAATAGCGCACAGGGCGATGCCCTTCAgg                                   >  2:3056905/1‑150 (MQ=255)
                                                                                                                      tgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGtt                                                                               <  1:1915025/104‑1 (MQ=255)
                                                                                                                      tgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGtt                                                                               >  2:1915025/1‑104 (MQ=255)
                                                                                                                      tgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGccc                                <  2:2426580/150‑1 (MQ=255)
                                                                                                                        tgAAAAGGCCTGCATCGACGCGGCCGCTTTGCGGCCCCTTCCGACCTGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                              <  2:3030054/150‑1 (MQ=255)
                                                                                                                        tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg                                                                                  >  1:528413/1‑99 (MQ=255)
                                                                                                                        tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg                                                                                  <  2:528413/99‑1 (MQ=255)
                                                                                                                        tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                              <  1:1537179/150‑1 (MQ=255)
                                                                                                                        tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTGGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATa                                                                     >  1:1049687/1‑112 (MQ=255)
                                                                                                                        tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTGGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATa                                                                     <  2:1049687/112‑1 (MQ=255)
                                                                                                                           aaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATc                                                                                                                             >  2:572789/1‑54 (MQ=255)
                                                                                                                                ccTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGccc                                <  1:891561/140‑1 (MQ=255)
                                                                                                                                ccTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGccc                                >  2:891561/1‑140 (MQ=255)
                                                                                                                                ccTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCg                      <  2:882544/150‑1 (MQ=255)
                                                                                                                                      tCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCGGTCATTACCCGGCAGAATCAGGTTTAGGGCAAGAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGt                >  2:2765874/1‑150 (MQ=255)
                                                                                                                                        aactgtgCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCAAAATGGCAAACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCa              >  1:1361628/7‑150 (MQ=255)
                                                                                                                                         aaGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAggg                                                                          <  2:1102927/90‑1 (MQ=255)
                                                                                                                                         aaGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAc             >  1:1775416/1‑150 (MQ=255)
                                                                                                                                            gggCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGc          >  1:2839181/1‑150 (MQ=255)
                                                                                                                                             ggCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGTCGAGGTCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCt         >  2:1781523/1‑150 (MQ=255)
                                                                                                                                              gCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATg                                                            <  1:972340/99‑1 (MQ=255)
                                                                                                                                              gCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATg                                                            >  2:972340/1‑99 (MQ=255)
                                                                                                                                               ccGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCCTTACACCGACAAATCACGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGc       <  1:2672995/150‑1 (MQ=255)
                                                                                                                                                   tttGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGAt   >  2:1170503/1‑150 (MQ=255)
                                                                                                                                                   tttGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGAt   <  1:2640398/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATc  <  2:2032412/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATC  >  NC_002947/863754‑864051

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