Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,686,332 C→G 7.5% R367G (CGC→GGC)  yejA → microcin C ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,686,3320CG 7.5% 248.3 / 10.4 80R367G (CGC→GGC) yejAmicrocin C ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base G (2/4);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

CATCCGCCGCCCGGTTTTCCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACAGCGCCAGGCTTACAAGTT  >  NC_002947/4686190‑4686477
                                                                                                                                              |                                                                                                                                                  
cATCCGCCGCCCGGTTTTCCAGGATGTGGGTGTGCGCCAGGCCCTGGGCGTGATGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCt                                                                                                                                             <  2:234703/150‑1 (MQ=255)
 aTCCGCCGCCCGGTTTTCCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAGGTAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCGGGCGGccat                                                                                                                                            <  2:1173530/150‑3 (MQ=255)
  tCCGCCGCCCGGTTTTCCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCGGGCGGccatt                                                                                                                                           <  1:1167723/150‑4 (MQ=255)
          ccGGTTTTTCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGcccc                                                                                                                                   <  2:1577071/150‑1 (MQ=255)
             gTTTTCCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAAc                                                                                                                                <  2:566596/150‑1 (MQ=255)
                ttCCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAg                                                                                                                             <  2:652348/150‑1 (MQ=255)
                 cccAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGc                                                                                                                            <  2:1210101/149‑1 (MQ=255)
                   cAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCGGGCGGccattt                                                                                                                                          >  2:2152667/1‑130 (MQ=255)
                   cAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCGGGCGGccattt                                                                                                                                          <  1:2152667/134‑5 (MQ=255)
                      gATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCCACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATAGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGaaa                                                                                                                       >  1:221498/1‑150 (MQ=255)
                          tgcgtgtgcgCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCCACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATcc                                                                                                                   >  2:771869/1‑150 (MQ=255)
                            cgtgtgcgCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAACTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTc                                                                                                                 <  1:1166535/150‑1 (MQ=255)
                            cgtgtgcgCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCCACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTc                                                                                                                 <  1:771869/150‑1 (MQ=255)
                             gtgtgcgCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCg                                                                                                                >  1:2154980/1‑150 (MQ=255)
                              tgtgcgCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGa                                                                                                               >  1:1776033/1‑150 (MQ=255)
                                   gcCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAAc                                                                                                                                >  1:75661/1‑128 (MQ=255)
                                   gcCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAAc                                                                                                                                <  2:75661/128‑1 (MQ=255)
                                    ccAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCt                                                                                                                           >  1:826634/1‑132 (MQ=255)
                                    ccAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCt                                                                                                                           <  2:826634/132‑1 (MQ=255)
                                        gCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACATCCTCAACAATATCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTgcgc                                                                                                    <  1:43261/150‑1 (MQ=255)
                                           cTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGc                                                                                                 >  2:1572696/1‑150 (MQ=255)
                                             gAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAATCACTGCGCGGCaa                                                                                               >  2:2126364/1‑150 (MQ=255)
                                              aGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCaaa                                                                                              <  1:1572696/150‑1 (MQ=255)
                                                          gACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCaa                                                                                  <  1:1914640/150‑1 (MQ=255)
                                                          gACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCaa                                                                                  <  1:2144675/150‑1 (MQ=255)
                                                           aCTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAg                                                                                 >  2:1614086/1‑150 (MQ=255)
                                                             tACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTa                                                                               >  2:680877/1‑150 (MQ=255)
                                                                gAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGa                                                                                                                         >  1:381016/1‑106 (MQ=255)
                                                                gAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGa                                                                                                                         <  2:381016/106‑1 (MQ=255)
                                                                     gACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCaa                                                                                  >  1:1321043/1‑139 (MQ=255)
                                                                     gACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCaa                                                                                  <  2:1321043/139‑1 (MQ=255)
                                                                        caacaaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGt                                                                    >  2:1794715/1‑150 (MQ=255)
                                                                         aacaaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTc                                                                   >  2:729194/1‑150 (MQ=255)
                                                                          acaaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCaaa                                                                                              >  2:514625/1‑122 (MQ=255)
                                                                          acaaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCaaa                                                                                              <  1:514625/122‑1 (MQ=255)
                                                                          acaaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAg                                                                      >  2:819967/1‑146 (MQ=255)
                                                                          acaaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAg                                                                      <  1:819967/146‑1 (MQ=255)
                                                                            aaGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTc                                                                <  1:214362/150‑1 (MQ=255)
                                                                               cagcTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCcac                                                             >  2:2239769/1‑150 (MQ=255)
                                                                                agcTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCcaca                                                            >  1:1592848/1‑150 (MQ=255)
                                                                                       cAACGGCGCCTACACCCGTACCGGCAGCGACTTCGAGAACTCGGAAATGGCCGTCCGCGGCCTGCCCAGCCCCAACGAGCTGAAGATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACGCGGt                                                     >  2:1091521/1‑150 (MQ=255)
                                                                                        aaCGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCGGGCGGccat                                                                                                                                            <  1:1816152/63‑3 (MQ=255)
                                                                                        aaCGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCGGGCGGccat                                                                                                                                            >  2:1816152/1‑61 (MQ=255)
                                                                                        aaCGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAg                                                                      >  2:1864105/1‑132 (MQ=255)
                                                                                        aaCGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAg                                                                      <  1:1864105/132‑1 (MQ=255)
                                                                                         aCGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCa                                                   >  1:1618675/1‑150 (MQ=255)
                                                                                                cTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTAACGAGGTCTTCCACAACCCGGTCAGTGACgg                                            >  1:222975/1‑150 (MQ=255)
                                                                                                 tACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGc                                           >  1:131778/1‑150 (MQ=255)
                                                                                                   cacCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGAc                                              >  1:200212/1‑145 (MQ=255)
                                                                                                   cacCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGAc                                              <  2:200212/145‑1 (MQ=255)
                                                                                                    acCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGcc                                                                                                                                          <  2:1777542/53‑1 (MQ=255)
                                                                                                    acCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGcc                                                                                                                                          >  1:1777542/1‑53 (MQ=255)
                                                                                                    acCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAACGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGCCAGTGACGGCAGc                                        >  2:2084990/1‑150 (MQ=255)
                                                                                                    acCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACg                                             <  2:1726862/145‑1 (MQ=255)
                                                                                                    acCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACg                                             >  1:1726862/1‑145 (MQ=255)
                                                                                                    acCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGACCGCCGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGc                                        >  1:1760332/1‑150 (MQ=255)
                                                                                                     cccGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCg                                       <  2:131778/150‑1 (MQ=255)
                                                                                                      ccGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCaa                                                                                  >  1:510664/1‑106 (MQ=255)
                                                                                                      ccGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCaa                                                                                  <  2:510664/106‑1 (MQ=255)
                                                                                                      ccGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTa                                                                           <  2:690219/113‑1 (MQ=255)
                                                                                                      ccGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTa                                                                           >  1:690219/1‑113 (MQ=255)
                                                                                                      ccGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACACTCACGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCgg                                      >  1:761290/1‑150 (MQ=255)
                                                                                                         tACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCAt                                   <  1:1164956/150‑1 (MQ=255)
                                                                                                          aCCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGcc                                                                                         <  2:77127/95‑1 (MQ=255)
                                                                                                          aCCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGcc                                                                                         >  1:77127/1‑95 (MQ=255)
                                                                                                            cGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGAt                                >  2:90580/1‑150 (MQ=255)
                                                                                                              gCAGCTACTTCGAAAATTCGGAAATGGCCGCCCGGGGCCTGCCCAGCACCAACGATCTGAAAATCCTCGAACCACTTCGCGGCAAAGTGCCCGAGTAATTATTTACCGACGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATcc                              <  2:221498/150‑1 (MQ=255)
                                                                                                                aGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCcgc                            >  2:2377281/1‑150 (MQ=255)
                                                                                                                    aCTTCGAGACCTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGc                        <  2:300162/150‑1 (MQ=255)
                                                                                                                      ttCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCaa                      >  2:120470/1‑150 (MQ=255)
                                                                                                                      ttCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCaa                      >  2:670691/1‑150 (MQ=255)
                                                                                                                        cGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACa                    <  2:1776033/150‑1 (MQ=255)
                                                                                                                         gataaCTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACgg                                            >  1:349547/4‑125 (MQ=255)
                                                                                                                         gagaACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACgg                                            <  2:349547/125‑1 (MQ=255)
                                                                                                                         gagaACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACAg                   >  1:2334611/1‑150 (MQ=255)
                                                                                                                                 ggAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGt                                                                                >  2:1775832/1‑81 (MQ=255)
                                                                                                                                 ggAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGt                                                                                <  1:1775832/81‑1 (MQ=255)
                                                                                                                                    aaTGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACAGCGCCAGGCTTa        <  1:120470/150‑1 (MQ=255)
                                                                                                                                       ggCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACTGCGCCAGGCTTACaa     <  2:761290/150‑1 (MQ=255)
                                                                                                                                        gccgccCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACAGCGCCAGGCTTACAAg    <  2:1850262/150‑1 (MQ=255)
                                                                                                                                         ccgcccgcGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAgc                                                                                            >  1:79375/1‑60 (MQ=255)
                                                                                                                                         ccgcccgcGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAgc                                                                                            <  2:79375/61‑2 (MQ=255)
                                                                                                                                          cgcccgcGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACAGCGCCAGGCTTACAAgtt  >  2:440993/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CATCCGCCGCCCGGTTTTCCAGGATGTGCGTGTGCGCCAGGCCCTGAGCCTGCTGCTGGACTACGAGTGGACCAACAAGCAGCTGTTCAACGGCGCCTACACCCGTACCGGCAGCTACTTCGAGAACTCGGAAATGGCCGCCCGCGGCCTGCCCAGCCCCAACGAGCTGAAAATCCTCGAACCACTGCGCGGCAAAGTGCCCGAGCAAGTATTTACCGAGGTCTTCCACAACCCGGTCAGTGACGGCAGCGGCATGATCCGCGAGCAACAGCGCCAGGCTTACAAGTT  >  NC_002947/4686190‑4686477

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: