Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,076,299 C→A 5.7% L338M (CTG→ATG) ‡ 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,2990CA 5.7% 433.0 / 8.0 133L338M (CTG→ATG) ‡mreBcell wall structural actin‑like protein MreB
Reads supporting (aligned to +/- strand):  ref base C (67/57);  new base A (3/5);  total (71/62)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

CTGGCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTAC  >  NC_002947/1076152‑1076447
                                                                                                                                                   |                                                                                                                                                    
cTGGCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTg                                                                                                                                                    >  1:602950/1‑150 (MQ=255)
   gCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCg                                                                                                                                                 <  1:2352432/150‑1 (MQ=255)
   gCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCg                                                                                                                                                 >  1:978415/1‑150 (MQ=255)
     gACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGt                                                                                                                                               >  1:1121738/1‑150 (MQ=255)
       cgaTCGTCCAGGCGGTGAAGAGCGCGCTCGAGCAAGCGCCGCCGGAACTCGCCTCGGACATCGCCGGGCGCGGGCTGGTGCTGACCGGTGGTGGTGCGCTGGTGCGTCAGTTGGAGAAACTGCTGTCCCAGGAAACCGGCATGCCGGTGa                                                                                                                                             <  2:882962/148‑1 (MQ=255)
       ccATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGa                                                                                                                                             <  1:625269/150‑1 (MQ=255)
         aTCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACGCGCCTCGGACATTGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGAtt                                                                                                                                           >  2:1892599/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)
                      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                                                                                                                        >  2:1631313/1‑144 (MQ=255)
                                  tCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGa                                                                                                                        <  1:1631313/144‑1 (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)
                                                  ggAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGc                                                                                                                                      <  2:1040266/114‑1 (MQ=255)
                                                  ggAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGc                                                                                                                                      >  1:1040266/1‑114 (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)
                                                    aaCTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTc                                                                                                <  1:2168346/150‑1 (MQ=255)
                                                       tCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTcgcg                                                                                             >  2:1443643/1‑150 (MQ=255)
                                                          ccTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGtagatcgg                                                                                          >  1:872231/1‑143 (MQ=255)
                                                            tCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTg                                                                                        >  2:1607024/1‑150 (MQ=255)
                                                                  aTCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAgatg                                                                                  <  2:651803/150‑1 (MQ=255)
                                                                    cGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCTCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAgatgat                                                                                <  1:1249948/150‑1 (MQ=255)
                                                                    cGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAgatgat                                                                                <  2:2378012/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)
                                                                               gCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGa                                                                                                                                   >  1:621321/1‑88 (MQ=255)
                                                                                ccTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgc                                                                    <  1:1443643/150‑1 (MQ=255)
                                                                                 cTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgcg                                                                   >  2:1349087/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)
                                                                                    gtgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCAcgcgcag                                                                <  1:1892599/150‑4 (MQ=255)
                                                                                    gtgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATg                                                                <  2:602950/150‑1 (MQ=255)
                                                                                    gtgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATg                                                                >  2:1358666/1‑150 (MQ=255)
                                                                                    gtgtTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCAGTCGCGCCCTGGAGATGATGGACAAGCACGCGATg                                                                >  2:1983773/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)
                                                                                       ttGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCGGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGAc                                                             >  1:305808/1‑150 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGAtt                                                                                                                                           >  2:1816538/1‑70 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGAtt                                                                                                                                           <  1:1816538/70‑1 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCt                                                           >  1:1902083/1‑150 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCt                                                           <  1:2114694/150‑1 (MQ=255)
                                                                                         gACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCATGGAGATGATGGACAAGCACGCGATGGACCt                                                           >  2:410157/1‑150 (MQ=255)
                                                                                           ccGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGc                                                         <  1:1930787/150‑1 (MQ=255)
                                                                                             ggtggtggtGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGctc                                                       >  1:1449419/1‑150 (MQ=255)
                                                                                                 gtggtgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCa                                                   <  1:1349087/150‑1 (MQ=255)
                                                                                                   ggtgCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAATCGGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGc                                                 <  1:747139/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)
                                                                                                          tgctgcGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAgcggg                                          >  1:1771447/1‑146 (MQ=255)
                                                                                                                gTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgct                                    <  2:72557/150‑1 (MQ=255)
                                                                                                                 tGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCTCGCCCTAGGGATGGTTGGCAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgctc                                   <  1:47183/150‑1 (MQ=255)
                                                                                                                    ccTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCcgctcgct                                >  1:670322/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:775407/1‑150 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTc                            >  2:2250117/1‑150 (MQ=255)
                                                                                                                        gACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTc                            <  2:462115/150‑1 (MQ=255)
                                                                                                                          cAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCAc                          <  2:305808/150‑1 (MQ=255)
                                                                                                                           aaaCTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACg                         >  1:2098246/1‑150 (MQ=255)
                                                                                                                            aaCTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                        <  1:775407/150‑1 (MQ=255)
                                                                                                                            aaCTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                        <  1:279018/150‑1 (MQ=255)
                                                                                                                                gctgGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTggcggc                                                                                                    >  2:760591/1‑70 (MQ=255)
                                                                                                                                    gCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCgg                <  1:189593/150‑1 (MQ=255)
                                                                                                                                     cccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGt               >  2:227502/1‑150 (MQ=255)
                                                                                                                                     cccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGt               <  1:173873/150‑1 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                        <  2:1519809/140‑1 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGa                        >  1:1519809/1‑140 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGtt              >  2:104548/1‑150 (MQ=255)
                                                                                                                                      ccAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGtt              <  1:754316/150‑1 (MQ=255)
                                                                                                                                       cAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAg                                              >  2:1425050/1‑117 (MQ=255)
                                                                                                                                        aGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTa            >  2:967856/1‑150 (MQ=255)
                                                                                                                                          gAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACa          <  2:2098246/150‑1 (MQ=255)
                                                                                                                                                  ccTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAgtac  >  2:1518253/1‑150 (MQ=255)
                                                                                                                                                   cTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGa                                               <  1:1389339/104‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGGCGACCATCGTCCAGGCGGTCAAGAGCGCCCTCGAGCAATCGCCGCCGGAACTCGCCTCGGACATCGCCGAGCGCGGCCTGGTGTTGACCGGTGGTGGTGCGCTGCTGCGTGACCTCGACAAACTGCTGGCCCAGGAAACCGGCCTGCCGGTGATTGTCGCCGAGGACCCGCTGACCTGCGTCGCCCGTGGCGGCGGTCGCGCCCTGGAGATGATGGACAAGCACGCGATGGACCTGCTCTCCAGCGAGTGATCCCGCTCGCTGTTCACGAGCGCCCGGTTTACAGGTAGTAC  >  NC_002947/1076152‑1076447

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: