Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,224,632 A→G 7.3% W16R (TGG→CGG) ‡ gltK ← glutamate/aspartate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,224,6320AG 7.3% 195.8 / 5.0 68W16R (TGG→CGG) ‡gltKglutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/27);  new base G (2/3);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGCGT  >  NC_002947/1224489‑1224779
                                                                                                                                               |                                                                                                                                                   
gCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAggg                                                                                                                                                >  2:353264/1‑149 (MQ=255)
    aaGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcatg                                                                                                                                            <  3:352726/149‑3 (MQ=255)
    aaGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTCCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcatg                                                                                                                                            <  5:349418/149‑3 (MQ=255)
                   gcgcCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGgatgat                                                                                                                             <  2:325165/149‑1 (MQ=255)
                       cGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCg                                                                                                                         <  6:397541/149‑1 (MQ=255)
                           cAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcat                                                                                                                                             <  6:36820/125‑2 (MQ=255)
                           cAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcat                                                                                                                                             >  5:36820/1‑124 (MQ=255)
                             ggTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGaaa                                                                                                                   <  4:38038/149‑1 (MQ=255)
                              gTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAt                                                                                                                  >  4:6580/1‑149 (MQ=255)
                               ttCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtc                                                                                                                 >  5:167997/1‑149 (MQ=255)
                                  gACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccat                                                                                                              >  5:14362/1‑149 (MQ=255)
                                   aCAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCAGCGAAAtcatcc                                                                                                             <  4:14573/149‑5 (MQ=255)
                                    cagcagTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccattt                                                                                                            >  1:319042/1‑149 (MQ=255)
                                         gTTTGCTTTACGACAGGCGCATCAGTTCCAGGATGGTGCCCAGCACGATACCGACGATGACTCAAAAGACCATCACTTTGAGCGACATGACCATGCCCTCCCAGAGGGAACGCAGGGCGGGGATGATTTAGCTGAAAtccatttcaatt                                                                                                       >  2:289343/1‑149 (MQ=255)
                                         gTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccatttccatt                                                                                                       >  5:62303/1‑149 (MQ=255)
                                               ttGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCAGCGAAAtcatc                                                                                                              >  7:222533/1‑133 (MQ=255)
                                               ttGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCAGCGAAAtcatc                                                                                                              <  8:222533/136‑4 (MQ=255)
                                                 gacgacAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGcc                                                                                               >  3:93207/1‑149 (MQ=255)
                                                 gacgacAGGCGCATCAGTGCCAGGATGGTCCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGcc                                                                                               >  1:166856/1‑149 (MQ=255)
                                                   cgacAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGcccc                                                                                             <  4:93207/149‑1 (MQ=255)
                                                      cAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCa                                                                                          <  6:167997/149‑1 (MQ=255)
                                                       aGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCAc                                                                                         >  8:35511/1‑149 (MQ=255)
                                                        ggCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACg                                                                                        <  3:6580/149‑1 (MQ=255)
                                                         gcgcATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACgg                                                                                       <  6:62303/149‑1 (MQ=255)
                                                             aTCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGaaa                                                                                    >  6:279170/1‑148 (MQ=255)
                                                                    ccAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcca                                                                            <  2:368081/149‑1 (MQ=255)
                                                                    ccAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggc                                                                              >  5:69523/1‑147 (MQ=255)
                                                                     cAGGATGGTGCCAATCACGATACCGCCTATGACTACCATGACCATCAGTTTTAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGTAAATcaggccag                                                                           <  7:73981/149‑1 (MQ=255)
                                                                          tGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACg                                                                      >  5:51805/1‑149 (MQ=255)
                                                                               cccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGAc                                                                 >  5:326813/1‑149 (MQ=255)
                                                                                ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATca                                                                                 >  7:123272/1‑132 (MQ=255)
                                                                                ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATca                                                                                 <  8:123272/132‑1 (MQ=255)
                                                                                ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                                <  2:95378/149‑1 (MQ=255)
                                                                                ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                                >  4:331860/1‑149 (MQ=255)
                                                                                  aGCACGATACCGCCTATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACttt                                                              >  2:218445/1‑149 (MQ=255)
                                                                                            cgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAc                                                     >  2:67178/1‑148 (MQ=255)
                                                                                              ccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCAt                                                  >  6:277191/1‑149 (MQ=255)
                                                                                              ccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCAt                                                  >  2:400413/1‑149 (MQ=255)
                                                                                                    aCTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCAt                                            >  8:74588/1‑149 (MQ=255)
                                                                                                      tCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTTCTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                                                        <  1:286785/135‑1 (MQ=255)
                                                                                                      tCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTTCTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                                                        >  2:286785/1‑135 (MQ=255)
                                                                                                        ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGa                                          <  1:400413/147‑1 (MQ=255)
                                                                                                        ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATCCGCATCag                                         >  7:121198/1‑148 (MQ=255)
                                                                                                          aTGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGagcag                                      <  7:35511/149‑1 (MQ=255)
                                                                                                               cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCAACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGAAGCATCa                                 >  8:371818/1‑149 (MQ=255)
                                                                                                                 tCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAgg                               <  1:162743/149‑1 (MQ=255)
                                                                                                                      ttGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCa                           <  7:74588/148‑1 (MQ=255)
                                                                                                                           ggTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATc                                                                                  <  8:281559/88‑1 (MQ=255)
                                                                                                                           ggTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATc                                                                                  >  7:281559/1‑88 (MQ=255)
                                                                                                                             tCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAgg                   <  2:166856/149‑1 (MQ=255)
                                                                                                                             tCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAgg                   <  6:51805/149‑1 (MQ=255)
                                                                                                                                 gACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTg                >  6:64784/1‑148 (MQ=255)
                                                                                                                                  aCCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGaa              <  1:353264/149‑1 (MQ=255)
                                                                                                                                   ccATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGACCATGTTCAGGGTGAAg             >  5:397878/1‑149 (MQ=255)
                                                                                                                                    cATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGc                              >  4:22814/1‑131 (MQ=255)
                                                                                                                                    cATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGc                              <  3:22814/131‑1 (MQ=255)
                                                                                                                                    cATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAg             >  4:188183/1‑148 (MQ=255)
                                                                                                                                     ctGCCCTCTGGGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAg             >  7:225140/2‑147 (MQ=255)
                                                                                                                                      tGCCCTCCCAGAGGGCAGGCAGTGCGGCGATGAATTCGCTGAAATCCATTTCCATTTACTTGCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGa         <  8:103314/149‑1 (MQ=255)
                                                                                                                                        cccTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAg             >  5:243746/1‑144 (MQ=255)
                                                                                                                                           tcacAGAGGGCAGGCATGGCGGGGATGATTTAGCTGAAATCCATTTCCTTTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATc       <  8:121198/144‑1 (MQ=255)
                                                                                                                                           tCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGc                              >  4:167269/1‑124 (MQ=255)
                                                                                                                                           tCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGc                              <  3:167269/124‑1 (MQ=255)
                                                                                                                                           tCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGAATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGc    >  6:291862/1‑149 (MQ=255)
                                                                                                                                             ccAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGc    <  5:64784/148‑1 (MQ=255)
                                                                                                                                             ccAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGCg   >  7:246161/1‑149 (MQ=255)
                                                                                                                                             ccAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGCg   >  5:418352/1‑149 (MQ=255)
                                                                                                                                              cAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGCGt  >  3:22207/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAGTAGATCAGCGT  >  NC_002947/1224489‑1224779

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