Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,893 Δ1 bp 9.2% coding (465/1836 nt) 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,8930G. 9.2% 202.9 / 15.1 65coding (465/1836 nt)glmSL‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base . (3/3);  total (32/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 = 8.19e-01

TCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATG  >  NC_002947/6171750‑6172042
                                                                                                                                               |                                                                                                                                                      
tCAGTAACCTGACGCAGCGCCAGCTGGTCGGCCGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGCGCTACCCCTGCGTGTTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                                 <  2:303983/150‑1 (MQ=255)
tCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                                 >  2:115015/1‑150 (MQ=255)
tCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                                 >  1:771946/1‑150 (MQ=255)
    tAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTACTGCAACCAGGCGATCCAGCTACTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCt                                                                                                                                             >  1:520267/1‑150 (MQ=255)
     aaCCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCAATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCtt                                                                                                                                            >  1:792043/1‑150 (MQ=255)
      accgaaCGCAGCGCCAGCTGGTCCGACGCCAGGAAGTCCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGGGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTc                                                                                                                                           <  2:413530/145‑1 (MQ=255)
        cTGACGCAGCGCCAGCTGTACCGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCTGCACGGCGCttacg                                                                                                                                         <  1:56431/150‑4 (MQ=255)
        cTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAc                                                                                                                                         >  2:775502/1‑150 (MQ=255)
         tGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAcc                                                                                                                                        >  2:763225/1‑150 (MQ=255)
           aCGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCAATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAccgc                                                                                                                                      <  2:792043/150‑1 (MQ=255)
               aGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCgccgcc                                                                                                                                  <  2:771946/150‑1 (MQ=255)
                   ccAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTaa                                                                                                                              <  1:1116248/150‑1 (MQ=255)
                    cAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATCACCAGCGGACTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAg                                                                                                                             <  1:1098438/150‑1 (MQ=255)
                           tCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCg                                                                                                                      >  1:962476/1‑150 (MQ=255)
                                cGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAg                                                                                                                 <  1:321710/150‑1 (MQ=255)
                                   cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAg                                                                                                                                                  >  1:399992/1‑114 (MQ=255)
                                   cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAg                                                                                                                                                  <  2:399992/114‑1 (MQ=255)
                                   cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTc                                                                                                              >  2:884857/1‑150 (MQ=255)
                                        aGTTCTCGCCAAGGCCCAGGACGATAACCAGCGGGTTATCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGaa                                                                                                         <  2:336100/150‑1 (MQ=255)
                                          ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATg                                                                                                       >  2:86507/1‑150 (MQ=255)
                                          ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATg                                                                                                       <  1:971519/150‑1 (MQ=255)
                                             tcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                                                                                    >  1:454628/1‑150 (MQ=255)
                                                 cAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCa                                                                                                >  1:57938/1‑150 (MQ=255)
                                                 cAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCa                                                                                                >  2:771058/1‑150 (MQ=255)
                                                  aaGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                                                                               >  2:474453/1‑150 (MQ=255)
                                                       ccAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGgtgt                                                                                          <  2:961874/150‑1 (MQ=255)
                                                          ggCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                       <  2:962476/150‑1 (MQ=255)
                                                           gCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggtg                                                                                      <  1:267132/150‑1 (MQ=255)
                                                            ccGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGtggtgg                                                                                     <  1:61034/150‑1 (MQ=255)
                                                                  aCCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                               <  1:316027/150‑1 (MQ=255)
                                                                   ccAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGt                                                                              >  2:183447/1‑150 (MQ=255)
                                                                             aCCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGAc                                                                    <  1:771058/150‑1 (MQ=255)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGTCGGAAGCGCCTGCACGGCGCttac                                                                                                                                          <  1:813114/76‑3 (MQ=25)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGTCGGAAGCGCCTGCACGGCGCttac                                                                                                                                          >  2:813114/1‑74 (MQ=25)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                               >  2:1004058/1‑135 (MQ=255)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                               <  1:1004058/135‑1 (MQ=255)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                                >  1:690874/1‑150 (MQ=255)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCTGCACGGCGCttacg                                                                                                                                         >  2:802210/1‑74 (MQ=38)
                                                                                 cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCTGCACGGCGCttacg                                                                                                                                         <  1:802210/77‑4 (MQ=38)
                                                                                     gTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTAt                                                            <  2:26447/150‑1 (MQ=255)
                                                                                       gctgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCg                                                          <  1:115015/150‑1 (MQ=255)
                                                                                         tgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGCGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGt                                                        <  1:1061959/150‑1 (MQ=255)
                                                                                           ccaccaGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCt                                                      <  1:884857/150‑1 (MQ=255)
                                                                                                     gTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATc                                                           >  2:1029310/1‑135 (MQ=255)
                                                                                                     gTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATc                                                           <  1:1029310/135‑1 (MQ=255)
                                                                                                          gctgctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACAt                                       >  1:754176/1‑150 (MQ=255)
                                                                                                                tGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCa                                 <  2:57938/150‑1 (MQ=255)
                                                                                                                  gCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                       <  2:906225/94‑1 (MQ=255)
                                                                                                                  gCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                       >  1:906225/1‑94 (MQ=255)
                                                                                                                   cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGc                              >  2:751564/1‑150 (MQ=255)
                                                                                                                      tGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCt                           >  2:1124948/1‑150 (MQ=255)
                                                                                                                       gATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCtt                          <  2:435418/150‑1 (MQ=255)
                                                                                                                       gATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCtt                          >  2:463322/1‑150 (MQ=255)
                                                                                                                        aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTa                                                             >  2:738272/1‑114 (MQ=255)
                                                                                                                        aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTa                                                             <  1:738272/114‑1 (MQ=255)
                                                                                                                           aGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGc                      >  1:709738/1‑150 (MQ=255)
                                                                                                                           aGCGCCAGGCCGTAAGAGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGc                      >  2:739354/1‑150 (MQ=255)
                                                                                                                            gcgcCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGct                     >  1:74234/1‑150 (MQ=255)
                                                                                                                                cAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGa                                                                     >  1:1018917/1‑98 (MQ=255)
                                                                                                                                cAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGa                                                                     <  2:1018917/98‑1 (MQ=255)
                                                                                                                                 aGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcg                <  1:474453/150‑1 (MQ=255)
                                                                                                                                   gCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                       <  1:1061833/77‑1 (MQ=255)
                                                                                                                                   gCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                       >  2:1061833/1‑77 (MQ=255)
                                                                                                                                    ccGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCgcgc             <  2:520267/150‑1 (MQ=255)
                                                                                                                                     gtgaaGCGCCTGCACGGCGCTTCACCGCCCCCTTACGGGCGTCGGTCAGGTCTGGACTGCTCTTCCGAGTGTGGTGGATCAGGTGGCCGATGCCTTCGGTATCGGTCTGCGACGTGAACACATCGCCCCGGCCCTTCAGCTCTTCgcgcc            >  2:424310/4‑149 (MQ=255)
                                                                                                                                               gTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATg  <  1:739354/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
TCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATG  >  NC_002947/6171750‑6172042

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