Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 413,108 C→G 6.9% R687G (CGC→GGC)  glnE → glutamate‑ammonia‑ligase adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947413,1080CG 6.9% 224.2 / 6.8 73R687G (CGC→GGC) glnEglutamate‑ammonia‑ligase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base G (3/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 = 9.28e-01

GCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC  >  NC_002947/412959‑413256
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                            >  2:890418/1‑92 (MQ=255)
                                                                 atgaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtct                                                                                                                                            <  1:890418/95‑4 (MQ=255)
                                                                       gTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      >  2:911132/1‑79 (MQ=255)
                                                                       gTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGccc                                                                                                                                                      <  1:911132/79‑1 (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                                                <  2:1838135/150‑1 (MQ=255)
                                                                                                      ccATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgg                                                >  1:1925633/1‑150 (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                                          >  1:30965/1‑137 (MQ=255)
                                                                                                                         ggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctgg                                          <  2:30965/137‑1 (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                                     >  1:547036/1‑128 (MQ=255)
                                                                                                                                       gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg                                     <  2:547036/128‑1 (MQ=255)
                                                                                                                                         ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgga                                               <  1:665940/112‑1 (MQ=255)
                                                                                                                                         ttggccGTGGCGGGCCACGGTCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCgga                                               >  2:665940/5‑116 (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)
                                                                                                                                                     |                                                                                                                                                    
GCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC  >  NC_002947/412959‑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: