Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,613,139 G→A 5.8% H489H (CAC→CAT)  ushA ← 5'‑nucleotidase/2',3'‑cyclic phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,613,1390GA 5.8% 209.0 / 2.2 69H489H (CAC→CAT) ushA5'‑nucleotidase/2',3'‑cyclic phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/34)
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.06e-01

CACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAATCTGCTTGCCTGT  >  NC_002947/1612992‑1613285
                                                                                                                                                   |                                                                                                                                                  
gaCGTTCTCCAACAAGGCTTCAAGGTCGGTGATGCCGGTTTCCTGGCGTCGCCTTCCTTATTTGAAAGTGGGCAAGCGACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGt                                                                                                                                                   <  2:251888/148‑1 (MQ=255)
      cTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGAc                                                                                                                                             <  6:216120/149‑1 (MQ=255)
         caccaAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGca                                                                                                                                          >  5:206858/1‑149 (MQ=255)
          accaAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGAGTCGCCTTCCATATTTGAAATTGGCCAAGCAACCGCCACCATCGGCCAGGAAAATTTTGCTGACGACGCGATAATGGGTCAACCGATCAACGGTTTAGCAGTGAAGACGca                                                                                                                                          >  6:253432/1‑148 (MQ=255)
            caAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCacac                                                                                                                                       >  8:354149/1‑149 (MQ=255)
            caAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCacac                                                                                                                                       >  4:58739/1‑149 (MQ=255)
                gCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTccc                                                                                                                                   <  2:173137/149‑1 (MQ=255)
                 cTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGc                                                                                                                                           >  1:361014/1‑140 (MQ=255)
                 cTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGc                                                                                                                                           <  2:361014/140‑1 (MQ=255)
                 cTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTcccc                                                                                                                                  <  8:203171/149‑1 (MQ=255)
                   tCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACt                                                                                                                                      >  6:82140/1‑143 (MQ=255)
                   tCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACt                                                                                                                                      <  5:82140/143‑1 (MQ=255)
                      aGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTcgac                                                                                                                             <  3:87712/149‑1 (MQ=255)
                        gTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGc                                                                                                                                           >  8:64819/1‑133 (MQ=255)
                        gTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGc                                                                                                                                           <  7:64819/133‑1 (MQ=255)
                           ggTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGacac                                                                                                                        <  3:58739/149‑1 (MQ=255)
                               aTGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCt                                                                                                                    <  5:194332/149‑1 (MQ=255)
                                  ccGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGc                                                                                                                  <  4:116053/148‑1 (MQ=255)
                                  ccGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGCTTTGCCGTGAAGACGCACACTACCCTCGACGACACGGCTGcc                                                                                                                 >  7:129211/1‑149 (MQ=255)
                                   cGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTcga                                                                                                                              <  6:310659/135‑1 (MQ=255)
                                   cGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGcc                                                                                                                 >  3:169250/1‑148 (MQ=255)
                                   cGGTGTCCTGGCCTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTcga                                                                                                                              >  5:310659/1‑135 (MQ=255)
                                        tCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACAATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGATCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAgg                                                                                                           >  5:166387/1‑149 (MQ=255)
                                             gCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGt                                                                                                       <  3:248368/148‑1 (MQ=255)
                                             gCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTg                                                                                                      <  4:169250/149‑1 (MQ=255)
                                                      tccttcTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCTACGGCTCAACGGGTTTCCCGTGAAGACGCACACTACCCTCGACGACACCGCTGCCGACAGCCCGTGAACCATcc                                                                                              >  1:108102/1‑148 (MQ=255)
                                                        cttcttTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCAc                                                                                           >  8:351030/1‑149 (MQ=255)
                                                         ttcttTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCAc                                                                                           >  1:175826/1‑148 (MQ=255)
                                                         ttcttTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACCGCTGCCGACAGGCCGTCAACCATCCCAc                                                                                           >  6:267569/1‑148 (MQ=255)
                                                         ttcttTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCATGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATc                                                                                               >  6:241743/1‑144 (MQ=255)
                                                         ttcttTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCATGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATc                                                                                               <  5:241743/144‑1 (MQ=255)
                                                           cttTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCAc                                                                                           <  2:175826/146‑1 (MQ=255)
                                                              tGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATaa                                                                                     >  1:348141/1‑149 (MQ=255)
                                                                aaaGTGGCCAAGCCACCGCCACCATCGGCCAGGAAATGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAAt                                                                                    <  6:206858/148‑1 (MQ=255)
                                                                        cAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGACCCATCCCACTCATAATGAAAACCCt                                                                           >  1:74032/1‑149 (MQ=255)
                                                                        cAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGATCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCt                                                                           <  5:253432/149‑1 (MQ=255)
                                                                        cAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGCTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGCCTGCCGACAGCCAGTGAACCATCCCACTCATAATGAAAAACCt                                                                           >  4:180330/1‑149 (MQ=255)
                                                                         aaGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCtt                                                                          <  7:99526/149‑1 (MQ=255)
                                                                         aaGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCATCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCtt                                                                          <  5:267569/149‑1 (MQ=255)
                                                                            ccaccgccaccATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGa                                                                       <  7:354149/149‑1 (MQ=255)
                                                                               ccgccaccATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCATGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATAc                                                                    >  5:108767/1‑149 (MQ=255)
                                                                                           gCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGTCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGAtt                                                        <  3:180330/149‑1 (MQ=255)
                                                                                            ccAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGAACACTCCCCTCGAAGACACGGCTCCCGACAGGCCGTGAACCATCCCACTCATAATGAAACCCCTTCGAAACTTTCAGCGGAtt                                                        >  5:174159/1‑148 (MQ=255)
                                                                                              aGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGCCTGCCGACAGGCCGAGAACCAACCCACTCATAATGAAAACCCATCGATACATGCAGCGGATTGa                                                      >  2:118017/1‑148 (MQ=255)
                                                                                                       ttGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCt                                            >  6:338668/1‑149 (MQ=255)
                                                                                                           tGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTc                                        >  2:142119/1‑149 (MQ=255)
                                                                                                               gacgCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGGCACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTc                                                                         >  4:133942/1‑112 (MQ=255)
                                                                                                               gacgCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGGCACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTc                                                                         <  3:133942/112‑1 (MQ=255)
                                                                                                                 cgcgATAATTGGTCAACGGCTCAACGGGTTTGCCGTGAAGACTCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGa                                  <  7:61800/149‑1 (MQ=255)
                                                                                                                  gcgATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTg                                                                 <  7:232407/117‑1 (MQ=255)
                                                                                                                  gcgATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTg                                                                 >  8:232407/1‑117 (MQ=255)
                                                                                                                    gATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCa                                                                                        >  4:227247/1‑92 (MQ=255)
                                                                                                                    gATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCa                                                                                        <  3:227247/92‑1 (MQ=255)
                                                                                                                      tAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGtt                             <  1:142119/149‑1 (MQ=255)
                                                                                                                         tggtgcAACGGCTCAACGGGTTTGCCGGGAAGACGCACACTACCCTCGCCGACACCGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAg                          <  8:129211/144‑1 (MQ=255)
                                                                                                                         tGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAg                          <  7:351030/149‑1 (MQ=255)
                                                                                                                         tGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAg                          <  5:338668/149‑1 (MQ=255)
                                                                                                                             tCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGtt                             >  6:350134/1‑142 (MQ=255)
                                                                                                                             tCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGtt                             <  5:350134/142‑1 (MQ=255)
                                                                                                                             tCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGt                      >  8:258428/1‑149 (MQ=255)
                                                                                                                              cAACGGCTCAACGGGTTTGCCGTGAAGACGCAGACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAAGTGCCCCTGTTCCATTGAAGGTTGAGCAGt                      <  1:118017/148‑1 (MQ=255)
                                                                                                                               aaCGGCTCATCGGGTTTGCCGTGAAGACGCAGACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCTTCCCACTCATAATGTAAACCCTTCGATACTTGCAGCGGATTGAAACTGGCCCTGTTCCATTGAAGGTTGAGCAGtt                     <  6:249830/148‑1 (MQ=255)
                                                                                                                                aCGGCTCAACGGGTTTGCCATGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCg                   <  6:108767/149‑1 (MQ=255)
                                                                                                                                 cGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGTAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCgg                  >  2:27219/1‑149 (MQ=255)
                                                                                                                                     tCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAAt              >  6:12758/1‑149 (MQ=255)
                                                                                                                                         cGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAATCTGc          >  4:46432/1‑149 (MQ=255)
                                                                                                                                          gggTTTGCCGTGAACACGAACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCAATCATAATGAAAACCATTAGATACTTGCAGCGCATTGAAACTGCCCCTGTTAAAATGAAGGTTGACAAATTCGGTAAAATGCt         >  8:233050/1‑149 (MQ=255)
                                                                                                                                             tttGCCGTGAAGAAGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAATCTGCTTGc      >  2:190558/1‑149 (MQ=255)
                                                                                                                                                gCCGTGAAGCCGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAATCTGCTTGCCTg   >  3:247867/1‑149 (MQ=255)
                                                                                                                                                   gTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAATCTGCTTGCCTGt  >  4:51581/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCTGACGACGCGATAATGGGTCAACGGCTCAACGGGTTTGCCGTGAAGACGCACACTCCCCTCGACGACACGGCTGCCGACAGGCCGTGAACCATCCCACTCATAATGAAAACCCTTCGATACTTGCAGCGGATTGAAACTGCCCCTGTTCCATTGAAGGTTGAGCAGTTCGGTAATCTGCTTGCCTGT  >  NC_002947/1612992‑1613285

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