Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 413,107 C→G 8.9% A686A (GCC→GCG glnE → glutamate‑ammonia‑ligase adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947413,1070CG 8.9% 231.4 / 11.0 76A686A (GCC→GCGglnEglutamate‑ammonia‑ligase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base G (2/5);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGG  >  NC_002947/412958‑413252
                                                                                                                                                      |                                                                                                                                                 
cgcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGcc                                                                                                                                                   >  1:1512278/1‑150 (MQ=255)
  ccATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGCCCGTGGCCCg                                                                                                                                                 >  1:8874/1‑150 (MQ=255)
       tCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACg                                                                                                                                            >  2:129358/1‑150 (MQ=255)
          aGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtc                                                                                                                                         <  1:129358/150‑3 (MQ=255)
             ttgCCCACATCCTGCGCGTCGCAGCCTAGCAAGACTGCGGCAATCCGGCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtctgg                                                                                                                                      <  2:429691/148‑6 (MQ=255)
             tGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGGCCGTGGCGGGCCACggtctgg                                                                                                                                      <  2:478557/150‑6 (MQ=255)
               gCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAccc                                                                                                                                    >  1:389685/1‑150 (MQ=255)
                 ccACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCaa                                                                                                                                  >  1:468871/1‑150 (MQ=255)
                  cacaGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAg                                                                                                                                 >  1:275744/1‑150 (MQ=255)
                  cacaGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAg                                                                                                                                 <  2:436325/150‑1 (MQ=255)
                                 gccgccTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGc                                                                                                                  <  2:1512278/150‑1 (MQ=255)
                                  ccgccTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGcc                                                                                                                 >  2:1054585/1‑150 (MQ=255)
                                  ccgccTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGcc                                                                                                                 >  1:699762/1‑150 (MQ=255)
                                              tCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCgg                                                                                                     <  1:1054585/150‑1 (MQ=255)
                                              tCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGATACCAAGGCGGGGGCGCCAGACCGTGGCCCGCCACGGCCAACCACAGCCACCCGACGGCAGCCTGTGCGTCCCgg                                                                                                     >  1:1593911/1‑150 (MQ=255)
                                                 gCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGtt                                                                                                  <  2:64242/150‑1 (MQ=255)
                                                   ggCAACCTGCCGTTGATGAAAGTCGGCCACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGTGTGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtc                                                                                                <  1:359927/150‑1 (MQ=255)
                                                      aaCCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc                                                                                             <  2:960989/150‑1 (MQ=255)
                                                      aaCCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc                                                                                             >  1:9855/1‑150 (MQ=255)
                                                       aCCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatca                                                                                            <  2:1030382/150‑1 (MQ=255)
                                                            ccGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATc                                                                                       <  2:1355864/150‑1 (MQ=255)
                                                             cGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                      <  2:1388265/150‑1 (MQ=255)
                                                             cGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                      <  2:699762/150‑1 (MQ=255)
                                                                    gaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACgg                                                                               >  2:1363379/1‑150 (MQ=255)
                                                                      aaGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCa                                                                             <  1:1056902/150‑1 (MQ=255)
                                                                              gACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggc                                                                     <  2:275744/150‑1 (MQ=255)
                                                                                 tACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggcggc                                                                  >  1:716331/1‑150 (MQ=255)
                                                                                 tACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggcggc                                                                  >  2:330347/1‑150 (MQ=255)
                                                                                   cctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGt                                                                                                   <  1:252398/115‑1 (MQ=255)
                                                                                   cctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGt                                                                                                   >  2:252398/1‑115 (MQ=255)
                                                                                    ctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTg                                                               >  1:103744/1‑150 (MQ=255)
                                                                                      gacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGa                                                             <  2:963001/150‑1 (MQ=255)
                                                                                          tggctggcTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGAt                                                                         <  2:1267689/134‑1 (MQ=255)
                                                                                          tggctggcTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGAt                                                                         >  1:1267689/1‑134 (MQ=255)
                                                                                          gggctggctGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGGCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCTTTATCGGTTACGGCTAGATGGGCGGCCTGGAGCTc                                                         <  2:669485/149‑1 (MQ=255)
                                                                                            gctggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTAc                                                                                 <  2:637214/124‑1 (MQ=255)
                                                                                            gctggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTAc                                                                                 >  1:637214/1‑124 (MQ=255)
                                                                                              tggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGcc                                                     >  2:960620/1‑150 (MQ=255)
                                                                                                gctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCa                                                                                                                    >  2:1121811/1‑85 (MQ=255)
                                                                                                gctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCa                                                                                                                    <  1:1121811/85‑1 (MQ=255)
                                                                                                 ctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACg                                                  <  1:1285596/150‑1 (MQ=255)
                                                                                                  tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                        >  1:855146/1‑143 (MQ=255)
                                                                                                  tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                        <  2:855146/143‑1 (MQ=255)
                                                                                                    aaGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGtt                                               <  1:1303131/150‑1 (MQ=255)
                                                                                                      gCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggc                                                                     <  2:621543/126‑1 (MQ=255)
                                                                                                      gCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGGTGCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCg                                             >  2:110836/1‑150 (MQ=255)
                                                                                                      gCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACACGGGTTTCATCATTATCGGTTACGGCAAGATGggc                                                                     >  1:621543/1‑126 (MQ=255)
                                                                                                          tCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacc                                         >  1:3797/1‑150 (MQ=255)
                                                                                                           ccTGGACCGGGTGATGTGCCGGGCCTGGCGCGAGACGGGGGCGCGGCGGGGGCAACGGAAGCGCAGGGACGGCAGCCTGTGCGGCTCGGGTTTCATCATGATCGGGTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacct                                        <  2:666596/150‑1 (MQ=255)
                                                                                                           ccTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacct                                        <  1:960620/150‑1 (MQ=255)
                                                                                                              cGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctgga                                     <  2:9855/150‑1 (MQ=255)
                                                                                                                aCCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCt                                                                                                                >  1:1212587/1‑73 (MQ=255)
                                                                                                                aCCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCt                                                                                                                <  2:1212587/73‑1 (MQ=255)
                                                                                                                   aGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGacctggacctgg                                <  2:468871/150‑1 (MQ=255)
                                                                                                                   aGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGCGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGacctggacctgg                                <  2:1182756/150‑1 (MQ=255)
                                                                                                                    ggTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTa                                                                                  <  2:190116/99‑1 (MQ=255)
                                                                                                                    ggTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTa                                                                                  >  1:190116/1‑99 (MQ=255)
                                                                                                                      tGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCt                             <  2:103744/150‑1 (MQ=255)
                                                                                                                        ctggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggac                                          >  1:506797/1‑135 (MQ=255)
                                                                                                                        ctggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggac                                          <  2:506797/135‑1 (MQ=255)
                                                                                                                        ctggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTc                           <  2:8874/150‑1 (MQ=255)
                                                                                                                        ctggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTc                           <  1:330347/150‑1 (MQ=255)
                                                                                                                            ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctggacc                                   <  2:1469502/138‑1 (MQ=255)
                                                                                                                            ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctggacc                                   >  1:1469502/1‑138 (MQ=255)
                                                                                                                              ctggcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCCCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCAt                     <  1:1363379/150‑1 (MQ=255)
                                                                                                                                 gcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggac                                          >  2:1374602/1‑126 (MQ=255)
                                                                                                                                 gcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggac                                          <  1:1374602/126‑1 (MQ=255)
                                                                                                                                    ctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                      <  2:822864/80‑1 (MQ=255)
                                                                                                                                    ctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                      >  1:822864/1‑80 (MQ=255)
                                                                                                                                      ggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCg              >  2:179813/1‑150 (MQ=255)
                                                                                                                                         gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATTACGGCGATc           >  1:1483667/1‑150 (MQ=255)
                                                                                                                                            aGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATc                        <  2:610388/134‑1 (MQ=255)
                                                                                                                                            aGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCAGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATc                        >  1:610388/1‑134 (MQ=255)
                                                                                                                                              gccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTa                                                                                         >  1:163647/2‑67 (MQ=38)
                                                                                                                                              gccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTa                                                                                         <  2:163647/66‑1 (MQ=38)
                                                                                                                                                cGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGc    <  1:179813/150‑1 (MQ=255)
                                                                                                                                                  tGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGAACTGGTCTTCATCCATGACGGCGATCCGCAGGCgg  >  1:421789/1‑150 (MQ=255)
                                                                                                                                                     cccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGc    <  1:319579/145‑1 (MQ=255)
                                                                                                                                                     cccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGc    >  2:319579/1‑145 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CGCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGG  >  NC_002947/412958‑413252

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