Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,582,559 A→T 5.6% L290* (TTG→TAG) ‡ 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,5590AT 5.6% 218.6 / 2.6 72L290* (TTG→TAG) ‡treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base A (36/31);  new base T (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

TTCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCT  >  NC_002947/4582412‑4582706
                                                                                                                                                   |                                                                                                                                                      
ttCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCaa                                                                                                                                                       <  8:61673/149‑1 (MQ=255)
    acaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                                    >  2:85268/1‑148 (MQ=255)
      aGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCTAAGGCCAt                                                                                                                                                 >  3:37450/1‑149 (MQ=255)
      aGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCAt                                                                                                                                                 <  3:73889/149‑1 (MQ=255)
      aGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCAt                                                                                                                                                 <  3:332307/149‑1 (MQ=255)
       ggTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCTATg                                                                                                                                                <  7:1811/149‑1 (MQ=255)
         tAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATgga                                                                                                                                              <  2:127791/149‑3 (MQ=255)
          aGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCAATGGCCAATGCAATGTa                                                                                                                                              >  7:245421/1‑148 (MQ=255)
            tCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACAt                                                                                                                                           >  4:131853/1‑149 (MQ=255)
               cgcTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCGAGGCCATGTACATgcg                                                                                                                                        >  7:276061/1‑149 (MQ=255)
                gcTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATgcgc                                                                                                                                       <  5:256962/149‑1 (MQ=255)
                 cTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGcgcg                                                                                                                                      <  5:295732/149‑1 (MQ=255)
                    gcgGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCa                                                                                                                                   <  7:342813/149‑1 (MQ=255)
                              aCCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGa                                                                                                                           >  7:116760/1‑147 (MQ=255)
                                cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                                    <  8:139317/136‑1 (MQ=255)
                                cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                                    >  7:139317/1‑136 (MQ=255)
                                     ccAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGaa                                                                                                                   <  1:85268/148‑1 (MQ=255)
                                      cAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGaaa                                                                                                                  <  5:188634/148‑1 (MQ=255)
                                      cAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGaaa                                                                                                                  <  6:277996/148‑1 (MQ=255)
                                       aGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGc                                                                                                                <  4:328332/149‑1 (MQ=255)
                                       aGGGTCAGTTAATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGc                                                                                                                <  4:190807/149‑1 (MQ=255)
                                        ggctCAGTTCATGATGCGTGCGCAACAAGATCGCCCACTGCCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGcc                                                                                                               <  8:175424/146‑1 (MQ=255)
                                         ggTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCa                                                                                                              >  7:84935/1‑149 (MQ=255)
                                          gTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCGTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCa                                                                                                              >  5:87940/1‑148 (MQ=255)
                                             aGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATg                                                                                                          <  5:313220/149‑1 (MQ=255)
                                             aGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATg                                                                                                          <  8:116760/149‑1 (MQ=255)
                                                tcatcaTGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCa                                                                                                       >  8:35252/1‑149 (MQ=255)
                                                  atcatGGATGAGCAAAAACATCCCCCACTGGAACTTGGCCCGGATAACAGGTGTCTGGCGAAGGATTTCTCTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGAACATCCACAGCATCAACGGAAACTGGAAAGCGATATCGCAt                                                                                                      >  4:126774/1‑148 (MQ=255)
                                                            cgcAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtc                                                                                           >  6:145546/1‑149 (MQ=255)
                                                              cAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtcgc                                                                                         >  5:131088/1‑149 (MQ=255)
                                                               aaaaaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                        >  6:221986/1‑149 (MQ=255)
                                                               aaaaaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCAAGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                        >  1:175270/1‑149 (MQ=255)
                                                                  aaGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑c                                                                                  >  3:256498/1‑149 (MQ=255)
                                                                   aGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑Cg                                                                                 >  4:321285/1‑149 (MQ=255)
                                                                        gCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGaa                                                                            <  4:101472/149‑1 (MQ=255)
                                                                        gCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGaa                                                                            <  8:84935/149‑1 (MQ=255)
                                                                          ccACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGcgcg                                                                                                                                      <  4:88807/92‑1 (MQ=255)
                                                                          ccACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGcgcg                                                                                                                                      >  3:88807/1‑92 (MQ=255)
                                                                           cACTGGCAGTTGGCGGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCTGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTa                                                                         >  6:92515/1‑149 (MQ=255)
                                                                                    ttGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCgt                                                                <  7:337074/149‑1 (MQ=255)
                                                                                       gCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTc                                                             <  4:256498/149‑1 (MQ=255)
                                                                                         cGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCTTGTc                                                             >  8:64577/1‑147 (MQ=255)
                                                                                         cGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCtt                                                           >  4:205268/1‑149 (MQ=255)
                                                                                          gggATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTc                                                          <  3:321285/149‑1 (MQ=255)
                                                                                          gggATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTc                                                          <  8:276061/149‑1 (MQ=255)
                                                                                                       tCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTg                                             >  2:144552/1‑149 (MQ=255)
                                                                                                           gcgcAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCtt                                         >  8:182925/1‑149 (MQ=255)
                                                                                                            cgcAGGATGTCGGTGATAGGGAAGCTGTCTTCCATGGCTAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTc                                        <  4:37450/149‑1 (MQ=255)
                                                                                                            cgcAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTc                                        <  5:221986/149‑1 (MQ=255)
                                                                                                               aGGATGTCTGTGATAGGGAAGCGGTCTTCCATGGCAAAGGCAATGTACATGCGAGGCATCAACTGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTTGTTGGCTTCg                                       <  8:245421/147‑1 (MQ=255)
                                                                                                                 gATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGc                                     >  6:191149/1‑147 (MQ=255)
                                                                                                                  aTGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCcag                                  <  3:131853/149‑1 (MQ=255)
                                                                                                                     tCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCcagca                                <  6:131088/148‑1 (MQ=255)
                                                                                                                      cGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACgg                                                                      <  6:251616/109‑1 (MQ=255)
                                                                                                                      cGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACgg                                                                      >  5:251616/1‑109 (MQ=255)
                                                                                                                         tGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATg                           <  1:144552/149‑1 (MQ=255)
                                                                                                                               ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAATGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCg                     >  8:98032/1‑149 (MQ=255)
                                                                                                                               ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCg                     <  3:126774/149‑1 (MQ=255)
                                                                                                                               ggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCg                     >  2:130196/1‑149 (MQ=255)
                                                                                                                                  aGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCgggg                  >  5:306478/1‑149 (MQ=255)
                                                                                                                                     ggTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAg               <  5:191149/149‑1 (MQ=255)
                                                                                                                                     gcatgtacATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAg               >  4:203218/8‑149 (MQ=255)
                                                                                                                                     gcatgtacATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAg               >  2:373594/8‑149 (MQ=255)
                                                                                                                                      gTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGt              >  4:295922/1‑149 (MQ=255)
                                                                                                                                         ttCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTgg           >  5:103172/1‑149 (MQ=255)
                                                                                                                                           ccttGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTg            <  6:87940/143‑1 (MQ=255)
                                                                                                                                             aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc       >  8:110622/1‑149 (MQ=255)
                                                                                                                                               ggCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGa     >  5:184075/1‑149 (MQ=255)
                                                                                                                                                gCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGa     >  1:264243/1‑148 (MQ=255)
                                                                                                                                                 ccAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGAtc   >  8:350073/1‑149 (MQ=255)
                                                                                                                                                 ccAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGAt    >  8:146952/1‑148 (MQ=255)
                                                                                                                                                  cAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTTTGTGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGa     >  1:303497/1‑149 (MQ=255)
                                                                                                                                                  cAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGAtct  >  3:138929/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TTCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTT‑‑‑CGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCT  >  NC_002947/4582412‑4582706

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: