Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,080,964 A→T 8.0% intergenic (‑144/+19) amaC ← / ← dpkA D‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,080,9640AT 8.0% 209.5 / 9.1 75intergenic (‑144/+19)amaC/dpkAD‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base T (3/3);  total (42/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 = 5.32e-01

AACACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCT  >  NC_002947/4080816‑4081113
                                                                                                                                                     |                                                                                                                                                     
aaCACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTaaa                                                                                                                                                      >  2:803488/1‑150 (MQ=255)
   acaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGTATGAGGGCTGGCTTTTTCGCGGGTAAAccc                                                                                                                                                   >  2:1018397/1‑150 (MQ=255)
   acaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAccc                                                                                                                                                   >  1:1640823/1‑150 (MQ=255)
     aGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACGGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcgc                                                                                                                                                 >  2:1282253/1‑150 (MQ=255)
        gTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGaa                                                                                                                                              <  1:1758766/150‑1 (MQ=255)
          gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  1:1018397/150‑1 (MQ=255)
          gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  1:898960/150‑1 (MQ=255)
           cAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTTTACCCGCGaaaaa                                                                                                                                           <  2:1083118/150‑4 (MQ=255)
           cAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:803488/150‑1 (MQ=255)
             gCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                         <  1:180406/150‑1 (MQ=255)
             gCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                         <  2:1244815/150‑1 (MQ=255)
             gCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGACTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                         <  1:1445508/150‑1 (MQ=255)
                 cTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCaaa                                                                                                                                     >  1:503001/1‑150 (MQ=255)
                    ggATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCa                                                                                                                                  <  2:1640823/150‑1 (MQ=255)
                      aTCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGTCTTTTTCGCGGGTAAACCCGCGACGGGGCCAAATCAGc                                                                                                                                >  1:511006/1‑150 (MQ=255)
                           gACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTAGTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagc                                                                                                                           >  2:642580/1‑150 (MQ=255)
                             cGGTACGCACAGTCCGCGTTAGACAAACCCCGATTGCGGGTTCAGGAACATCCACTATTGTGCACGGTTACGGATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagcag                                                                                                                         <  1:417528/150‑1 (MQ=255)
                                   ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATACACTATTGTGCAAGGTTGTCTATGCATCATTGCATGACCGGGATGAGGGCTTGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCttt                                                                                                                   >  1:977034/1‑150 (MQ=255)
                                           aGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACGGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCtt                                                                                                           <  1:1282253/150‑1 (MQ=255)
                                             tgtTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                         <  2:503001/150‑1 (MQ=255)
                                             tgtTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                         <  1:276758/150‑1 (MQ=255)
                                             ggtTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                         >  2:777385/2‑150 (MQ=255)
                                               ttGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCaa                                                                                                       <  1:777385/150‑1 (MQ=255)
                                                tGAGCATACCCGGAGTGCGGTTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                                      <  1:532162/150‑1 (MQ=255)
                                                 gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGAAGACCGTCGAGCAGTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACt                                                                                                     >  2:927748/1‑150 (MQ=255)
                                                  aGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                    >  1:1338299/1‑150 (MQ=255)
                                                  aGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                    >  1:1667026/1‑150 (MQ=255)
                                                   gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTcc                                                                                                   <  2:443517/150‑1 (MQ=255)
                                                    cATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCt                                                                                                  <  2:113160/150‑1 (MQ=255)
                                                         ccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAg                                                                                                                                 >  1:1823607/1‑114 (MQ=255)
                                                         ccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAg                                                                                                                                 <  2:1823607/114‑1 (MQ=255)
                                                         ccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                             <  2:1667026/150‑1 (MQ=255)
                                                          gggAGTTCGGTTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCgg                                                                                            <  2:679868/149‑1 (MQ=255)
                                                           ggAGTGCGGGTTCAGGAGCATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:1534695/102‑1 (MQ=255)
                                                           ggAGTGCGGGTTCAGGAGCATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           >  2:1534695/1‑102 (MQ=255)
                                                                 cGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTc                                                                                                                  <  1:1004002/121‑1 (MQ=255)
                                                                 cGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGAGATGAGGGCTGGCTTTTTCGCGGGGAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTc                                                                                                                  >  2:1004002/1‑121 (MQ=255)
                                                                   ggTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGcta                                                                                   >  1:1212929/1‑148 (MQ=255)
                                                                     ttCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAcc                                                                                 >  1:1174488/1‑150 (MQ=255)
                                                                      tCAGGAACATACACTATTGTGCGAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAccc                                                                                <  1:76508/150‑1 (MQ=255)
                                                                          gAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCCGCTCTTTCTGGCCTTTCAACTCCTGCTCGGTCACCGCCatacgttt                                                                            >  2:312078/1‑143 (MQ=255)
                                                                          gAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctc                                                                            >  1:1325449/1‑150 (MQ=255)
                                                                          gAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGATTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCtt                                                                                                           >  2:1710396/1‑119 (MQ=255)
                                                                          gAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGATTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCtt                                                                                                           <  1:1710396/119‑1 (MQ=255)
                                                                               tCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctcct                                                                       >  1:481572/1‑150 (MQ=255)
                                                                                    tATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                            >  2:1134773/1‑92 (MQ=255)
                                                                                    tATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                            <  1:1134773/92‑1 (MQ=255)
                                                                                      ttGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCATTTTCGCGGGTAAACCCGCGAAGGGGCCAATTCAGCCAAGCAGCTCTTTCAGTCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCAc                                                                <  2:979431/150‑1 (MQ=255)
                                                                                            aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCt                                                                                                            >  1:1820293/1‑100 (MQ=255)
                                                                                            aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCt                                                                                                            <  2:1820293/100‑1 (MQ=255)
                                                                                            aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcc                                                                           >  2:1529709/1‑133 (MQ=255)
                                                                                            aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcc                                                                           <  1:1529709/133‑1 (MQ=255)
                                                                                            aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgcg                                                          >  1:516976/1‑150 (MQ=255)
                                                                                             aGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGTTGGATTTCTTCGCGGGCAACCGTCGAATGTGCCTAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCgcgc                                                         >  1:514583/1‑150 (MQ=255)
                                                                                                       tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGTAAATGTTTTCGTTTCAACTCCTGCTCGGTCACGGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgc                                               >  1:1772811/1‑150 (MQ=255)
                                                                                                       tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgc                                               >  1:1266558/1‑150 (MQ=255)
                                                                                                       tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgc                                               >  2:1340273/1‑150 (MQ=255)
                                                                                                       tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgc                                               >  1:561685/1‑150 (MQ=255)
                                                                                                          aaCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGGAAACCCGGGAAGGGGCCAAGTCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCCCCTCCTCGGCCACCTCGGGCTCACGCTAGcgcc                                            >  1:1131758/1‑150 (MQ=255)
                                                                                                                 cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                         >  2:776363/1‑82 (MQ=255)
                                                                                                                 cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                         <  1:776363/82‑1 (MQ=255)
                                                                                                                  aTGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGccc                                    >  1:290120/1‑150 (MQ=255)
                                                                                                                       cGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCAt                               >  2:1615429/1‑150 (MQ=255)
                                                                                                                          gATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACAGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcg                            >  1:604844/1‑150 (MQ=255)
                                                                                                                           aTGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                           <  1:642580/150‑1 (MQ=255)
                                                                                                                            tGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCAGCCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGcgcg                          >  1:1260403/1‑150 (MQ=255)
                                                                                                                                ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAg                      >  2:1027131/1‑150 (MQ=255)
                                                                                                                                  ctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGcc                    >  1:1896469/1‑150 (MQ=255)
                                                                                                                                       tttttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAcc               >  1:1682539/1‑150 (MQ=255)
                                                                                                                                        cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAc                >  2:1436064/4‑149 (MQ=255)
                                                                                                                                        cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAc                <  1:1436064/146‑1 (MQ=255)
                                                                                                                                         ttttGCGGGTAACCCCGCAACGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACACTTGTTTAATCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc             <  1:594013/150‑1 (MQ=255)
                                                                                                                                         ccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCa                                                      >  2:303381/3‑110 (MQ=255)
                                                                                                                                         ccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCa                                                      <  1:303381/108‑1 (MQ=255)
                                                                                                                                          cttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc             >  1:791888/2‑150 (MQ=255)
                                                                                                                                              gcgGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTGCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGc         >  2:66412/1‑150 (MQ=255)
                                                                                                                                                     aaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCt  <  1:927452/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCT  >  NC_002947/4080816‑4081113

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: