Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,811,485 A→T 9.2% intergenic (+50/+78) ispD → / ← PP_1615 2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,811,4850AT 9.2% 381.1 / 5.2 152intergenic (+50/+78)ispD/PP_16152‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (66/65);  new base T (2/15);  total (72/80)
Fisher's exact test for biased strand distribution p-value = 3.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

TCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCG  >  NC_002947/1811336‑1811633
                                                                                                                                                     |                                                                                                                                                    
tCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGa                                                                                                                                                      <  2:3968383/150‑1 (MQ=255)
tCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGa                                                                                                                                                      <  1:662077/150‑1 (MQ=255)
 cAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGAt                                                                                                                                                     <  2:2464975/150‑1 (MQ=255)
 cAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGAt                                                                                                                                                     <  1:1395917/150‑1 (MQ=255)
 cAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGAt                                                                                                                                                     <  1:684860/150‑1 (MQ=255)
     cgccgcGGCTGGTGGGAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCa                                                                                                                                                 <  2:5335308/150‑1 (MQ=255)
      gccgcGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCAt                                                                                                                                                <  2:1205962/150‑1 (MQ=255)
        cgcgGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCg                                                                                                                                              <  2:3518929/150‑1 (MQ=255)
            gCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCggtgg                                                                                                                                          <  1:1290856/150‑1 (MQ=255)
            gCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCggtgg                                                                                                                                          <  2:788560/150‑1 (MQ=255)
              tggtggAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGtt                                                                                                                                        <  1:1198354/150‑1 (MQ=255)
              tggtggAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGtt                                                                                                                                        <  2:3816826/150‑1 (MQ=255)
              tggtggAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGtt                                                                                                                                        >  2:5177070/1‑150 (MQ=255)
               ggtggAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTg                                                                                                                                       <  1:1719075/150‑1 (MQ=255)
                 gggAAGGGCGCAGTGACAACATCAAGGTGACCCGGCCCGCAGCTCTCGAATGGTTGCGCCAGCTCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAATCCCGCTCCTACACTGGATTGCATCGGTGGTTGGa                                                                                                                                     <  2:2036430/149‑1 (MQ=255)
                  ggAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGac                                                                                                                                    <  2:1968130/150‑1 (MQ=255)
                  ggAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGGTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGac                                                                                                                                    <  2:2707095/150‑1 (MQ=255)
                   gacgggCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTTACCCCCCGCCTGCAATGGCTTTCCTAGGGGGTTGGcca                                                                                                                                   >  1:61703/4‑150 (MQ=255)
                   gAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGCTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAAAGGATTGCATCGGTGGTTGGaca                                                                                                                                   >  2:4999020/1‑150 (MQ=255)
                      gggCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAg                                                                                                                                >  2:4161447/1‑150 (MQ=255)
                          gcAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGAcca                                                                                                                            <  1:1141770/150‑2 (MQ=255)
                          gcAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACcg                                                                                                                            <  1:3546936/150‑1 (MQ=255)
                            aGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCGCCTACAATGGATTGCATCGGTGGTTGGACACAGACcgtc                                                                                                                          >  1:4313242/1‑150 (MQ=255)
                                acaacaTAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAAAGGATTGCATCGGTGGTTTGACACAGACCGTCGTAg                                                                                                                      >  1:3513268/1‑150 (MQ=255)
                                  aacaTAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCATTTCGCGGGTAAACCCCCACCAACAAATGAATTAAAAGGTGGTTTGACAAACAACGTAGTAGGa                                                                                                                    >  1:2710908/1‑150 (MQ=255)
                                      tAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCgg                                                                                                                <  2:4393833/150‑1 (MQ=255)
                                          ggTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGttt                                                                                                            <  2:3520586/150‑1 (MQ=255)
                                          ggTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGttt                                                                                                            <  2:4742295/150‑1 (MQ=255)
                                            tGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTAc                                                                                                          <  2:3871523/150‑1 (MQ=255)
                                             gACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTAcc                                                                                                         <  1:3441199/150‑1 (MQ=255)
                                              aCCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGGGGTTGGACACAGACCGTCGTAGGAGCGGGTTTAccc                                                                                                        >  2:3004007/1‑150 (MQ=255)
                                                ccGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                      <  1:2032318/150‑1 (MQ=255)
                                                ccGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                      >  2:3775148/1‑150 (MQ=255)
                                                ccGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGGTGGACACAGACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                      >  2:5302943/1‑150 (MQ=255)
                                                ccGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCAGCGGTGGTTGGACACAGACCGGCGTAGGAGCGGGTTTACCcgc                                                                                                      >  2:5448670/1‑150 (MQ=255)
                                                 cGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCCCTCCTACAAAGGAGTGCATCGGTGGTTGGAACCAGACCGTCGTAGGAGCGGGTTTAACcgcg                                                                                                     >  2:51448/1‑150 (MQ=255)
                                                   gCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                   <  2:1911468/150‑1 (MQ=255)
                                                     ccGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGGTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTGGTAGGAGCGGTTTTACCCGCGAaga                                                                                                 <  2:3428931/150‑1 (MQ=255)
                                                     ccGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGGAGGAGCGGGTTTAACCGCGAaga                                                                                                 >  2:4607180/1‑150 (MQ=255)
                                                     ccGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAAACCGCTCCTACAAAGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                 >  1:1692083/1‑150 (MQ=255)
                                                       gAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                               >  2:2940065/1‑150 (MQ=255)
                                                       gAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                               >  2:2612879/1‑150 (MQ=255)
                                                        aaGATCTCTATTGGTTGCGCCAGCGCTGGGCGGTTAAACGGTGCGCGTTCTCCTGGGCTGGCATTTCCGCGGTTAAACCCGCTCCTAAACGGGTTTGCATGGGTGTTTGCACACAGCCCGTCTTTGGAGGGGTTTTACCCGCGAAGAGGc                                                                                              <  2:2872552/150‑1 (MQ=255)
                                                        aaGATCTCGACTGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTCCCAGGGATGGCATCGGTGGTTGGACACAGACGGTCGTAGGAGGGGTTTTACCCGCGAAGAGGc                                                                                              <  2:3971724/150‑1 (MQ=38)
                                                           aTCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAAAGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAg                                                                                           >  2:2693548/1‑150 (MQ=255)
                                                              tcGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACAAAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGGCAGTAc                                                                                        >  2:2090535/1‑150 (MQ=255)
                                                              tcGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTACTCCTGGGCTGGCCTTTTCGCGGGTAACCCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACCGCCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTAc                                                                                        <  2:2259803/150‑1 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACcgtcgt                                                                                                                        >  2:440603/1‑116 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACcgtcgt                                                                                                                        <  1:440603/116‑1 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGtt                                                                                                             >  2:1362071/1‑127 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGtt                                                                                                             <  1:1362071/127‑1 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGACGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACca                                                                                      <  2:3628846/150‑1 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGGTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCACTACca                                                                                      >  2:1160017/1‑150 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTAAAATGGTTTGCATAGGTGGTTGGACACAGACCGGCGTAGGAGCGGGTTTAACCGCGAAGGGGCCCGTACca                                                                                      >  2:4586806/1‑150 (MQ=39)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCGCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACca                                                                                      >  1:4050714/1‑150 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCGCCGACAATGGATTGCATCGGTGGTTGGACACACACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACca                                                                                      >  1:1270407/1‑150 (MQ=255)
                                                                gAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCGGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACca                                                                                      <  2:887060/150‑1 (MQ=255)
                                                                  aTGGTTGCGCCAGCGCGGGGCGGGTAAACGCTGAGCGTTCTCGTGGGCGGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGTACcagc                                                                                    <  2:5156322/150‑1 (MQ=255)
                                                                   tGGTTGCGCCCGCGCGGGGCGGGTAAACGCGGAGCGTTCTCCGGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTCCAATGGTTTGCTCGGGTGGTGGGACACAGACCGTCGTAGGAGGGGTTTTACCCGCGAAGAGGCCAGTACcagca                                                                                   <  2:2533367/150‑1 (MQ=255)
                                                                   tGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACcagcc                                                                                   >  1:4089444/1‑149 (MQ=255)
                                                                    ggTTGCGCCAGCGCTGGGCGGGTAAACGCTTAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGTTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACcagcag                                                                                  <  2:2097926/150‑1 (MQ=255)
                                                                     gTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCGACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGGAGGAGCGGGGTTACCCGCGAAGAGGCCAGTACCagccgc                                                                                 >  2:2561249/1‑150 (MQ=255)
                                                                      gtGCGCCAGCGCTGGGCGGGTAAACGCTGCGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAgcagca                                                                                <  2:1795782/149‑1 (MQ=255)
                                                                            ccgcgcTGTGCGGGTAAACGCTGCGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGTTTTGCATGGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATccc                                                                          <  1:4931292/148‑1 (MQ=38)
                                                                             aGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGa                                                                                                                               >  1:4399758/1‑96 (MQ=255)
                                                                              gcgcTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCGGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTGGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGt                                                                        <  2:4719606/150‑1 (MQ=255)
                                                                              gcgcTGGGCGGGAAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGTTTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGt                                                                        <  2:4719636/150‑1 (MQ=255)
                                                                                   gggCGGGTAACCGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGcc                                                                   <  2:3079720/150‑1 (MQ=255)
                                                                                     gCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCCCTCCTACAAGGGATTGCATCGGTGGGTGGACACAGACCGTCGTTGGGGCGGGGTTTCCCGCGAAGAGGGCAGGACCAGCAGCACATCCCGGCAGCCGa                                                                 >  1:4811824/1‑150 (MQ=255)
                                                                                      caggagaaCGCTCAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGCTTGCATCGGGGGTTGGACCCAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACCTCCCGGCCGCCGaa                                                                >  2:3876032/7‑150 (MQ=37)
                                                                                         gTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGTTTGCATCGGTGGTTTGACACATAGCGTGGTTGGTGCGGGTTTTACCGCGTAGTTGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACg                                                             >  2:728711/1‑150 (MQ=255)
                                                                                                 tGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACCATGGATTGCATCGGTGGTTGGACACAGACCGGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCCGCACCTCCCGTCAGCCGAAACGAGACTCCg                                                     >  2:4812436/1‑150 (MQ=38)
                                                                                                 tGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGTGCGGGTTTACCCGCGCAGAGGCCAGTACCAGCAGCACATCCCGGCAGCCGCAACGATACTCCg                                                     >  1:4893317/1‑150 (MQ=255)
                                                                                                 ggAGCTTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCGGCTCCTACAATGGTTTGCATCGGTGGTTGGACACAGACCGTCGTCGGAGGGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCg                                                     <  1:4017546/149‑1 (MQ=39)
                                                                                                        tctcCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACCAGGTATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAAccc                                              <  1:5302943/150‑1 (MQ=255)
                                                                                                            gtGGGCTGGCCTTTTCGCGGGTAAACCCGCGCCTCCATTGGTTTGCATGGGTGGTTGGACACAGCCCGTCGTAGGCGGGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGcc                                          <  1:2653037/149‑1 (MQ=38)
                                                                                                            cTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGTTTGCATGGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGcc                                          <  1:452792/150‑1 (MQ=255)
                                                                                                               ggctggctTTTTCGGGGGTAAACCCGCTCCTACATTGTTTTGCATGGGTGGTTGCACACTGACGGTGTTTGGTGGGGTTTAACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAAc                                       <  1:3775148/150‑1 (MQ=255)
                                                                                                                   ggCCTTTGCGCGGGAACCCCCGCTCATCCAATGGATGGCTGGGGTGTTGGCACACAGCCCGGCGTCGGCGGGGGTTAACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCcct                                   <  1:2968871/150‑1 (MQ=255)
                                                                                                                     ccTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGtt                                                                                                                                        >  1:4449348/1‑47 (MQ=255)
                                                                                                                     ccTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGtt                                                                                                                                        <  2:4449348/47‑1 (MQ=255)
                                                                                                                     cattttCGGGGGAAAACCCTCTACTACTTTGGTTTGCTTGGGTTTTTGTACACAGACCGTCGTAGGAGCGTTTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCcctcc                                 <  1:3120930/148‑1 (MQ=255)
                                                                                                                       ttttCGCGGGTAAACCCGCTCCTACAATGGATTTCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTAACCGCGAAGAGGCCAGTACCAGCAGCACAACCCGTCAGCCCAAACGATACACCGGCAACCCCCCCAACCccaccgt                               >  2:5111244/1‑148 (MQ=38)
                                                                                                                         ttCGCGGGTAAACCCGCTCCTACAATGTATTGCATGGGTGGTTGGACACAGACCGTCGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCa                             <  2:1100142/150‑1 (MQ=255)
                                                                                                                          tCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAg                            <  2:4030447/150‑1 (MQ=255)
                                                                                                                           cgcgGGTAAACCCGCTCATACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAgg                           <  2:61703/150‑1 (MQ=255)
                                                                                                                            gcgGGTTTACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTGCAGGt                          >  2:1705588/1‑150 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTc                                                                       >  1:4651303/1‑105 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTc                                                                       <  2:4651303/105‑1 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGt                          >  2:4589487/1‑150 (MQ=255)
                                                                                                                                 tAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATcc                     >  1:3920962/1‑150 (MQ=255)
                                                                                                                                 tAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATcc                     >  2:4332450/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATccacca                 >  2:1390985/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATccacca                 >  1:3059858/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                              <  1:343057/72‑1 (MQ=255)
                                                                                                                                             tACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGt         >  2:386883/1‑150 (MQ=255)
                                                                                                                                             tACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGt         >  2:1838223/1‑150 (MQ=255)
                                                                                                                                               cAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATGCTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTAc       >  1:312043/1‑150 (MQ=255)
                                                                                                                                               cAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTAc       >  1:4983922/1‑150 (MQ=255)
                                                                                                                                                    gATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCg  <  2:4050714/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGGCGCCGCGGCTGGTGGAAGGGCGCAGTGACAACATAAAGGTGACCCGGCCCGAAGATCTCGAATGGTTGCGCCAGCGCTGGGCGGGTAAACGCTGAGCGTTCTCCTGGGCTGGCCTTTTCGCGGGTAAACCCGCTCCTACAATGGATTGCATCGGTGGTTGGACACAGACCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCG  >  NC_002947/1811336‑1811633

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