Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,894,573 A→G 5.5% V258A (GTG→GCG)  cysW ← sulfate ABC transporter permease CysW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,894,5730AG 5.5% 214.6 / 2.3 73V258A (GTG→GCG) cysWsulfate ABC transporter permease CysW
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GAAGCTGTTGAAGCGCTTGCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAGCACCACGCCGT  >  NC_002947/5894425‑5894720
                                                                                                                                                    |                                                                                                                                                   
gAAGCTGTTGATGCGCTTGCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCa                                                                                                                                                     <  8:108702/149‑1 (MQ=255)
  aGCTGTTGAAGCGCTTGCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACg                                                                                                                                                   <  3:382140/149‑1 (MQ=255)
   gCTGTTGAAGCGCTTGCTGACGTTACGATCTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGc                                                                                                                                                  >  1:225090/1‑149 (MQ=255)
       ttGAAGCGCTTGCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGaa                                                                                                                                              <  7:13701/149‑1 (MQ=255)
              gcTTGCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCcgc                                                                                                                                        >  8:290757/1‑148 (MQ=255)
                  gCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGt                                                                                                                                   >  7:314785/1‑149 (MQ=255)
                   cTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTg                                                                                                                                  <  7:45868/149‑1 (MQ=255)
                    tGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTg                                                                                                                                                >  6:227337/1‑134 (MQ=255)
                    tGACGTTACGAACCTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTg                                                                                                                                                <  5:227337/134‑1 (MQ=255)
                      actgaaCGAACATCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGtt                                                                                                                               <  8:315141/144‑1 (MQ=255)
                      aCGTTACTAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGtt                                                                                                                               <  1:135261/149‑1 (MQ=255)
                      aCGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTgg                                                                                                                                 <  6:183151/147‑1 (MQ=255)
                      aCGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGt                                                                                                                                <  8:26893/148‑1 (MQ=255)
                        gTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACATGCTGGCCACGATGAAGGACGCGTCGTGGTTg                                                                                                                              >  4:41882/1‑148 (MQ=255)
                        gTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGa                                                                                                                                      <  6:223301/140‑1 (MQ=255)
                        gTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGa                                                                                                                                      >  5:223301/1‑140 (MQ=255)
                           aCGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTAc                                                                                                                           <  3:265018/148‑1 (MQ=255)
                           aCGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACt                                                                                                                          >  5:367484/1‑149 (MQ=255)
                            cGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTc                                                                                                                         <  5:353908/149‑1 (MQ=255)
                             gAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCg                                                                                                                        >  5:33873/1‑149 (MQ=255)
                             gAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCg                                                                                                                        >  4:82172/1‑149 (MQ=255)
                             gAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCg                                                                                                                        >  1:328434/1‑149 (MQ=255)
                               aCTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGtt                                                                                                                      <  3:82172/149‑1 (MQ=255)
                               aCTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGtt                                                                                                                      <  4:171037/149‑1 (MQ=255)
                                cTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTg                                                                                                                     <  7:290757/149‑1 (MQ=255)
                                 ttCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTg                                                                                                                     >  5:167410/1‑148 (MQ=255)
                                     atcgaCATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTagag                                                                                                                <  1:235169/149‑1 (MQ=255)
                                            tGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTcc                                                                                                         <  1:285434/149‑1 (MQ=255)
                                                  aTTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATg                                                                                                    >  3:88438/1‑148 (MQ=255)
                                                   ttCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATg                                                                                                    >  8:89750/1‑147 (MQ=255)
                                                    tcctccGCGGCGCTGTGGCGCAGGTGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTAc                                                                                                                           <  7:226919/123‑1 (MQ=255)
                                                    tcctccGCGGCGCTGTGGCGCAGGTGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTAc                                                                                                                           >  8:226919/1‑123 (MQ=255)
                                                     cctccGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAg                                                                                                 <  6:339204/148‑1 (MQ=255)
                                                      ctccGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCg                                                                                               >  7:46405/1‑149 (MQ=255)
                                                       tccGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCgg                                                                                              >  3:267845/1‑149 (MQ=255)
                                                             gcgcTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGgt                                                                                        <  2:328434/149‑1 (MQ=255)
                                                               gcTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGgtgt                                                                                      <  6:367484/149‑1 (MQ=255)
                                                                   tgGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTg                                                                                    <  6:68178/147‑1 (MQ=255)
                                                                        cAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCGCGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGc                                                                             >  7:239838/1‑149 (MQ=255)
                                                                            cGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACg                                                                         >  8:67426/1‑149 (MQ=255)
                                                                             ggTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACg                                                                         >  5:8671/1‑148 (MQ=255)
                                                                               ttAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCGCGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGAt                                                                      <  8:239838/149‑1 (MQ=255)
                                                                                          ctcGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGac                                                            <  6:33873/148‑1 (MQ=255)
                                                                                             gCTCCACTGCTTGAGCAGGAGGATGAAGAGCGCCAGGATCAGCAACATGCTGCCCACGCTGAAGGCCGCGTCGTGGTTGTACTCGTTGTAGAGCATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGAcacca                                                        <  3:41882/149‑1 (MQ=255)
                                                                                             gCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGAcacca                                                        >  4:50136/1‑149 (MQ=255)
                                                                                               tCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACaccacc                                                      >  1:13884/1‑149 (MQ=255)
                                                                                                 cACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCGCGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGccc                                                               <  1:89748/138‑1 (MQ=255)
                                                                                                 cACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCGCGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGccc                                                               >  2:89748/1‑138 (MQ=255)
                                                                                                 cACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGa                                                    >  2:112105/1‑149 (MQ=255)
                                                                                                  aCTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTg                                                                                                                     <  8:9953/83‑1 (MQ=255)
                                                                                                  aCTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTg                                                                                                                     >  7:9953/1‑83 (MQ=255)
                                                                                                  aCTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGa                                                    >  7:265203/1‑148 (MQ=255)
                                                                                                         cagcagcAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACcgc                                              <  6:12997/146‑1 (MQ=255)
                                                                                                             agcagGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCaaa                                        >  3:345789/1‑149 (MQ=255)
                                                                                                                  gATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGc                                   >  7:118714/1‑149 (MQ=255)
                                                                                                                  gATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGc                                   <  4:345789/149‑1 (MQ=255)
                                                                                                                     gAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCGGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCa                                >  2:189973/1‑149 (MQ=255)
                                                                                                                       agagCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATg                              <  8:314785/149‑1 (MQ=255)
                                                                                                                       agagCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATg                              >  3:290985/1‑149 (MQ=255)
                                                                                                                       agagCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATg                              <  8:118714/149‑1 (MQ=255)
                                                                                                                         agCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGc                            <  6:8671/149‑1 (MQ=255)
                                                                                                                               aGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCcgcgc                      <  4:267845/149‑1 (MQ=255)
                                                                                                                               aGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCcgcgc                      <  4:290985/149‑1 (MQ=255)
                                                                                                                                 gATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACTTGGTTGTAGTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCgg                    <  8:46405/149‑1 (MQ=255)
                                                                                                                                 gATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCgg                    <  2:225090/149‑1 (MQ=255)
                                                                                                                                   tCAGCAACCGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTa                  >  5:337690/1‑149 (MQ=255)
                                                                                                                                     aGCAACAGGCTGGCCACGCTGAGGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTaca                >  7:247930/1‑149 (MQ=255)
                                                                                                                                     aGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTac                 <  7:122712/148‑1 (MQ=255)
                                                                                                                                      gCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTac                 >  8:122712/1‑147 (MQ=255)
                                                                                                                                        aaCAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAg             >  3:60538/1‑149 (MQ=255)
                                                                                                                                        aaCAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAg             <  6:337690/149‑1 (MQ=255)
                                                                                                                                           aGGCTGGCCACGCTGAAGTCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAGcac          <  4:88438/149‑1 (MQ=255)
                                                                                                                                                  ccACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAGCACCACGCCg   >  6:398518/1‑149 (MQ=255)
                                                                                                                                                   cACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAGCACCACGCCGt  <  6:288972/149‑1 (MQ=255)
                                                                                                                                                    aCGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAGCACCACGCCGt  >  5:152016/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGCTGTTGAAGCGCTTGCTGACGTTACGAACTTCGATCGACATGAATCATTCCTCCGCGGCGCTGTGGCGCAGGCGGTTAATACGGTTCTCGCTCCACTGCTTGAGCAGCAGGATGAAGAGCGCCAGGATCAGCAACAGGCTGGCCACGCTGAAGGCCGCGACGTGGTTGTACTCGTTGTAGAGGATCTCCACATGCAGCGGCAAGGTGTTGGTCACGCCACGGATGTGGCCCGACACCACCGACACCGCGCCAAACTCGCCCATGGCCCGCGCGGTACACAGCACCACGCCGT  >  NC_002947/5894425‑5894720

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: