Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,080,964 A→T 10.4% 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,9640AT10.4% 260.3 / 20.3 96intergenic (‑144/+19)amaC/dpkAD‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (47/39);  new base T (7/3);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

GAACACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGC  >  NC_002947/4080815‑4081112
                                                                                                                                                     |                                                                                                                                                    
gAACACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTaa                                                                                                                                                      >  2:1626042/1‑150 (MQ=255)
  acacAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAc                                                                                                                                                    >  2:2261958/1‑150 (MQ=255)
    acaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAccc                                                                                                                                                  >  1:1583123/1‑150 (MQ=255)
       gCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcgcg                                                                                                                                               <  1:334898/150‑1 (MQ=255)
         gTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGaa                                                                                                                                             <  2:2321075/150‑1 (MQ=255)
           gCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                           <  1:1382916/150‑1 (MQ=255)
                  cTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCaaa                                                                                                                                    <  1:2261958/150‑1 (MQ=255)
                   ttGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                           >  1:939578/1‑142 (MQ=255)
                   ttGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                           <  2:939578/142‑1 (MQ=255)
                      gATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGAGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGAAGAAAGCTTGGCAGTTAAACCCGCGAAGGGGCCAAATCAg                                                                                                                                >  2:1380268/1‑150 (MQ=255)
                            gACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCTGGATGAGGGCTGGCTTTTTTGCGGGTAAACCCGCGAAGGGGCCAAAt                                                                                                                                   >  1:1603804/1‑141 (MQ=255)
                            gACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCTGGATGAGGGCTGGCTTTTTTGCGGGTAAACCCGCGAAGGGGCCAAAt                                                                                                                                   <  2:1603804/141‑1 (MQ=255)
                             aCGGTACGGAAAGGCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagca                                                                                                                         <  1:1691809/150‑1 (MQ=255)
                               ggTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagcagc                                                                                                                       <  2:1677095/150‑1 (MQ=255)
                                gTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCGCTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGct                                                                                                                      >  2:744846/1‑150 (MQ=255)
                                    ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCa                                                                                                                                 <  1:935748/135‑1 (MQ=255)
                                    ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCa                                                                                                                                 >  2:935748/1‑135 (MQ=255)
                                    ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCttt                                                                                                                  >  1:2355832/1‑150 (MQ=255)
                                    ggAAAGTCAGTGTTGAGCATACCAGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTAAGCAACATTGCATGACCAGGATAAGGGCTGGCTTTTTCACTGGTAAACCCGCGCACGGGCCAAATCAGCCAAGAAGCTAttt                                                                                                                  >  2:181024/1‑150 (MQ=255)
                                         gTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGAAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCAGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                             <  1:788877/150‑1 (MQ=255)
                                            aGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCtt                                                                                                          <  2:1583123/150‑1 (MQ=255)
                                                ttGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCaa                                                                                                      >  2:481578/1‑150 (MQ=255)
                                                   aGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                   <  2:2263103/150‑1 (MQ=255)
                                                    gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgggg                                                                                                                                         >  2:366887/1‑111 (MQ=255)
                                                    gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgggg                                                                                                                                         <  1:366887/111‑1 (MQ=255)
                                                    gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTcc                                                                                                  >  1:587278/1‑150 (MQ=255)
                                                          ccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCAGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                            >  1:268441/1‑150 (MQ=255)
                                                            ggAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCaa                                                                                                                            <  2:1437302/116‑1 (MQ=255)
                                                            ggAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCaa                                                                                                                            >  1:1437302/1‑116 (MQ=255)
                                                             gAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTc                                                                                         >  1:1250625/1‑150 (MQ=255)
                                                                tGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCACCTCCTGCTCGGTCAcc                                                                                      >  2:1793126/1‑150 (MQ=255)
                                                                    ggTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCa                                                                                  <  1:528127/150‑1 (MQ=255)
                                                                     gTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAc                                                                                 >  2:1033039/1‑150 (MQ=255)
                                                                       tCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAg                                                                                                                                            <  1:1381681/89‑1 (MQ=255)
                                                                       tCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAg                                                                                                                                            >  2:1381681/1‑89 (MQ=255)
                                                                           gAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTCAACCCGCGAAGGGGCCAAACAGCCAAGCAGGCCCTTCCGGCCTTTTAACCCCTGGTCCGCCCCCGCCAcccccccc                                                                          >  1:2372571/1‑150 (MQ=255)
                                                                               aTCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctcc                                                                       <  2:818910/150‑1 (MQ=255)
                                                                                tCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctccc                                                                      >  1:2264371/1‑149 (MQ=255)
                                                                                  cACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGTTCACCGCCACCCCCTCCTCCTCg                                                                    <  1:1066072/150‑1 (MQ=255)
                                                                                   aCTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCgg                                                                   >  1:245512/1‑150 (MQ=255)
                                                                                    cTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGACGGCCGAGTTGTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGc                                                                  <  1:744846/150‑1 (MQ=255)
                                                                                     tATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGcc                                                                 >  1:650702/1‑150 (MQ=255)
                                                                                            cAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGCAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          <  2:2031370/150‑1 (MQ=255)
                                                                                             aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAACCCGGCGAGGGGGCCAAATCGGCCAAGCAGCGTTTTCGGGCCTTTCAACCCCTGCCCGGGCACCGCAGCCCCCTCCTCCCCAGCcccccgcgg                                                         >  2:657663/1‑145 (MQ=255)
                                                                                             aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgcg                                                         >  2:2399194/1‑150 (MQ=255)
                                                                                             aaGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgcg                                                         >  1:1054294/1‑150 (MQ=255)
                                                                                                    tCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACgg                                                  >  1:1622184/1‑150 (MQ=255)
                                                                                                    tCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACgg                                                  >  2:2391279/1‑150 (MQ=255)
                                                                                                      tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCCCGGCCACCTCGCGCTCACGGTa                                                <  2:587278/150‑1 (MQ=255)
                                                                                                      tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTCCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCt                                                                                              >  2:79717/1‑104 (MQ=255)
                                                                                                      tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTCCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCt                                                                                              <  1:79717/104‑1 (MQ=255)
                                                                                                       aTGCAACATTGCAGGGCGGGGATGAGGGCTGGCTGTTGCGCGGGTAATCCCGCGAAGGGGCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgc                                              <  2:2372571/150‑1 (MQ=255)
                                                                                                        tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                           <  1:1926504/73‑1 (MQ=255)
                                                                                                        tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                           >  2:1926504/1‑73 (MQ=255)
                                                                                                        tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCAc                                                    <  2:1324285/144‑1 (MQ=255)
                                                                                                        tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCAc                                                    >  1:1324285/1‑144 (MQ=255)
                                                                                                         gCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgcg                                             >  1:27534/1‑150 (MQ=255)
                                                                                                           aaCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGcgcc                                           >  2:370247/1‑150 (MQ=255)
                                                                                                           aaCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAACCAGCCAAGCAGCTCTTTCAGGCCTTTCAACCCCTGCTCGGTAACCGCCACCCCGTCCTCCCCGACCACCTCGCGCTCCCGGTCGcgcc                                           >  2:803474/1‑150 (MQ=255)
                                                                                                            aCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGcgccg                                          >  1:111064/1‑150 (MQ=255)
                                                                                                             cATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCCCCCCGGCCACCTCGCGCTCACGGTAGcgccgc                                         <  1:1626042/150‑1 (MQ=255)
                                                                                                              aTTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCt                                        >  1:353460/1‑150 (MQ=255)
                                                                                                               ttGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTc                                       <  2:1250625/150‑1 (MQ=255)
                                                                                                                  cATGACCGTGATGAGGGCTGGCTTTTTCGCGGGTAAAGCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCCCCCCCTCCTCCTCGGCCCCCTCCCGCTCACGGTAGCGCCGCGCGcc                                    <  2:268441/150‑1 (MQ=255)
                                                                                                                  cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGcc                                    <  2:245512/150‑1 (MQ=255)
                                                                                                                  cATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGcc                                    >  1:820442/1‑150 (MQ=255)
                                                                                                                   aTGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGccc                                   >  1:1674869/1‑150 (MQ=255)
                                                                                                                   aTGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCCCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGccc                                   >  2:1288845/1‑150 (MQ=255)
                                                                                                                       ccGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                    >  1:882842/1‑113 (MQ=255)
                                                                                                                       ccGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                    <  2:882842/113‑1 (MQ=255)
                                                                                                                       ccGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCa                               >  2:2027584/1‑150 (MQ=255)
                                                                                                                       ccGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCa                               >  1:850276/1‑150 (MQ=255)
                                                                                                                            aTGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                          <  2:2449508/150‑1 (MQ=255)
                                                                                                                            aTGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                          <  2:2339236/150‑1 (MQ=255)
                                                                                                                                 ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcct                                                                            <  1:601300/95‑1 (MQ=255)
                                                                                                                                 ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcct                                                                            >  2:601699/1‑95 (MQ=255)
                                                                                                                                 ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcct                                                                            <  1:601699/95‑1 (MQ=255)
                                                                                                                                 ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcct                                                                            >  2:601300/1‑95 (MQ=255)
                                                                                                                                  gctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGc                    >  1:2097938/1‑150 (MQ=255)
                                                                                                                                     ggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGcccc                 <  2:820442/150‑1 (MQ=255)
                                                                                                                                        tttttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          <  1:1004829/106‑1 (MQ=255)
                                                                                                                                        tttttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          >  2:1004829/1‑106 (MQ=255)
                                                                                                                                        cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAcc              >  1:842939/4‑150 (MQ=255)
                                                                                                                                         ccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTc                                                            >  1:734582/3‑103 (MQ=255)
                                                                                                                                         ccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTc                                                            <  2:734582/101‑1 (MQ=255)
                                                                                                                                          cttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                            >  2:220090/2‑70 (MQ=38)
                                                                                                                                          cttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                            <  1:220090/69‑1 (MQ=37)
                                                                                                                                          cttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc            >  1:84206/2‑150 (MQ=255)
                                                                                                                                          cttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc            >  2:1839048/2‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCt                                                                                                           >  1:2133807/1‑54 (MQ=25)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCt                                                                                                           <  2:2133807/54‑1 (MQ=25)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGcc           >  1:2421436/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGcc           <  2:1674869/150‑1 (MQ=255)
                                                                                                                                              gcgGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGc        >  1:172145/1‑150 (MQ=255)
                                                                                                                                                gggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCCGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGTGAACACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAt      <  1:1241100/150‑1 (MQ=255)
                                                                                                                                                   tAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGcgccgc                                         <  1:2045240/112‑1 (MQ=255)
                                                                                                                                                   tAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGcgccgc                                         >  2:2045240/1‑112 (MQ=255)
                                                                                                                                                    aaaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCATCCTCGGCCACATCGCGCTGACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGc  >  1:281069/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGC  >  NC_002947/4080815‑4081112

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