Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,861 A→T 8.1% L147L (CTT→CTA) ‡ glmS ← L‑glutamine/D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,171,8610AT 8.1% 182.1 / 5.3 62L147L (CTT→CTA) ‡glmSL‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (2/3);  total (29/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 = 7.28e-01

GTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTTGTG  >  NC_002947/6171720‑6172007
                                                                                                                                              |                                                                                                                                                  
gTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTTGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                                                                                                  <  7:109511/145‑1 (MQ=255)
gTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTaa                                                                                                                                                   <  1:148808/143‑1 (MQ=255)
gTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTaa                                                                                                                                                   >  2:148808/1‑143 (MQ=255)
gTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGt                                                                                                                                             >  2:21814/1‑149 (MQ=255)
 tCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTc                                                                                                                                            <  1:49628/149‑1 (MQ=255)
  cGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGTCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAggcggcg                                                                                                                                            <  1:81141/149‑6 (MQ=255)
     aCGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTc                                                                                                                                        <  5:341127/149‑1 (MQ=255)
      cGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCa                                                                                                                                       <  4:141228/149‑1 (MQ=255)
       gCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAg                                                                                                                                      <  7:27508/149‑1 (MQ=255)
       gCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGCGGCGGtgaa                                                                                                                                      <  2:184053/149‑4 (MQ=255)
       gCCAGGAAGTTCTCGCCAAGGCCCAGGCCGAAAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCGGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAg                                                                                                                                      <  2:294157/149‑1 (MQ=255)
        ccAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGCGGCGGtgaag                                                                                                                                     <  8:174905/149‑5 (MQ=255)
         cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGt                                                                                                                                    >  6:184702/1‑149 (MQ=255)
         cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTc                                                                                                                                   >  2:67951/1‑149 (MQ=255)
                ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAAt                                                                                                                             >  3:265693/1‑149 (MQ=255)
                ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGcggcgg                                                                                                                                          <  8:233053/136‑1 (MQ=255)
                ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGcggcgg                                                                                                                                          >  7:233053/1‑136 (MQ=255)
                     gCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctct                                                                                                                        >  8:17433/1‑149 (MQ=255)
                      ccAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGCTCAGGTCTGGAATGCTCtt                                                                                                                       >  8:8554/1‑149 (MQ=255)
                       cAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCt                                                                                                                                  >  5:330328/1‑137 (MQ=255)
                       cAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCt                                                                                                                                  <  6:330328/137‑1 (MQ=255)
                              cAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGgtgt                                                                                                               >  8:96932/1‑149 (MQ=255)
                                 gCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                                            >  8:46042/1‑149 (MQ=255)
                                     aTAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGAt                                                                                                        <  1:21814/149‑1 (MQ=255)
                                         ccAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                                                    >  4:156961/1‑149 (MQ=255)
                                           aGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCg                                                                                                                                           <  4:177838/108‑1 (MQ=255)
                                           aGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCg                                                                                                                                           >  3:177838/1‑108 (MQ=255)
                                              gggCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTg                                                                                                  <  6:218262/146‑1 (MQ=255)
                                              gggCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTg                                                                                                  >  5:218262/1‑146 (MQ=255)
                                              gggCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGAc                                                                                               <  5:184702/149‑1 (MQ=255)
                                                    ccACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGa                                                                                          <  1:67951/147‑1 (MQ=255)
                                                             gctgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATc                                                                                >  6:304001/1‑149 (MQ=255)
                                                              ctgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCg                                                                               <  5:309094/149‑1 (MQ=255)
                                                                gcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCgg                                                                              >  4:222085/1‑148 (MQ=255)
                                                                   accaGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTg                                                                          >  4:51026/1‑149 (MQ=255)
                                                                   accaGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTg                                                                          >  5:141848/1‑149 (MQ=255)
                                                                            tCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGaa                                                                 <  3:156961/149‑1 (MQ=255)
                                                                            tCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGaa                                                                 <  7:96932/149‑1 (MQ=255)
                                                                              cGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAg                                                                <  7:17433/148‑1 (MQ=255)
                                                                               ggCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGa                                                               >  1:205896/1‑148 (MQ=255)
                                                                                 ctgctTGGCAATGATCAGCTCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACAt                                                            <  4:245784/149‑1 (MQ=255)
                                                                                         cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAgg                                                    >  4:263889/1‑149 (MQ=255)
                                                                                          aCTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGTTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGc                                                   <  1:94444/149‑1 (MQ=255)
                                                                                          aCTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGc                                                   >  3:75426/1‑149 (MQ=255)
                                                                                               tCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGaa                                                                      <  2:98755/125‑1 (MQ=255)
                                                                                               tCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGaa                                                                      >  1:98755/1‑125 (MQ=255)
                                                                                                cAGCGCCAGCCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTc                                              <  4:265693/148‑1 (MQ=255)
                                                                                                 aGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCa                                             <  7:8554/148‑1 (MQ=255)
                                                                                                 aGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAg                                            <  3:222085/149‑1 (MQ=255)
                                                                                                       aGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTc                                      <  3:51026/149‑1 (MQ=255)
                                                                                                         gCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcgc                                    <  7:140796/149‑1 (MQ=255)
                                                                                                           cGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGcc                                                  <  6:136821/133‑1 (MQ=255)
                                                                                                           cGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGcc                                                  >  5:136821/1‑133 (MQ=255)
                                                                                                                      tGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATg                       <  7:46042/149‑1 (MQ=255)
                                                                                                                             gcTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTtctc                 >  1:70558/1‑148 (MQ=255)
                                                                                                                              cTTCACGGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGTTGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCga               <  3:263889/149‑1 (MQ=255)
                                                                                                                              cTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCga               <  2:70558/149‑1 (MQ=255)
                                                                                                                               ttCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgat              >  1:26792/1‑149 (MQ=255)
                                                                                                                                 cACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgatag            >  5:204087/1‑147 (MQ=255)
                                                                                                                                 cACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCga               <  6:204087/146‑1 (MQ=255)
                                                                                                                                   ccgacgCCCTTAAGGCGGCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGatgatg          >  5:193655/1‑149 (MQ=255)
                                                                                                                                           cTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATg                       >  4:175805/1‑129 (MQ=255)
                                                                                                                                            ttAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTtgtg  >  7:271418/1‑149 (MQ=255)
                                                                                                                                             tAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTtgtg  >  2:229577/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTTGTG  >  NC_002947/6171720‑6172007

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: