Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,080,964 A→T 6.6% 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 6.6% 265.8 / 4.3 87intergenic (‑144/+19)amaC/dpkAD‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base T (4/2);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCT  >  NC_002947/4080816‑4081113
                                                                                                                                                    |                                                                                                                                                     
aaCACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTaaa                                                                                                                                                      >  2:44221/1‑150 (MQ=255)
 acacAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAc                                                                                                                                                     >  1:852325/1‑150 (MQ=255)
  cacaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAcc                                                                                                                                                    <  1:94881/150‑1 (MQ=255)
  cacaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAcc                                                                                                                                                    <  1:1212088/150‑1 (MQ=255)
   acaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAccc                                                                                                                                                   >  1:831554/1‑150 (MQ=255)
     aGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcgc                                                                                                                                                 <  2:434733/150‑1 (MQ=255)
        gTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGaa                                                                                                                                              <  1:1879107/150‑1 (MQ=255)
        gTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGaa                                                                                                                                              <  2:1066581/150‑1 (MQ=255)
          gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  2:1473784/150‑1 (MQ=255)
          gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  2:1907974/150‑1 (MQ=255)
          gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            >  1:1222136/1‑150 (MQ=255)
          gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  2:1135086/150‑1 (MQ=255)
           cAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:2329257/150‑1 (MQ=255)
           cAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  2:726540/150‑1 (MQ=255)
             gCTCTTTGAATCAGGTTGGTACTGAAATTCAGTGTTGAGCATACCGGGAGTGCGCTTTGAGGAACATCTACTACTGTGCAACGTTCTCTATGCAACATTGCCTTTCCGGGATGAGGGCTGGCTTTTTCGCGGGTAAGCCAGCGAAGTGGc                                                                                                                                         <  2:598846/150‑1 (MQ=255)
                 cTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGACGGGGCCaaa                                                                                                                                     >  1:935944/1‑150 (MQ=255)
                   gggATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATc                                                                                                                                   <  1:338520/149‑1 (MQ=255)
                      aTCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGTGTGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGc                                                                                                                                >  1:530401/1‑150 (MQ=255)
                           gACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagc                                                                                                                           >  1:1388772/1‑150 (MQ=255)
                           gACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagc                                                                                                                           >  1:2165201/1‑150 (MQ=255)
                               gTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGct                                                                                                                       <  2:831554/150‑1 (MQ=255)
                                tACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGctc                                                                                                                      >  1:1580174/1‑150 (MQ=255)
                                   ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:614180/126‑1 (MQ=255)
                                   ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           >  2:614180/1‑126 (MQ=255)
                                   ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                         >  2:1844300/1‑128 (MQ=255)
                                   ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                         <  1:1844300/128‑1 (MQ=255)
                                   ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATAAGACAAGCAGCTCttt                                                                                                                   >  2:1377536/1‑150 (MQ=255)
                                    gAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTc                                                                                                                  >  1:2334098/1‑150 (MQ=255)
                                         tCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGcc                                                                                                             >  2:1594831/1‑150 (MQ=255)
                                         tCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGcc                                                                                                             >  1:93547/1‑150 (MQ=255)
                                                 gAGCATACCCTGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACt                                                                                                     >  2:1609450/1‑150 (MQ=255)
                                                   gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTATCCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTcc                                                                                                   <  1:1377536/150‑1 (MQ=255)
                                                    cATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                            >  2:142325/1‑124 (MQ=255)
                                                    cATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                            <  1:142325/124‑1 (MQ=255)
                                                      tACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGc                                                                                                <  2:2334098/150‑1 (MQ=255)
                                                        cccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                    >  2:1807379/1‑144 (MQ=255)
                                                        cccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                    <  1:1807379/144‑1 (MQ=255)
                                                         ccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTTTTTTTTTTTCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                             >  2:2222955/1‑150 (MQ=255)
                                                         ccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                             >  2:1880156/1‑150 (MQ=255)
                                                         ccGGAGTGCGGGTCAAGCACAACCCGCTCTTCTGCAATGTTGTCGATGCAACATTGCATGATGCGGTAGCAGGCGTGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                             <  1:44221/150‑1 (MQ=255)
                                                          cGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGATGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCgg                                                                                            >  2:1212814/1‑150 (MQ=255)
                                                              gTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCAc                                                                                        >  2:729128/1‑150 (MQ=255)
                                                                   ggTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCa                                                                                   >  1:1160717/1‑150 (MQ=255)
                                                                     ttCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAcc                                                                                 >  1:2339149/1‑150 (MQ=255)
                                                                       cAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAcccc                                                                               <  2:981541/150‑1 (MQ=255)
                                                                        aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                              <  1:829191/118‑1 (MQ=255)
                                                                        aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                              >  2:829191/1‑118 (MQ=255)
                                                                        aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                                      <  1:1363010/126‑1 (MQ=255)
                                                                        aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                                      >  2:1363010/1‑126 (MQ=255)
                                                                              aTCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                             <  1:85907/129‑1 (MQ=255)
                                                                              aTCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                             >  2:85907/1‑129 (MQ=255)
                                                                                 cACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                                <  2:1669317/107‑1 (MQ=255)
                                                                                 cACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                                >  1:1669317/1‑107 (MQ=255)
                                                                                 cACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                     >  1:1168008/1‑150 (MQ=255)
                                                                                         tgCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTc                                                             <  1:511676/150‑1 (MQ=255)
                                                                                            aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgcg                                                          >  2:1204460/1‑150 (MQ=255)
                                                                                                    cTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGt                                                  <  2:852325/150‑1 (MQ=255)
                                                                                                     tagacAACCTTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTa                                                 <  1:1176057/146‑1 (MQ=255)
                                                                                                     tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCaa                                                                                                       >  1:809176/1‑96 (MQ=255)
                                                                                                     tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCaa                                                                                                       <  2:809176/96‑1 (MQ=255)
                                                                                                     tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAc                                                                                  >  2:1166337/1‑117 (MQ=255)
                                                                                                     tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAc                                                                                  <  1:1166337/117‑1 (MQ=255)
                                                                                                        gCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                                <  2:2178574/84‑1 (MQ=255)
                                                                                                        gCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                                >  1:2178574/1‑84 (MQ=255)
                                                                                                          aaCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGAGCCAAATCAGCCAAGCAACTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGcgcc                                            >  1:1589602/1‑150 (MQ=255)
                                                                                                             aTTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCt                                         >  2:1777023/1‑150 (MQ=255)
                                                                                                                gCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGc                                      <  2:2339149/150‑1 (MQ=255)
                                                                                                                 cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGc                                                                                                                                >  1:2157753/1‑59 (MQ=255)
                                                                                                                 cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGc                                                                                                                                <  2:2157753/59‑1 (MQ=255)
                                                                                                                 cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                     >  1:2009728/1‑118 (MQ=255)
                                                                                                                 cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                     <  2:2009728/118‑1 (MQ=255)
                                                                                                                     aCCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCGTTTTGATTTTGCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGc                                 <  1:2222955/150‑1 (MQ=255)
                                                                                                                          gATGAGGGCTGGCTTTTTTTTTTTTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcg                            >  2:1196099/1‑150 (MQ=255)
                                                                                                                               gggCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCa                       <  2:1388772/150‑1 (MQ=255)
                                                                                                                                ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAATCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCGCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAg                      >  1:976909/1‑150 (MQ=255)
                                                                                                                                   tggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGccc                   <  1:1204460/150‑1 (MQ=255)
                                                                                                                                      ccccttCGCGGGTTTACCCGCGAAAGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCgcgc                                                         <  2:810586/105‑1 (MQ=39)
                                                                                                                                      ccccttCGCGGGTTTACCCGCGAAAGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCgcgc                                                         >  1:810586/5‑109 (MQ=39)
                                                                                                                                       cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAcc               >  2:550644/4‑150 (MQ=255)
                                                                                                                                         tttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc             <  2:1222136/150‑1 (MQ=255)
                                                                                                                                          ttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGcc            >  1:2345633/1‑150 (MQ=255)
                                                                                                                                             gcgGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGc         >  2:332480/1‑150 (MQ=255)
                                                                                                                                              cgGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCTTGCGCGTCAGCCCCACCGCCTGCa        <  1:2131557/150‑1 (MQ=255)
                                                                                                                                               gggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAt       <  1:1880156/150‑1 (MQ=255)
                                                                                                                                               gggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAt       >  2:792169/1‑150 (MQ=255)
                                                                                                                                               gggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCACTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCAACCTCACGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCAACCGCCTGCAt       >  1:40923/1‑150 (MQ=255)
                                                                                                                                                ggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAtg      <  2:1580174/150‑1 (MQ=255)
                                                                                                                                                 gTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAtgt     <  1:1777023/150‑1 (MQ=255)
                                                                                                                                                   aaaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGc   <  1:1196099/150‑1 (MQ=255)
                                                                                                                                                    aaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCt  >  1:1401622/1‑150 (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 < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: