Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,227 T→A 9.5% intergenic (+34/+34) ilvD → / ← PP_5129 dihydroxy‑acid dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,2270TA 9.5% 181.3 / 8.9 63intergenic (+34/+34)ilvD/PP_5129dihydroxy‑acid dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (2/4);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCCAGCCGAGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTACCCCTCC  >  NC_002947/5851079‑5851374
                                                                                                                                                    |                                                                                                                                                   
gccagccGAGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTg                                                                                                                                                    >  1:379803/1‑150 (MQ=255)
 ccagccGAGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTgg                                                                                                                                                   >  1:518790/1‑150 (MQ=255)
       gAGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCCGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCg                                                                                                                                             <  1:1110911/150‑1 (MQ=255)
        aGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCAGGGCGCCgg                                                                                                                                            <  2:1146740/150‑1 (MQ=255)
          gtccgtccACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCAGGGCGCCGGtt                                                                                                                                          <  1:1023486/150‑1 (MQ=255)
                ccACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCAGGGCGCCGGtttttttt                                                                                                                                    <  2:699804/150‑1 (MQ=255)
                  aCGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGa                                                                                                                                  >  2:1291563/1‑150 (MQ=255)
                  aCGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGa                                                                                                                                  <  1:371882/150‑1 (MQ=255)
                    gCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGa                                                                                                                                  >  2:655664/1‑148 (MQ=255)
                    gCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGa                                                                                                                                  <  1:655664/148‑1 (MQ=255)
                      aaGGTAACCCCCGCGCTCAAGGCCGATGCCCTGCGGTCCACCCGTGCAGACCAGGGTGCACTCCACTGCGGCTGGATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGt                                                                                                                              <  1:738120/150‑1 (MQ=255)
                      aaGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGt                                                                                                                              >  1:685906/1‑150 (MQ=255)
                             ccaccGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGa                                                                                                                       <  2:6148/150‑1 (MQ=255)
                             ccaccGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGa                                                                                                                       <  1:713269/150‑1 (MQ=255)
                              caccGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCAGGGCGcc                                                                                                                                              <  2:574479/126‑1 (MQ=255)
                              caccGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCAGGGCGcc                                                                                                                                              >  1:574479/1‑126 (MQ=255)
                               cgatctcTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACACGGATGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACg                                                                                                                     <  1:513038/144‑1 (MQ=255)
                               accGCGCTCAAGTCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCTTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGTCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACg                                                                                                                     <  1:1245086/150‑1 (MQ=255)
                                      tCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTg                                                                                                              >  1:1290507/1‑150 (MQ=255)
                                      tCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTg                                                                                                              >  1:1290639/1‑150 (MQ=255)
                                         aGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAg                                                                                                           >  1:1051457/1‑150 (MQ=255)
                                              tATGCCCTGCTGGCCACCAGTGCCGACAAAGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAgg                                                                                                      >  2:1138834/1‑150 (MQ=255)
                                                    ctgctgGCAACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGCCGATTACTTAAGCGAGGTGGTTg                                                                                                >  2:660656/1‑150 (MQ=255)
                                                     tgctgGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGCAACCCCTGGAAGGACTATGAGGTGCTTTTGGTGGGAAATAAGATCCGGCTCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTgg                                                                                               <  2:518790/150‑1 (MQ=255)
                                                        tgGCCACCAGTGCCGACAAGGGTGCTGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTc                                                                                            <  2:297763/150‑1 (MQ=255)
                                                          gCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCgg                                                                                          >  2:711227/1‑150 (MQ=255)
                                                           ccaccaGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGa                                                                                         <  1:491672/150‑1 (MQ=255)
                                                            caccaGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGa                                                                                         >  2:135853/1‑149 (MQ=255)
                                                            caccaGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGa                                                                                         <  1:135853/149‑1 (MQ=255)
                                                                aGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGCAGCGAGGTGGCTGGGTCGGATCCaa                                                                                    <  2:379803/150‑1 (MQ=255)
                                                                       aCAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGa                                                                             <  1:343293/150‑1 (MQ=255)
                                                                        cAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGAc                                                                            <  2:1301746/150‑1 (MQ=255)
                                                                         aaGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCCGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACa                                                                           >  1:1167817/1‑150 (MQ=255)
                                                                              tgcggtgcgTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAg                                                                                                       >  2:460714/1‑117 (MQ=255)
                                                                              tgcggtgcgTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAg                                                                                                       <  1:460714/117‑1 (MQ=255)
                                                                                cggtgcgTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTg                                                                    <  1:1138834/150‑1 (MQ=255)
                                                                                    gcgTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGgaga                                                                <  2:685906/150‑1 (MQ=255)
                                                                                       tAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCa                                                                                                                 <  2:343700/98‑1 (MQ=255)
                                                                                       tAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCa                                                                                                                 >  1:343700/1‑98 (MQ=255)
                                                                                          ccaGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGAtt                                                          >  2:12251/3‑150 (MQ=255)
                                                                                          cAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATc                                                                                       <  1:850969/121‑1 (MQ=255)
                                                                                          cAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATc                                                                                       >  2:850969/1‑121 (MQ=255)
                                                                                           aaGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTc                                                                                            >  2:1313656/1‑115 (MQ=255)
                                                                                           aaGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTc                                                                                            <  1:1313656/115‑1 (MQ=255)
                                                                                            aGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTcc                                                        >  2:297305/1‑150 (MQ=255)
                                                                                              gCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAg                                                                                                       >  1:565481/1‑101 (MQ=255)
                                                                                              gCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAg                                                                                                       <  2:565481/101‑1 (MQ=255)
                                                                                                aaTGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGt                                                    <  2:1290639/150‑1 (MQ=255)
                                                                                                aaTGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGt                                                    >  2:817561/1‑150 (MQ=255)
                                                                                                  tGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGc                                                  <  1:297305/150‑1 (MQ=255)
                                                                                                             ctctGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGa                                       >  1:140926/1‑150 (MQ=255)
                                                                                                              tctGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGaa                                      <  1:660656/150‑1 (MQ=255)
                                                                                                                tGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAAcc                                    <  2:976045/150‑2 (MQ=255)
                                                                                                                     tGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGCCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCTAAGTGCCTCGTCATCGAACAGTTTg                               <  1:1291563/150‑1 (MQ=255)
                                                                                                                               gTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGag                     <  2:140926/150‑1 (MQ=255)
                                                                                                                                  aaaTAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGagcag                  <  2:1167817/150‑1 (MQ=255)
                                                                                                                                     tAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTTAAAAACACTAAGAGTATTGAGGGGTCTGGGGCGTATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCgg               >  1:27217/1‑150 (MQ=255)
                                                                                                                                      aaGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCgg               <  2:1233284/149‑1 (MQ=255)
                                                                                                                                      aaGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCgg               >  1:1233284/1‑149 (MQ=255)
                                                                                                                                           ccGGCGCCCAGGGCGCCGGTTTTTTTTGATTGTTCATTTACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTa         >  2:688648/1‑150 (MQ=255)
                                                                                                                                            cGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGagca                   <  2:17582/139‑1 (MQ=255)
                                                                                                                                            cGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGagca                   >  1:17582/1‑139 (MQ=255)
                                                                                                                                                  ccTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTACCcctcc  <  1:817561/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCAGCCGAGGTCCGTCCACGCAAGGTAACCACCGCGCTCAAGGCCTATGCCCTGCTGGCCACCAGTGCCGACAAGGGTGCGGTGCGTAACAAGGCAATGCTGGAAGGCCTCTGAGGTGCTTTTAGCGTGAAATAAGAACCGGCGCCCTGGGCGCCGGTTTTTTTTGATTGTTCATTGACGGTCAGTGAAGCGAGGTGGCTGGGTCGGATCGAAACATTGACAGCTTTTGGAGAGCGATTCCAAGTGCCTCGTCATCGAACAGTTTGCTTTTTCGAGCAGCGGCACCTACCCCTCC  >  NC_002947/5851079‑5851374

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