Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,861 A→T 5.3% 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 5.3% 227.8 / 2.8 75L147L (CTT→CTAglmSL‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base T (2/2);  total (35/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 = 3.62e-01

AGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTTGTG  >  NC_002947/6171715‑6172007
                                                                                                                                                   |                                                                                                                                                  
aGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAg                                                                                                                                                  <  7:226814/149‑1 (MQ=255)
 gCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAgg                                                                                                                                                 <  6:37988/149‑1 (MQ=255)
 gCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAgg                                                                                                                                                 <  7:125253/149‑1 (MQ=255)
        ggACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCg                                                                                                                                           <  1:299448/148‑1 (MQ=255)
              cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGt                                                                                                                                    >  3:253903/1‑149 (MQ=255)
                 gAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTg                                                                                                                                 <  4:57651/149‑1 (MQ=255)
                 gAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTg                                                                                                                                 <  1:87325/149‑1 (MQ=255)
                      tctcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATg                                                                                                                            >  5:137894/1‑149 (MQ=255)
                       ctcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGc                                                                                                                           >  4:251613/1‑149 (MQ=255)
                          gCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                                                                                                         <  1:160560/148‑1 (MQ=255)
                               ggCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAgg                                                                                                                   <  5:88906/149‑1 (MQ=255)
                               ggCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAgg                                                                                                                   <  4:207928/149‑1 (MQ=255)
                                  ccAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                                                                                                    <  5:289791/145‑1 (MQ=255)
                                  ccAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                                                                                                    >  6:289791/1‑145 (MQ=255)
                                  ccAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGgtg                                                                                                                >  7:27294/1‑149 (MQ=255)
                                   cAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCCCCTTAAGGGCGTCGGTCAGGACTGGAATGCTCTTCAGGgtgt                                                                                                               >  7:222962/1‑149 (MQ=255)
                                    aGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGtgtg                                                                                                              >  5:173064/1‑149 (MQ=255)
                                       ccGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggtg                                                                                                           <  8:27294/149‑1 (MQ=255)
                                         gATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGt                                                                                                                                             <  2:238062/113‑1 (MQ=255)
                                         gATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGt                                                                                                                                             >  1:238062/1‑113 (MQ=255)
                                           tAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGAt                                                                                                        >  5:274613/1‑148 (MQ=255)
                                                aGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGCGGCGGtgaa                                                                                                                                      >  5:207968/1‑110 (MQ=255)
                                                aGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGCGGCGGtgaa                                                                                                                                      <  6:207968/113‑4 (MQ=255)
                                                 gCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTgg                                                                                                 <  6:137894/149‑1 (MQ=255)
                                                 gCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTg                                                                                                  >  3:267581/1‑148 (MQ=255)
                                                   gggCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGAc                                                                                               >  4:47144/1‑149 (MQ=255)
                                                     gCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGa                                                                                             >  1:183482/1‑149 (MQ=255)
                                                      cTACCAGTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGAACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGATTCACCGCCGCCTTAAGTGCGTCGGTCAGGTCTGGAATTCTGTTCAGTGTGTGGTGGATCAGGTGGACGAt                                                                                            <  8:90941/149‑1 (MQ=255)
                                                      cTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGAt                                                                                            <  4:253903/149‑1 (MQ=255)
                                                      cTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGAt                                                                                            <  6:274613/149‑1 (MQ=255)
                                                      cTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGAt                                                                                            >  4:127569/1‑149 (MQ=255)
                                                          cACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACt                                                                                        >  6:4367/1‑149 (MQ=255)
                                                            cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGTATGCTCTTCAGGGTGTGGTGGATCAg                                                                                                     <  1:9243/134‑1 (MQ=255)
                                                            cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAg                                                                                                     >  2:9243/1‑134 (MQ=255)
                                                               cGTGCTGCCACCAGGCGGTCCGGCTGCTTGTCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                                                     <  6:173064/147‑1 (MQ=255)
                                                                gTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                                                     <  4:136588/146‑1 (MQ=255)
                                                                gTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                                                     >  3:136588/1‑146 (MQ=255)
                                                                   ctgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCg                                                                               <  5:187995/149‑1 (MQ=255)
                                                                    tgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCgg                                                                              >  7:159009/1‑149 (MQ=255)
                                                                        accaGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTg                                                                          >  5:100050/1‑149 (MQ=255)
                                                                                    ggCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGa                                                               >  5:176464/1‑148 (MQ=255)
                                                                                    ggCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGAc                                                              <  3:127569/149‑1 (MQ=255)
                                                                                       tgctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCt                                                                           <  7:69775/133‑1 (MQ=255)
                                                                                       tgctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCt                                                                           >  8:69775/1‑133 (MQ=255)
                                                                                         ctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGc                                                         >  5:140863/1‑149 (MQ=255)
                                                                                             gCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCa                                                      <  6:176464/148‑1 (MQ=255)
                                                                                              cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACt                                                                                        <  2:298657/113‑1 (MQ=255)
                                                                                              cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACt                                                                                        >  1:298657/1‑113 (MQ=255)
                                                                                              cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAgg                                                    <  2:183482/149‑1 (MQ=255)
                                                                                              cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAgg                                                    <  8:159009/149‑1 (MQ=255)
                                                                                               aCTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGc                                                   >  5:53640/1‑149 (MQ=255)
                                                                                                   aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCt                                                >  7:128746/1‑148 (MQ=255)
                                                                                                    tCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTc                                              <  3:251613/149‑1 (MQ=255)
                                                                                                     cAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTc                                              <  8:222962/148‑1 (MQ=255)
                                                                                                      aGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAg                                            <  6:53640/149‑1 (MQ=255)
                                                                                                       gcgcCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGc                                           >  5:237497/1‑149 (MQ=255)
                                                                                                       gcgcCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGc                                           >  3:302281/1‑149 (MQ=255)
                                                                                                             ggCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcg                                     <  6:237497/149‑1 (MQ=255)
                                                                                                              gCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcgc                                    >  1:16774/1‑149 (MQ=255)
                                                                                                               ccGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcgcg                                   <  2:16774/149‑1 (MQ=255)
                                                                                                                cGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTGAGGGTGTGGTGGATCAGGTGGCCGATGACTTCGGTATCGGTCTGCGATGTGAAGACATAGCCGAGGCCCGGCAGCTCTTCgggc                                  >  4:225815/1‑149 (MQ=255)
                                                                                                                 tgaaGCGCCTGCACGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcgcg                                   >  7:34252/3‑147 (MQ=255)
                                                                                                                 gTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGt                                                                    <  6:213480/114‑1 (MQ=255)
                                                                                                                 gTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGt                                                                    >  5:213480/1‑114 (MQ=255)
                                                                                                                  tAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAg                                <  4:267581/149‑1 (MQ=255)
                                                                                                                         cGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCa                         <  1:234433/149‑1 (MQ=255)
                                                                                                                                 gcgcTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATg                       >  6:182982/1‑144 (MQ=255)
                                                                                                                                 gcgcTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATg                       <  5:182982/144‑1 (MQ=255)
                                                                                                                                  cgcTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCTGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTtctc                 <  5:221339/149‑1 (MQ=255)
                                                                                                                                     ttCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgat              >  5:124159/1‑149 (MQ=255)
                                                                                                                                       cACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgatg             <  6:140863/148‑1 (MQ=255)
                                                                                                                                       cACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgatg             <  3:47144/148‑1 (MQ=255)
                                                                                                                                        aCCGACGCCCTTAAGGCGGCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTAt                                                                                 >  7:283228/1‑79 (MQ=38)
                                                                                                                                        aCCGACGCCCTTAAGGCGGCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTAt                                                                                 <  8:283228/79‑1 (MQ=38)
                                                                                                                                                 ttAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTtgtg  >  3:224384/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGCCGTTGTG  >  NC_002947/6171715‑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: