Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 891,587 C→A 5.5% intergenic (+40/+40) yiaD → / ← pta OmpA/MotB domain‑containing protein/phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947891,5870CA 5.5% 243.1 / 3.3 73intergenic (+40/+40)yiaD/ptaOmpA/MotB domain‑containing protein/phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/39);  new base A (2/2);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CCAATGGTCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCTTGCGCAAACCCTGC  >  NC_002947/891450‑891732
                                                                                                                                          |                                                                                                                                                 
ccAATGGTCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTcc                                                                                                                                        <  6:163422/149‑1 (MQ=255)
ccAATGGTCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTcc                                                                                                                                        >  8:354443/1‑149 (MQ=255)
ccAATGGTCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTcc                                                                                                                                        <  1:47355/149‑1 (MQ=255)
    tGGTCGTGAGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCACGATCCcac                                                                                                                                     >  8:343905/1‑148 (MQ=255)
      gTCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACaa                                                                                                                                   <  3:252026/148‑1 (MQ=255)
       tCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAgg                                                                                                                                 <  4:350389/149‑1 (MQ=255)
       tCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAgg                                                                                                                                 <  2:199736/149‑1 (MQ=255)
       tCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAgg                                                                                                                                 <  2:193385/149‑1 (MQ=255)
          tGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCtt                                                                                                                              <  7:311060/149‑1 (MQ=255)
          tGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCtt                                                                                                                              <  2:343158/149‑1 (MQ=255)
             gcAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGc                                                                                                                            >  8:126383/1‑148 (MQ=255)
                 aaCCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCtt                                                                                                                              >  7:239986/1‑142 (MQ=255)
                 aaCCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCtt                                                                                                                              <  8:239986/142‑1 (MQ=255)
                     gccgTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTc                                                                                                                                         <  2:1936/127‑1 (MQ=255)
                     gccgTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTc                                                                                                                                         >  1:1936/1‑127 (MQ=255)
                        gtgtGGAAGTGAACCTGAAGCCGATTCCGGGCCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCctct                                                                                                                >  6:350786/1‑149 (MQ=255)
                           tgGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGc                                                                                                                    <  4:219977/142‑1 (MQ=255)
                           tgGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGc                                                                                                                    >  3:219977/1‑142 (MQ=255)
                             gAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCttt                                                                                                           <  8:267609/149‑1 (MQ=255)
                               aGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCca                                                                                                                                      <  7:40948/120‑1 (MQ=255)
                               aGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCca                                                                                                                                      >  8:40948/1‑120 (MQ=255)
                                           gCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAaagaa                                                                                                                      <  2:308841/124‑1 (MQ=255)
                                           gCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAaagaa                                                                                                                      >  1:308841/1‑124 (MQ=255)
                                           gCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCt                                                                                             >  3:276729/1‑149 (MQ=255)
                                           gCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCt                                                                                             >  2:10224/1‑149 (MQ=255)
                                                   cGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGaa                                                                                     <  5:113854/149‑1 (MQ=255)
                                                    gggTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAAt                                                                                    <  8:347843/149‑1 (MQ=255)
                                                         agcagTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGaa                                                                                     >  5:359878/1‑143 (MQ=255)
                                                         agcagTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGaa                                                                                     <  6:359878/143‑1 (MQ=255)
                                                            agTATGAGCAGCAAGGTGAGGTTCAGCAGTATCTCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGTTGGCTGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCa                                                                            <  3:51027/149‑1 (MQ=255)
                                                             gTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAaagaa                                                                                                                      >  6:15089/1‑106 (MQ=255)
                                                             gTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAaagaa                                                                                                                      <  5:15089/106‑1 (MQ=255)
                                                             gTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAg                                                                           >  2:150366/1‑149 (MQ=255)
                                                                tGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAggg                                                                         <  1:22228/148‑1 (MQ=255)
                                                                     agcaAGGTGAGGTTCAGCAGTATCGCTGATTGACCTTAAGGGCCGGTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGGTCCCACAAGGCTTGCAAAGAAGGCTCTGCTTTTAGGGCGCCACGCTGGCCTTAATCTCCGTCAGGGCTATgg                                                                   <  3:29540/149‑1 (MQ=255)
                                                                        aaGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTAt                                                                <  7:27057/149‑1 (MQ=255)
                                                                         aGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATg                                                               >  4:365312/1‑149 (MQ=255)
                                                                              aGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACg                                                           <  5:154487/148‑1 (MQ=255)
                                                                              aGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGa                                                          <  7:354443/149‑1 (MQ=255)
                                                                                gTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCt                                                                                        >  8:67754/1‑117 (MQ=255)
                                                                                gTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCt                                                                                        <  7:67754/117‑1 (MQ=255)
                                                                                 ttCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGatat                                                       <  4:276729/149‑1 (MQ=255)
                                                                                   cagcagTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATc                                                      <  7:126383/148‑1 (MQ=255)
                                                                                   cagcagTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATc                                                      <  2:251580/148‑1 (MQ=255)
                                                                                    agcagTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCg                                                     <  5:350786/148‑1 (MQ=255)
                                                                                              cTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGt                                                    <  1:62925/139‑1 (MQ=255)
                                                                                              cTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGt                                                    >  2:62925/1‑139 (MQ=255)
                                                                                              cTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTgcgc                                          <  1:10224/149‑1 (MQ=255)
                                                                                                gATTGACCTTAAGGGGCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGcc                                                                                <  2:60309/109‑1 (MQ=255)
                                                                                                gATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGcc                                                                                >  1:60309/1‑110 (MQ=255)
                                                                                                 aTTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCACGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGccc                                        <  7:343905/148‑1 (MQ=255)
                                                                                                      ccTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTcgcg                                  >  4:210679/1‑149 (MQ=255)
                                                                                                      ccTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTcgcg                                  >  3:291225/1‑149 (MQ=255)
                                                                                                       cTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCgg                                 >  6:178453/1‑149 (MQ=255)
                                                                                                        ttAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCgg                                 >  3:82337/1‑148 (MQ=255)
                                                                                                         tAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGc                                >  3:385091/1‑148 (MQ=255)
                                                                                                          aaGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGc                                <  2:24140/147‑1 (MQ=255)
                                                                                                          aaGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGc                                <  1:150366/147‑1 (MQ=255)
                                                                                                           aGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCaa                              <  4:82337/148‑1 (MQ=255)
                                                                                                            gggCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTgc                                            <  1:103564/133‑1 (MQ=255)
                                                                                                            gggCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTgc                                            >  2:103564/1‑133 (MQ=255)
                                                                                                            gggCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGt                            >  8:286208/1‑149 (MQ=255)
                                                                                                                cGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCGGGCCTGAATCGCCGTCAGGGCGATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGtt                        >  3:12285/1‑149 (MQ=255)
                                                                                                                    ttGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCa                               >  7:292031/1‑138 (MQ=255)
                                                                                                                    ttGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCa                               <  8:292031/138‑1 (MQ=255)
                                                                                                                       cGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGc                 <  4:291225/149‑1 (MQ=255)
                                                                                                                       aagccttgtgGGAGCGGGACTGCCCGCCCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTc                           <  6:421558/133‑1 (MQ=38)
                                                                                                                       aagccttgtgGGAGCGGGACTGCCCGCCCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTc                           >  5:421558/7‑139 (MQ=38)
                                                                                                                        ggCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCt                >  8:141881/1‑149 (MQ=255)
                                                                                                                          ccttgtgGGAGCGGGACTGCCCGCCCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCaa                              <  7:381602/130‑1 (MQ=255)
                                                                                                                          ccttgtgGGAGCGGGACTGCCCGCCCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCaa                              >  8:381602/4‑133 (MQ=255)
                                                                                                                           tctGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCTTgc             >  3:377281/1‑149 (MQ=255)
                                                                                                                                 gggggCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCTTGCGCaaa        >  6:170499/1‑149 (MQ=255)
                                                                                                                                    ggcgggcAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCTTGCGCAAAccc     >  8:13356/1‑149 (MQ=255)
                                                                                                                                       gggcAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCTTGCGCAAACCCTgc  <  4:385091/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CCAATGGTCGTGCGCAGAACCGCCGTGTGGAAGTGAACCTGAAGCCGATTCCGGGTCAGCAGTATGAGCAGCAAGGTCAGGTTCAGCAGTATCCCTGATTGACCTTAAGGGCCGCTTTGCGGCTCTGTGGGGGCGGGCAGTCCCGCTCCCACAAGGCTTGCAAAGAAGCCTCTGCTTTTAGCGCGCCACGCTGGCCTGAATCGCCGTCAGGGCTATGGTATGCACGATATCGTCAACCTGCGCCCCTCGCGGCAAGTCGTTCACCGGCTTGCGCAAACCCTGC  >  NC_002947/891450‑891732

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