Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,076,300 T→G 5.7% L338R (CTG→CGG) ‡ mreB → cell wall structural actin‑like protein MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,076,3000TG 5.7% 412.0 / 9.6 136L338R (CTG→CGG) ‡mreBcell wall structural actin‑like protein MreB
Reads supporting (aligned to +/- strand):  ref base T (69/59);  new base G (4/4);  total (73/63)
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.92e-01

Only 100 of 137 total aligned reads displayed.

CTGGCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGT  >  NC_002947/1076152‑1076449
                                                                                                                                                    |                                                                                                                                                     
cTGGCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTg                                                                                                                                                      >  1:602950/1‑150 (MQ=255)
   gCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCg                                                                                                                                                   >  1:978415/1‑150 (MQ=255)
     gACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGt                                                                                                                                                 >  1:1121738/1‑150 (MQ=255)
     gACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGt                                                                                                                                                 >  1:621728/1‑150 (MQ=255)
       cgaTCGTCCAGGCGGTGAAGAGCGCGCTCGAGCAAGCGCCGCCGGAACTCGCCTCGGACATCGCCGGGCGCGGGCTGGTGCTGACCGGTGGTGGTGCGCTGGTGCGTCAGTTGGAGAAACTGCTGTCCCAGGAAACCGGCATGCCGGTGa                                                                                                                                               <  2:882962/148‑1 (MQ=255)
       ccATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGa                                                                                                                                               <  1:625269/150‑1 (MQ=255)
        cATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGAt                                                                                                                                              >  1:72557/1‑150 (MQ=255)
          tCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTg                                                                                                                                            >  1:2367988/1‑150 (MQ=255)
           cGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCAGGCCGGtttcctg                                                                                                                                           <  2:452499/150‑7 (MQ=255)
                 ggCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTTATTTTCGCCGa                                                                                                                                     >  2:284083/1‑150 (MQ=255)
                 ggCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGa                                                                                                                                     <  2:557238/150‑1 (MQ=255)
                    ggTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGa                                                                                                                                  >  1:2378012/1‑150 (MQ=255)
                     gTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGAc                                                                                                                                 >  2:2141396/1‑150 (MQ=255)
                      tCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCAGGCCGGtttc                                                                                                                                              <  1:2004486/136‑4 (MQ=255)
                      tCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCAGGCCGGtttc                                                                                                                                              >  2:2004486/1‑133 (MQ=255)
                               cccTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCt                                                                                                                       <  1:776073/150‑1 (MQ=255)
                                ccTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTg                                                                                                                      >  1:1467536/1‑150 (MQ=255)
                                  tCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGa                                                                                                                          <  1:1631313/144‑1 (MQ=255)
                                  tCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGa                                                                                                                          >  2:1631313/1‑144 (MQ=255)
                                  tCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCg                                                                                                                    >  2:1930787/1‑150 (MQ=255)
                                       cAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGAGAGATTGTGACAAACTGCTGGCCCAGGAAACAGGCGGGCCGGTGATGGTCGCCGAGGACCCGCTGACCTGCGTGGcc                                                                                                               >  1:539666/1‑150 (MQ=255)
                                        aaTCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGccc                                                                                                              >  2:2114694/1‑150 (MQ=255)
                                            gccgccGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTgg                                                                                                          >  1:462115/1‑150 (MQ=255)
                                            gccgccGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTgg                                                                                                          <  1:2377292/150‑1 (MQ=255)
                                                ccgGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTggcggc                                                                                                      <  1:966267/150‑1 (MQ=255)
                                                  ggAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGc                                                                                                                                        >  1:1040266/1‑114 (MQ=255)
                                                  ggAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGc                                                                                                                                        <  2:1040266/114‑1 (MQ=255)
                                                  ggAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGCCCCGCTGACCTGCGTCGCCCGTGGcggcgg                                                                                                    <  1:2257106/150‑1 (MQ=255)
                                                  ggAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGcggcgg                                                                                                    <  2:839582/150‑1 (MQ=255)
                                                  ccgatctgCCTCGGACAAAGCCGAGCGCCGACTCTTGTTGACCGGTGTTGTTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGcggcgg                                                                                                    <  2:872231/143‑1 (MQ=255)
                                                    agcTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACCACATGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTc                                                                                                  <  1:161457/148‑1 (MQ=255)
                                                      cTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTcgc                                                                                                <  1:1857872/150‑1 (MQ=255)
                                                       tCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGAGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTcgcg                                                                                               >  2:1925508/1‑150 (MQ=255)
                                                         gCCTCGGACATCGCCGAGCGCGGCGTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGcc                                                                                             >  1:508690/1‑150 (MQ=255)
                                                                  aTCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAgatg                                                                                    <  2:651803/150‑1 (MQ=255)
                                                                   tCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAgatga                                                                                   <  2:2347890/150‑1 (MQ=255)
                                                                    cGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCTCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAgatgat                                                                                  <  1:1249948/150‑1 (MQ=255)
                                                                     gCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGatgatg                                                                                 <  2:539666/150‑1 (MQ=255)
                                                                              ggCCTGGTGTTGACCGGGGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAc                                                                        >  2:1556764/1‑150 (MQ=255)
                                                                               gCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGa                                                                                                                                     <  2:621321/88‑1 (MQ=255)
                                                                                ccTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgc                                                                      <  1:1443643/150‑1 (MQ=255)
                                                                                ccTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgc                                                                      >  2:1712360/1‑150 (MQ=255)
                                                                                 cTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgcg                                                                     >  2:189593/1‑150 (MQ=255)
                                                                                  tGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGa                                                                    >  2:47183/1‑150 (MQ=255)
                                                                                  tGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGa                                                                    <  1:2141396/150‑1 (MQ=255)
                                                                                  tGGTGTTGACCGGGGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGa                                                                    >  2:145409/1‑150 (MQ=255)
                                                                                    gtgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgcgcag                                                                  <  1:1892599/150‑4 (MQ=255)
                                                                                    gtgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATg                                                                  >  2:1358666/1‑150 (MQ=255)
                                                                                     tgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCAGGCCGGtt                                                                                                                                                >  1:1296062/1‑70 (MQ=255)
                                                                                     tgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCAGGCCGGtt                                                                                                                                                <  2:1296062/71‑2 (MQ=255)
                                                                                      gtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGAccc                                                                                                                               <  1:681116/87‑1 (MQ=255)
                                                                                       ttGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGAc                                                               <  2:2367988/150‑1 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGAtt                                                                                                                                             >  2:1816538/1‑70 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCt                                                             <  1:2114694/150‑1 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCATGGAGATGATGGACAAGCACGCGATGGACCt                                                             >  2:410157/1‑150 (MQ=255)
                                                                                          aCCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTg                                                            >  2:209077/1‑150 (MQ=255)
                                                                                           ccGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGc                                                           <  1:1930787/150‑1 (MQ=255)
                                                                                               tggtggtgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGc                                                                                                                     <  2:1528063/88‑1 (MQ=255)
                                                                                                 gtggtgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCa                                                     <  1:1349087/150‑1 (MQ=255)
                                                                                                   ggtgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAATCGGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGc                                                   <  1:747139/150‑1 (MQ=255)
                                                                                                     tgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGa                                                 <  1:145409/150‑1 (MQ=255)
                                                                                                     tgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGa                                                 <  1:1358666/150‑1 (MQ=255)
                                                                                                     tgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCAGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGa                                                 <  1:1983773/150‑1 (MQ=255)
                                                                                                      gcgcTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAg                                                <  2:1449419/150‑1 (MQ=255)
                                                                                                       cgcTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGt                                               <  1:2036658/150‑1 (MQ=255)
                                                                                                        gctgctGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGAc                                                               >  1:42594/1‑133 (MQ=255)
                                                                                                        gctgctGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGAc                                                               <  2:42594/133‑1 (MQ=255)
                                                                                                          tgctgcGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAgcggg                                            >  1:1771447/1‑146 (MQ=255)
                                                                                                                gTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgct                                      <  2:72557/150‑1 (MQ=255)
                                                                                                                 tGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgctc                                     >  1:376501/1‑150 (MQ=255)
                                                                                                                     cTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTg                                 >  2:959996/1‑150 (MQ=255)
                                                                                                                     cTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTg                                 <  1:1607024/150‑1 (MQ=255)
                                                                                                                     cTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGAGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGGGATGATGGACAAGCACGCGAGGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTg                                 >  1:1447436/1‑150 (MQ=255)
                                                                                                                      tCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGt                                >  1:1341980/1‑150 (MQ=255)
                                                                                                                      tCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGt                                >  2:237936/1‑150 (MQ=255)
                                                                                                                      tCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGt                                >  1:300822/1‑150 (MQ=255)
                                                                                                                      tCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGt                                >  1:1027611/1‑150 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCtctc                                                       <  2:2323924/125‑1 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCtctc                                                       >  1:2323924/1‑125 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTc                              >  2:2250117/1‑150 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTc                              <  2:462115/150‑1 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTc                              >  2:775407/1‑150 (MQ=255)
                                                                                                                         aCAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCa                             >  1:2218072/1‑150 (MQ=255)
                                                                                                                           aaaCTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACg                           >  1:2098246/1‑150 (MQ=255)
                                                                                                                            aaCTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                          <  1:775407/150‑1 (MQ=255)
                                                                                                                                gctgGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTggcggc                                                                                                      <  1:760591/70‑1 (MQ=255)
                                                                                                                                gctgGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTggcggc                                                                                                      >  2:760591/1‑70 (MQ=255)
                                                                                                                                    gCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCgg                  <  1:189593/150‑1 (MQ=255)
                                                                                                                                     cccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGt                 <  1:173873/150‑1 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                          >  1:1519809/1‑140 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGtt                >  2:586482/1‑150 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGtt                >  2:1179775/1‑150 (MQ=255)
                                                                                                                                       cAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAg                                                >  2:1425050/1‑117 (MQ=255)
                                                                                                                                       cAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGttt               >  1:1906107/1‑150 (MQ=255)
                                                                                                                                        aGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTa              >  2:280421/1‑150 (MQ=255)
                                                                                                                                        aGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTa              >  2:967856/1‑150 (MQ=255)
                                                                                                                                          gAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACa            <  1:284083/150‑1 (MQ=255)
                                                                                                                                                  ccTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAgtac    >  2:1518253/1‑150 (MQ=255)
                                                                                                                                                   cTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGa                                                 >  2:1389339/1‑104 (MQ=255)
                                                                                                                                                   cTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGa                                                 <  1:1389339/104‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTACGT  >  NC_002947/1076152‑1076449

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: