Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,809,640 T→A 5.3% intergenic (+54/‑118) accA → / → tilS acetyl‑CoA carboxylase carboxyltransferase subunit alpha/tRNA lysidine(34) synthetase TilS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,809,6400TA 5.3% 945.8 / 28.5 281intergenic (+54/‑118)accA/tilSacetyl‑CoA carboxylase carboxyltransferase subunit alpha/tRNA lysidine(34) synthetase TilS
Reads supporting (aligned to +/- strand):  ref base T (119/147);  new base A (6/9);  total (125/156)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 283 total aligned reads displayed.

GCCGATCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTTGGCTCAACGCCCC  >  NC_002947/1809491‑1809789
                                                                                                                                                     |                                                                                                                                                     
gCCGATCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTGGTGAGTGATGGCACCGGCtt                                                                                                                                                       >  2:2922451/1‑150 (MQ=255)
   gATCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGc                                                                                                                                                    <  2:1901854/150‑1 (MQ=255)
    aTCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGcc                                                                                                                                                   >  2:926313/1‑150 (MQ=255)
    aTCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGcc                                                                                                                                                   <  1:3530892/150‑1 (MQ=255)
     tCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCg                                                                                                                                                  <  1:4140024/150‑1 (MQ=255)
     tCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCg                                                                                                                                                  <  1:2000688/150‑1 (MQ=255)
      cTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGTTGGCACCGGCTTCGCCgg                                                                                                                                                 >  2:5108955/1‑150 (MQ=255)
         gTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCGAAGCCGgtgc                                                                                                                                              <  1:4668352/150‑2 (MQ=255)
          tGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCGAAGCCGgtgcc                                                                                                                                             <  1:1012106/150‑3 (MQ=255)
          tGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCGAAGCCGgtgcc                                                                                                                                             <  1:296090/150‑3 (MQ=255)
            cagcagCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCGAAGCCGgtgccat                                                                                                                                           <  2:3015880/150‑5 (MQ=255)
             agcagcTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTcgc                                                                                                                                          <  1:5540850/150‑1 (MQ=255)
              gcagcTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                                                                                         <  2:2793960/150‑1 (MQ=255)
              gcagcTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTGGTGATTGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                                                                                         <  2:3375038/150‑1 (MQ=255)
                agcTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCggg                                                                                                                                       >  2:238878/1‑150 (MQ=255)
                     gACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAc                                                                                                                                  <  2:2109158/150‑1 (MQ=255)
                          gCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCt                                                                                                                             >  2:5143794/1‑150 (MQ=255)
                               gCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTTTACCCGcgaacac                                                                                                                        <  2:5329849/150‑7 (MQ=255)
                               gCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTAc                                                                                                                        <  1:2731019/150‑1 (MQ=255)
                                cAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACa                                                                                                                       <  1:5431022/150‑1 (MQ=255)
                                 aaGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAg                                                                                                                      <  1:1463881/150‑1 (MQ=255)
                                 aaGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAg                                                                                                                      <  1:2544467/150‑1 (MQ=255)
                                 aaGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAg                                                                                                                      >  1:3786986/1‑150 (MQ=255)
                                  aGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                     <  1:447315/150‑1 (MQ=255)
                                  aGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                     <  1:2185936/150‑1 (MQ=255)
                                    ttCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTa                                                                                                                   >  2:3414353/1‑150 (MQ=255)
                                     tCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTAt                                                                                                                  >  1:3287655/1‑150 (MQ=255)
                                      cGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTAtt                                                                                                                 <  2:283639/150‑1 (MQ=255)
                                       gACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAAGTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTGGGGATTGATGGCACCGGCTTCGCCGGTGGTCGCGGGTAAACCCGCTCCTACAGGTCTTg                                                                                                                >  1:1257020/1‑150 (MQ=38)
                                        aCAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTa                                                                                                                   >  2:4778454/1‑146 (MQ=255)
                                          aaCGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAAGTCTGATTTGCCCCGGGGCCGGTTTGCGGCCCGTGTGATTGATGGCACCGGCTTCGCCGGGGTTCGCGGGTAAACCCGCTCCTACAGGTAGTgcgc                                                                                                             >  1:4070771/1‑150 (MQ=38)
                                            cGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTg                                                                                                           <  1:3082362/150‑1 (MQ=255)
                                            cGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTg                                                                                                           >  2:1394519/1‑150 (MQ=255)
                                               tGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGccc                                                                                                        <  2:3470315/150‑1 (MQ=255)
                                                gcgcTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCt                                                                                                       >  2:742154/1‑150 (MQ=255)
                                                gcgcTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCt                                                                                                       <  1:119205/150‑1 (MQ=255)
                                                gcgcTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCt                                                                                                       <  2:3891406/150‑1 (MQ=255)
                                                 cgcTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                                                                                        <  1:2853991/116‑1 (MQ=255)
                                                 cgcTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                                                                                        >  2:2853991/1‑116 (MQ=255)
                                                   cTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGa                                                                                                    <  2:2062214/150‑1 (MQ=255)
                                                   cTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGa                                                                                                    <  1:3450776/150‑1 (MQ=255)
                                                   cTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGa                                                                                                    >  1:3355731/1‑150 (MQ=255)
                                                    tACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAg                                                                                                   >  2:5216284/1‑150 (MQ=255)
                                                       tGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCGTTGAGGCt                                                                                                >  2:1340391/1‑150 (MQ=255)
                                                           ccGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTgcg                                                                                            >  2:1619243/1‑150 (MQ=255)
                                                             gtcgtTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGcc                                                                                          >  1:292850/1‑150 (MQ=255)
                                                             gtcgtTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGcc                                                                                          <  2:271742/150‑1 (MQ=255)
                                                                       gcgcCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTgg                                                                                <  2:1241693/150‑1 (MQ=255)
                                                                       gcgcCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTgg                                                                                <  2:5544960/150‑1 (MQ=255)
                                                                       gcgcCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCGGCCCTTGAGGCTTGCGCCGTACGTGTgg                                                                                <  1:2882874/150‑1 (MQ=255)
                                                                         gcCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGa                                                                              <  1:655302/150‑1 (MQ=255)
                                                                         gcCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGa                                                                              >  1:4368537/1‑150 (MQ=255)
                                                                            tgatgaGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCg                                                                           >  1:4551685/1‑150 (MQ=255)
                                                                             gatgaGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGGTTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCGGCTCTTGAGGCTTGCGCCGTACGTGTGGGAGCgg                                                                          <  2:1938166/150‑1 (MQ=255)
                                                                             gatgaGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCgg                                                                          <  1:1708059/150‑1 (MQ=255)
                                                                               tgaGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGc                                                                        >  2:2649602/1‑150 (MQ=255)
                                                                                gaGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCa                                                                       <  2:2534342/150‑1 (MQ=255)
                                                                                 aGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCAt                                                                      <  2:3266776/150‑1 (MQ=255)
                                                                                  gCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATg                                                                     >  1:224820/1‑150 (MQ=255)
                                                                                  gCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCAAGCCTCAAGGGCAGCGCAATACCTGTAGGAGCgggttta                                                                     <  1:1380736/150‑5 (MQ=25)
                                                                                      cGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGccc                                                                  >  1:4772616/1‑149 (MQ=255)
                                                                                       ggCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCcgc                                                                <  1:2396393/150‑1 (MQ=255)
                                                                                         ccTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTgg                                                                                >  2:2227343/1‑132 (MQ=255)
                                                                                          ctctAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGaa                                                             >  1:5298874/1‑150 (MQ=255)
                                                                                           tctAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAGCCCGCTCCCACACGTACGGCGCAAGCCTCAAGGGCAGCGCAATACCTGTAGGAGCgggtttaccc                                                                  >  1:1911252/1‑137 (MQ=9)
                                                                                           tctAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAAt                                                            <  2:505227/150‑1 (MQ=255)
                                                                                           tctAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAAt                                                            <  1:787928/150‑1 (MQ=255)
                                                                                             tAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGc                                                          >  1:368607/1‑150 (MQ=255)
                                                                                               aTTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCtt                                                        >  1:4813561/1‑150 (MQ=255)
                                                                                               aTTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCtt                                                        <  1:3905834/150‑1 (MQ=255)
                                                                                                 tCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTc                                                      >  2:4596971/1‑150 (MQ=255)
                                                                                                  cTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCt                                                     <  2:3990624/150‑1 (MQ=255)
                                                                                                  cTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCt                                                     <  2:5276869/150‑1 (MQ=255)
                                                                                                     aTTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAAt                                                            <  1:2668592/140‑1 (MQ=255)
                                                                                                      tttGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCTAATGCTTTCTGGaa                                                 <  2:1478634/150‑1 (MQ=255)
                                                                                                      tttGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGaa                                                 <  2:1109428/150‑1 (MQ=255)
                                                                                                         gCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAAAcg                                              >  1:2432351/1‑150 (MQ=255)
                                                                                                            ccGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATcgagc                                           >  2:869329/1‑150 (MQ=255)
                                                                                                            ccGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATcgagc                                           <  1:3559803/150‑1 (MQ=255)
                                                                                                            ccGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATcgagc                                           >  1:238041/1‑150 (MQ=255)
                                                                                                               gggCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATcgagcgag                                        <  1:2155358/150‑1 (MQ=255)
                                                                                                                 gCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAggg                                      <  2:4220213/150‑1 (MQ=255)
                                                                                                                  ccGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGc                                     <  1:422719/150‑1 (MQ=255)
                                                                                                                     ttttGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAg                                  <  1:3238952/150‑1 (MQ=255)
                                                                                                                     ttttGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAg                                  <  2:2499847/150‑1 (MQ=255)
                                                                                                                      tttGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAga                                 <  2:4070810/150‑1 (MQ=255)
                                                                                                                       ttGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAgat                                <  2:5413168/150‑1 (MQ=255)
                                                                                                                         gCGGCCCTTGTGATTGATGGCACCGGCTTCTCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGATGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAgatga                              <  1:936246/150‑1 (MQ=255)
                                                                                                                             cccTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCaa                          <  2:768314/150‑1 (MQ=255)
                                                                                                                              ccTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGTGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGTCTAGATGATCAAc                         <  2:4697433/150‑1 (MQ=255)
                                                                                                                              ccTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAAc                         <  1:5194165/150‑1 (MQ=255)
                                                                                                                                ttGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCt                       <  1:3488842/150‑1 (MQ=255)
                                                                                                                                   tgattgatGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCAc                    >  2:2691055/1‑150 (MQ=255)
                                                                                                                                     attgatGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTc                  >  2:431368/1‑150 (MQ=255)
                                                                                                                                      ccgcgaacACCGGCGAAGCCGGTGCTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTcc                 >  1:4233405/8‑150 (MQ=39)
                                                                                                                                          tGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCGGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTTgg             <  2:3725056/150‑1 (MQ=255)
                                                                                                                                          gaacACCGGCGAAGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTTgg             >  2:4347463/4‑150 (MQ=255)
                                                                                                                                            acACCGGCGAAGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAAAGCATTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTAGGCt           >  2:2332809/2‑150 (MQ=39)
                                                                                                                                             cACCGGCGAAGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTTGGCTc          >  2:2006507/1‑150 (MQ=255)
                                                                                                                                              aCCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTTGGCTCa         >  1:4076358/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGATCTGGTGCAGCAGCTGGACATGCTCGGCAAGTTCGACAACGATGCGCTACTGGCCCGTCGTTACGAGCGCCTGATGAGCTACGGCCTCTAATTCTGATTTGCCCCGGGGCCGTTTTGCGGCCCTTGTGATTGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGGTATTGCGCTGCCCTTGAGGCTTGCGCCGTACGTGTGGGAGCGGGCATGCCCGCGAATGCTTTCTGGAATCGAGCGAGGGCTAGATGATCAACCTCACTCCTTGGCTCAACGCCCC  >  NC_002947/1809491‑1809789

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