Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,582,560 A→T 6.1% L290M (TTG→ATG)  treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,582,5600AT 6.1% 245.6 / 4.0 82L290M (TTG→ATG) treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base A (39/38);  new base T (3/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTC  >  NC_002947/4582414‑4582701
                                                                                                                                                  |                                                                                                                                             
ccACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAgg                                                                                                                                             >  4:20130/1‑149 (MQ=255)
ccACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAg                                                                                                                                              <  1:135293/148‑1 (MQ=255)
 cacaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAgg                                                                                                                                             <  8:15552/148‑1 (MQ=255)
 cacaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                            >  8:159882/1‑149 (MQ=255)
 cacaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                            >  1:297368/1‑149 (MQ=255)
 cacaGGTAGTCGCGCTCGCGGTCGGGGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAgg                                                                                                                                             <  1:131786/148‑1 (MQ=255)
  acaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                            >  6:177581/1‑148 (MQ=255)
   caGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCa                                                                                                                                          <  1:10239/149‑1 (MQ=255)
     ggTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATg                                                                                                                                        <  3:20130/149‑1 (MQ=255)
      gTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATg                                                                                                                                        >  3:288210/1‑148 (MQ=255)
       tAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATgga                                                                                                                                      <  4:123092/149‑3 (MQ=255)
       tAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTa                                                                                                                                      <  2:58407/149‑1 (MQ=255)
       tAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTa                                                                                                                                      <  7:159882/149‑1 (MQ=255)
          tCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATggaaga                                                                                                                                   <  3:259384/149‑6 (MQ=255)
            gcgcTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATgc                                                                                                                                 >  3:176950/1‑149 (MQ=255)
                tCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATg                                                                                                                                        <  4:271699/138‑1 (MQ=255)
                tCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATg                                                                                                                                        >  3:271699/1‑138 (MQ=255)
                tCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCgg                                                                                                                             >  2:325742/1‑149 (MQ=255)
                     gtcggtGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATc                                                                                                                         >  7:237853/1‑148 (MQ=255)
                        ggtGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAAc                                                                                                                      <  4:18218/148‑1 (MQ=255)
                            aCCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACg                                                                                                                     <  5:175130/145‑1 (MQ=255)
                            aCCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACg                                                                                                                     >  6:175130/1‑145 (MQ=255)
                            aCCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGa                                                                                                                   >  1:145315/1‑147 (MQ=255)
                              cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGt                                                                                                               <  2:297368/149‑1 (MQ=255)
                                tttCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTgg                                                                                                             <  2:172071/149‑1 (MQ=255)
                                tttCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTg                                                                                                              >  1:40085/1‑148 (MQ=255)
                                   ccAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                            <  1:126297/131‑1 (MQ=255)
                                   ccAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                            >  2:126297/1‑131 (MQ=255)
                                   ccAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGaaa                                                                                                                 >  4:26367/1‑142 (MQ=255)
                                   ccAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGaaa                                                                                                                 <  3:26367/142‑1 (MQ=255)
                                     aGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACAt                                                                                                                                   <  2:328688/122‑1 (MQ=255)
                                     aGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACAt                                                                                                                                   >  1:328688/1‑122 (MQ=255)
                                     aGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGc                                                                                                        >  4:239490/1‑149 (MQ=255)
                                            gTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATgg                                                                                                 <  4:288210/149‑1 (MQ=255)
                                                atcatGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTAc                                                                                                                                     <  7:22900/109‑1 (MQ=255)
                                                atcatGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTAc                                                                                                                                     <  7:27954/109‑1 (MQ=255)
                                                atcatGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTAc                                                                                                                                     >  8:22900/1‑109 (MQ=255)
                                                atcatGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTAc                                                                                                                                     >  8:27954/1‑109 (MQ=255)
                                                  catGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATtcg                                                                                           >  1:163333/1‑149 (MQ=255)
                                                             aaaaaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                >  5:216508/1‑149 (MQ=255)
                                                             aaaaaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                >  2:210209/1‑149 (MQ=255)
                                                             aaaaaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                >  5:313637/1‑149 (MQ=255)
                                                                aaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTc                                                                             <  3:239490/149‑1 (MQ=255)
                                                                aaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTc                                                                             <  3:144677/149‑1 (MQ=255)
                                                                    tCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGTGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGcc                                                                          >  3:280371/1‑148 (MQ=255)
                                                                        ccACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGt                                                                     <  4:176950/149‑1 (MQ=255)
                                                                         cACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTa                                                                    >  1:266002/1‑149 (MQ=255)
                                                                           cTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACg                                                                  <  1:145058/149‑1 (MQ=255)
                                                                              gCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACgg                                                                 <  2:163333/147‑1 (MQ=255)
                                                                                    ggCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCgtgt                                                         >  6:166351/1‑149 (MQ=255)
                                                                                    ggCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCgtg                                                          <  2:145315/148‑1 (MQ=255)
                                                                                     gCCGGGCTGTCCGGGGGCGGGCGCAGGATGTCGGTGATAGGGACGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCgtgt                                                         <  8:237853/148‑1 (MQ=255)
                                                                                          gATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGTGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAg                                                   <  4:280371/149‑1 (MQ=255)
                                                                                          gATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGcc                                                                                                                                           <  6:322430/61‑1 (MQ=255)
                                                                                          gATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGcc                                                                                                                                           >  5:322430/1‑61 (MQ=255)
                                                                                          gATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAg                                                   <  2:266002/149‑1 (MQ=255)
                                                                                             ttCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACgg                                                                 <  8:110194/132‑1 (MQ=255)
                                                                                             ttCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACgg                                                                 >  7:110194/1‑132 (MQ=255)
                                                                                             ttCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGc                                                 >  1:191670/1‑148 (MQ=255)
                                                                                             ttCCGAGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGc                                                 >  6:91491/1‑148 (MQ=255)
                                                                                                 ggggTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCTTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTc                                                     <  7:280002/140‑1 (MQ=255)
                                                                                                 ggggTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCATCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTc                                                     >  8:280002/1‑140 (MQ=255)
                                                                                                     tCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTg                                        >  1:18473/1‑149 (MQ=255)
                                                                                                           gcAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCg                                  <  6:216508/149‑1 (MQ=255)
                                                                                                                     ggTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCgg                                 <  6:207140/140‑1 (MQ=255)
                                                                                                                     ggTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCgg                                 >  5:207140/1‑140 (MQ=255)
                                                                                                                       tGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATg                      <  1:325742/149‑1 (MQ=255)
                                                                                                                         aTAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGc                     >  2:241239/1‑148 (MQ=255)
                                                                                                                             ggAAGCGGTCTTCCATGGCCAAGGCCATGTTCATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCg                >  7:77641/1‑149 (MQ=255)
                                                                                                                             ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATg                      <  5:39002/143‑1 (MQ=255)
                                                                                                                             ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATg                      >  6:39002/1‑143 (MQ=255)
                                                                                                                             ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATc                 >  1:316629/1‑148 (MQ=255)
                                                                                                                             ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATc                 <  2:316629/148‑1 (MQ=255)
                                                                                                                             ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCg                <  2:40085/149‑1 (MQ=255)
                                                                                                                                    catgtacATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGt         >  4:323393/7‑149 (MQ=255)
                                                                                                                                     atgtacATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGtt        >  6:273452/6‑149 (MQ=255)
                                                                                                                                       ttCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTgg      >  4:228192/1‑149 (MQ=255)
                                                                                                                                       ttCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTgg      >  8:273682/1‑149 (MQ=255)
                                                                                                                                         ccATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCGTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGc     <  6:313637/148‑1 (MQ=255)
                                                                                                                                         acATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTCGGCg    >  3:291181/2‑149 (MQ=255)
                                                                                                                                          cATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGt   >  5:67666/1‑149 (MQ=255)
                                                                                                                                           aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc  <  2:18473/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTC  >  NC_002947/4582414‑4582701

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: