Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 413,107 C→G 5.6% 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 5.6% 224.3 / 4.9 73A686A (GCC→GCGglnEglutamate‑ammonia‑ligase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base G (2/2);  total (37/36)
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.60e-01

CGCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC  >  NC_002947/412958‑413256
                                                                                                                                                     |                                                                                                                                                     
cgcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGcc                                                                                                                                                       <  2:1733470/150‑1 (MQ=255)
 gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCNGAAATCAGAGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      >  1:1227755/1‑150 (MQ=255)
 gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  2:397334/150‑1 (MQ=255)
 gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  2:801665/150‑1 (MQ=255)
  ccATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCg                                                                                                                                                     >  1:703680/1‑150 (MQ=255)
             tGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAc                                                                                                                                          <  2:723990/150‑1 (MQ=255)
                 ccACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCaa                                                                                                                                      <  1:360207/150‑1 (MQ=255)
                 ccACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGAATACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCaa                                                                                                                                      >  2:468673/1‑150 (MQ=255)
                 ccACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAACACCGCGACTACCTGACCGGTCTTGCTTAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCaa                                                                                                                                      >  1:725430/1‑150 (MQ=255)
                  cacaGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAg                                                                                                                                     >  1:1253856/1‑150 (MQ=255)
                   acaGCCTGCGCGTGGCCGCCACGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTTACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAgc                                                                                                                                    <  1:336066/150‑1 (MQ=255)
                           cgcgTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGAc                                                                                                                            >  1:481056/1‑150 (MQ=255)
                               tGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCa                                                                                                                        <  2:703680/150‑1 (MQ=255)
                                   cgccTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCTACATAAGGGGAGCGACGGCAGCCt                                                                                                                    <  2:1142498/150‑1 (MQ=255)
                                   cgccTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCt                                                                                                                    <  2:857323/150‑1 (MQ=255)
                                       tCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGc                                                                                                                >  1:201201/1‑150 (MQ=255)
                                          gAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGAc                                                                                                             <  2:201201/150‑1 (MQ=255)
                                              tCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCgg                                                                                                         <  2:1459613/150‑1 (MQ=255)
                                               cAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAcc                                                                                                                                         >  1:1326315/1‑117 (MQ=255)
                                               cAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAcc                                                                                                                                         <  2:1326315/117‑1 (MQ=255)
                                                aGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGt                                                                                                       <  1:5767/150‑1 (MQ=255)
                                                 tcGGCAACCTGCCTTTGATGAAAGTCAGCGACTACCTGACCTAGCTGGCTGAAGCCATCCTCGACTAGGTGCTGCCCCTGGCCGGGCGCCAGACCGTGGCCCGCCACGGCCAAGCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGtt                                                                                                      <  2:336723/149‑1 (MQ=255)
                                                 gCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGtt                                                                                                      <  2:1768579/150‑1 (MQ=255)
                                                   ggCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAACCATCCTCGACCAAGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGGTtc                                                                                                    >  1:159329/1‑150 (MQ=255)
                                                    gCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtca                                                                                                   >  1:1792525/1‑150 (MQ=255)
                                                     cAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcat                                                                                                  <  2:1227755/150‑1 (MQ=255)
                                                     aaaCCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACTCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcat                                                                                                  <  1:555378/149‑1 (MQ=255)
                                                         cTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGAGCGAGGACGTAGCCGGCCTCGGCCAGCCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCAtt                                                                                              >  1:228786/1‑150 (MQ=255)
                                                            ccGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCaa                                                                                                                                           >  2:1469154/1‑102 (MQ=255)
                                                            ccGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCaa                                                                                                                                           <  1:1469154/102‑1 (MQ=255)
                                                             cGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                          <  2:481056/150‑1 (MQ=255)
                                                                  atgaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtct                                                                                                                                            <  1:890418/95‑4 (MQ=255)
                                                                  atgaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtct                                                                                                                                            >  2:890418/1‑92 (MQ=255)
                                                                        gTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  1:911132/79‑1 (MQ=255)
                                                                        gTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      >  2:911132/1‑79 (MQ=255)
                                                                         tCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCGGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGa                                                                              >  1:1772840/1‑150 (MQ=255)
                                                                             cGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAc                                                                                                                                          <  2:1438206/86‑1 (MQ=255)
                                                                             cGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAc                                                                                                                                          >  1:1438206/1‑86 (MQ=255)
                                                                                 tACCTGACCTGGCTGGCTGAAGCCATCCTCGTCCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggcggc                                                                      >  2:1916471/1‑150 (MQ=255)
                                                                                 gacctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggcggc                                                                      >  2:487074/2‑150 (MQ=255)
                                                                                  acctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGcc                                                                     <  2:1253856/150‑1 (MQ=255)
                                                                                    ctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcat                                                                                                  >  2:378208/1‑119 (MQ=255)
                                                                                    ctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcat                                                                                                  <  1:378208/119‑1 (MQ=255)
                                                                                       acctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAg                                                                <  1:751878/150‑1 (MQ=255)
                                                                                           ggctggctGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                            <  2:1792525/150‑1 (MQ=255)
                                                                                           ggctggctGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGCCTACGGCAAGATGGGCGGCCTGGAGCTCg                                                            <  1:544142/150‑1 (MQ=255)
                                                                                              tggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc                                                                                                 <  1:391617/110‑1 (MQ=255)
                                                                                              tggctgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc                                                                                                 >  2:391617/1‑110 (MQ=255)
                                                                                                  tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACgg                                                     >  1:1324051/1‑150 (MQ=255)
                                                                                                       ccATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGc                                                                                                                <  2:1128074/86‑1 (MQ=255)
                                                                                                       ccATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGc                                                                                                                >  1:1128074/1‑86 (MQ=255)
                                                                                                       ccATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgg                                                >  1:1925633/1‑150 (MQ=255)
                                                                                                       ccATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgg                                                <  2:1838135/150‑1 (MQ=255)
                                                                                                        cATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgga                                               >  2:553189/1‑150 (MQ=255)
                                                                                                          tCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacc                                             >  1:218108/1‑150 (MQ=255)
                                                                                                            cTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctg                                           >  2:1649911/1‑150 (MQ=255)
                                                                                                             tCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctgg                                          >  1:988758/1‑150 (MQ=255)
                                                                                                                  cAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTACATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggaactg                                     >  1:27490/1‑150 (MQ=255)
                                                                                                                     gTGCTGGCCGTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTc                                  <  2:1428097/150‑1 (MQ=255)
                                                                                                                     gTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc                                                                                                 >  2:1713450/1‑87 (MQ=255)
                                                                                                                     gTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc                                                                                                 <  1:1713450/87‑1 (MQ=255)
                                                                                                                     gTGCAGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTAAGGCAAGATGGGCGGCCTGGAGCGCGGCCACGGTTCGGACCTGGACCTGGTc                                  >  2:513421/1‑150 (MQ=255)
                                                                                                                         ctgccctggccTGGCGCCAGACCGTGGCCCCCGCCGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGCGCGGCCTGGAGCTCGGCCACGGTTCGAACCTGGACCTGGTCTTCa                              <  2:1891941/148‑1 (MQ=255)
                                                                                                                          ggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctgg                                          <  2:30965/137‑1 (MQ=255)
                                                                                                                          ggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctgg                                          >  1:30965/1‑137 (MQ=255)
                                                                                                                           gccctggccTGGCGCCAGACCGTGGATATCCACGGCCAACCCAAGCGCAGCGACGGCAGCATGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATc                            >  1:1321666/1‑150 (MQ=255)
                                                                                                                            ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTTCGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCAtag                           >  1:1891941/1‑148 (MQ=255)
                                                                                                                            ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATcc                           <  2:725430/150‑1 (MQ=255)
                                                                                                                                    tggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCGCGGCCACGGTTCGGACCTGGACCTGGCCTTCATCCATGACGGc                   >  2:764518/1‑150 (MQ=255)
                                                                                                                                        gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg                                     <  2:547036/128‑1 (MQ=255)
                                                                                                                                        gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg                                     >  1:547036/1‑128 (MQ=255)
                                                                                                                                          ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgga                                               >  2:665940/5‑116 (MQ=255)
                                                                                                                                          ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgga                                               <  1:665940/112‑1 (MQ=255)
                                                                                                                                           aGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGc            >  2:1904723/1‑150 (MQ=255)
                                                                                                                                            gACCGTGGCCCGCCACGGCCAACCCAAGCGCCGCGCGATCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCa           >  1:1125381/1‑150 (MQ=255)
                                                                                                                                             aCCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAg          <  1:747613/150‑1 (MQ=255)
                                                                                                                                              ccGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAgg         <  1:468673/150‑1 (MQ=255)
                                                                                                                                                gTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctggac                                        >  1:377701/1‑117 (MQ=255)
                                                                                                                                                gTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctggac                                        <  2:377701/117‑1 (MQ=255)
                                                                                                                                                     ccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                            <  1:378988/92‑1 (MQ=255)
                                                                                                                                                     ccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                            >  2:378988/1‑92 (MQ=255)
                                                                                                                                                     ccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAc  >  2:657373/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC  >  NC_002947/412958‑413256

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