Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,562,785 A→C 11.0% intergenic (‑156/+175) PP_4888 ← / ← purA methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,562,7850AC11.0% 675.2 / 12.9 241intergenic (‑156/+175)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (118/89);  new base C (3/31);  total (121/120)
Fisher's exact test for biased strand distribution p-value = 5.80e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

CCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCG  >  NC_002947/5562638‑5562934
                                                                                                                                                   |                                                                                                                                                     
ccTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACg                                                                                                                                                     >  2:1547284/1‑150 (MQ=255)
  tCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGc                                                                                                                                                   >  2:317061/1‑150 (MQ=255)
     tgtTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTg                                                                                                                                                >  1:2692181/1‑150 (MQ=255)
     tgtTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGACAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTg                                                                                                                                                >  2:2377876/1‑150 (MQ=255)
       ttttCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTgcg                                                                                                                                              <  1:2343229/150‑1 (MQ=255)
          tCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGACATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGccc                                                                                                                                           <  1:3844250/150‑1 (MQ=255)
            tGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCt                                                                                                                                                 >  1:1507026/1‑142 (MQ=255)
             ggCGCCACTATTTTGATGGCATCTACTGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCTGGCGCATTCGCGGGTTAACCCGCTCCTACAACGGCCTGCGCCCAAc                                                                                                                                        >  1:622025/1‑150 (MQ=255)
                gcCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCAt                                                                                                                                     >  1:4935955/1‑150 (MQ=255)
                gcCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCAt                                                                                                                                     >  2:2817447/1‑150 (MQ=255)
                      aTTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgt                                                                                                                               <  2:1842962/150‑1 (MQ=255)
                      aTTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgt                                                                                                                               <  2:4675710/150‑1 (MQ=255)
                      aTTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgt                                                                                                                               <  1:3892092/150‑1 (MQ=255)
                        tttGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgtgt                                                                                                                             <  1:4834978/150‑1 (MQ=255)
                        tttGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgtgt                                                                                                                             <  1:2447743/150‑1 (MQ=255)
                        tttGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgtgt                                                                                                                             <  2:5263702/150‑1 (MQ=255)
                            aTGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGccc                                                                                                                         <  2:3241993/150‑1 (MQ=255)
                             tGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAAACCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCca                                                                                                                        >  1:3321150/1‑150 (MQ=255)
                                 aTCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATa                                                                                                                    <  2:1413753/150‑1 (MQ=255)
                                  tCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAg                                                                                                                   >  2:3567519/1‑150 (MQ=255)
                                      cAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGc                                                                                                               >  2:4157850/1‑150 (MQ=255)
                                       aGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCg                                                                                                              >  1:5602622/1‑150 (MQ=255)
                                        gAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCgg                                                                                                             <  2:3539538/150‑1 (MQ=255)
                                        gAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCgg                                                                                                             <  1:4157850/150‑1 (MQ=255)
                                         aaaCTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCggg                                                                                                            <  1:2196119/150‑1 (MQ=255)
                                          atcttcttACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGACTAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTAGCGGGTAAACCCGCTCCTACAACGGCATGAGCACAACCATCGTTGTGTGCCCTCATAGGAGCGGGt                                                                                                           <  2:2542291/148‑1 (MQ=255)
                                             ttcttACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTa                                                                                                        <  2:2505811/150‑1 (MQ=255)
                                               cttACCGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATTGGGTTTCGTGTGCAGGCCCATTCGCGTGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTAcc                                                                                                      <  2:224822/150‑1 (MQ=255)
                                               cttACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTAcc                                                                                                      >  1:5163154/1‑150 (MQ=255)
                                                ttACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAGACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCgg                                                                                                             >  2:3081930/1‑142 (MQ=255)
                                                   cAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCcgcg                                                                                                  <  2:3638755/150‑1 (MQ=255)
                                                    aGCAAGGGTTGTGCCGCAGGCATCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAATACCGCTCCTAAAACGGCAAGCGCACACCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGa                                                                                                 >  2:2581771/1‑150 (MQ=255)
                                                     gCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGaa                                                                                                >  2:1454391/1‑150 (MQ=255)
                                                      cAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGact                                                                                               >  1:5206942/1‑148 (MQ=255)
                                                      cAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAAt                                                                                               >  1:3051231/1‑150 (MQ=255)
                                                       aaGGGTTGTGCCGCAGGCAGCGGGGCAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCGGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATa                                                                                              <  2:1945568/150‑1 (MQ=255)
                                                       aaGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATa                                                                                              <  2:4384971/150‑1 (MQ=255)
                                                       aaGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATa                                                                                              >  1:3167546/1‑150 (MQ=255)
                                                       aaGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACAGAGGAGCGGGTTTACCCGCGAATa                                                                                              >  2:4901832/1‑150 (MQ=255)
                                                       aaGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCCCACCTACAACGGCCTGCGCCCAACCAACGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATa                                                                                              >  1:168891/1‑150 (MQ=255)
                                                            ttGTGCCGCTGGCAGCGGGGAACTGCGCCACGGTCTAGCATTGGGGTTTGGTGTGCAGGCCCTTTCGCGGGTAACCCCGCTCCTACACCGGCCTGCGCCCAACCATCGTTGGGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTa                                                                                         <  2:5184659/150‑1 (MQ=38)
                                                                gCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGCGCGGGTTTACCCGCGAATACGGTAgcgg                                                                                     >  1:1250265/1‑150 (MQ=255)
                                                                 ccGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAgcggc                                                                                    >  2:890642/1‑150 (MQ=255)
                                                                 ccGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGGTGTGTGCCCACCTAGGAGCGGGGTTACCCGCGCATACGGTAgcggc                                                                                    >  1:3583915/1‑150 (MQ=255)
                                                                           cGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGTTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCAACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTg                                                                          <  2:4052741/150‑1 (MQ=255)
                                                                           cGGGGAAATGCGACATGGTCTAGCATTGGGGTTTCGGGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACACCGGCCTGCGCCCACCCATCGTTTTGTGCCAACATAGGAGCGGGTTTCCCCGCGAATACGGTAGCGGCGGCAACGGTg                                                                          <  2:964419/150‑1 (MQ=38)
                                                                            ggggAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTAACCGCGAATACGGTAGCGGCGGCAACGGTgg                                                                         >  1:3536119/1‑150 (MQ=255)
                                                                             gtgAAATGCGCCAAGGTCTAGCAATGGGGTTTCTTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGa                                                                        <  2:5105401/148‑1 (MQ=255)
                                                                             ggtaaaTGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGa                                                                        <  2:4442074/147‑1 (MQ=255)
                                                                             gggAAATGCGCCCCGGTCTAGCATTGGGTTTTCGTGTGCTGCCCCATTCGCGGCTAAACCCGCTCCTACACCGGCCTGCGCCCAACCATGGTTGTGTGCCCAAATCGGCGGGGGTTAACCCGCGAATACGGTAGCGGCGGCAACGGTGGa                                                                        <  1:135436/150‑1 (MQ=37)
                                                                             gggAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCCCCTAGGAGCGGGTTTACCCGCGAATACGGGAGCGGCGGCACCGGTGGa                                                                        >  2:3039537/1‑150 (MQ=255)
                                                                              ggAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCGACCCCGGCCTGCGCCCCACCCTCGGTGTGTGCCCACAGAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTggct                                                                       >  1:3868355/1‑148 (MQ=38)
                                                                                    gcgcCACGGTCTAGCATTGGGGTTTGGTGTGCCGGCCCATTGGCGGGTAACCCCGCTCCTCCACCGGCCGGCGCCCAACCCTGGTTGTGGGCCCACATCGGCGGGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCgg                                                                 <  1:805942/150‑1 (MQ=25)
                                                                                     cgcCACGGTCTAGCATTGGGGTTTGGTGTGCAGGCCCATGGGGGGGTAACCCCGCTCCAACACCGGCCGGCGCCCAACCATGGTTGTGGGCCCACATCGGCGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGt                                                                <  1:1098681/150‑1 (MQ=25)
                                                                                       ccAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGCATACGGTAGCGGCGGCAACGGTGGATTACCGGTgg                                                              >  2:218672/1‑150 (MQ=255)
                                                                                       ccAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTgg                                                              >  2:5044839/1‑150 (MQ=255)
                                                                                        cAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCGACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCCCAGAGGAGCGGGTTTAACCGCGACTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGa                                                             >  1:2920450/1‑150 (MQ=38)
                                                                                            gTCTAGCATTGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATAtt                                                         <  1:1855105/150‑1 (MQ=255)
                                                                                             tCTAGCATTGGGGTTTGGTGTGCAGGCCCATTCGGGGGTAACCCCGCTCATACACCGGCCGGCGCCAACCCTGCGTTGTGTCCCCACATCGGCGCGGTTTAACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTg                                                        <  1:4792347/150‑1 (MQ=25)
                                                                                             tCTAGCATTGGGGTTTCGTGTGCCGGCCCATTCGCGGGTAACCCCGCCCCTACACCGGCCGGCGCCCAACCATCGTTGTGTGCCCACACCGGCGGGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTg                                                        <  1:1547284/150‑1 (MQ=38)
                                                                                             tCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTg                                                        <  2:4856763/150‑1 (MQ=255)
                                                                                               tAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGTGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGc                                                      >  2:1262675/1‑150 (MQ=255)
                                                                                                  cAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACCACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAg                                                   <  2:1943483/150‑1 (MQ=255)
                                                                                                    aTGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGca                                                 >  1:428357/1‑150 (MQ=255)
                                                                                                     tGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGACTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGCTATTGGCCAGccc                                                >  2:1079755/1‑148 (MQ=255)
                                                                                                     tGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGACTACGGTAGCGGCGGCAACGGGGGATTACCGGTGGATATTGGCCAGcac                                                >  2:3412145/1‑150 (MQ=255)
                                                                                                      ggggTTTCGTGTGCAGGCCCTTTCGCGGGTAAACCCGCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGcaca                                               <  2:4521221/150‑1 (MQ=255)
                                                                                                      ggggTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGcacc                                               >  1:4522863/1‑149 (MQ=255)
                                                                                                        ggTTTCGTGTGCTGGCCCATTCGCGGTTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGa                                             <  2:1181457/150‑1 (MQ=255)
                                                                                                        ggTTTCGTGTGCTGGCCAATTAGCGTGTAAACCCGCTCATACACCGCCCTGCGCCCAACCATGGTTTTTTCCACACATAGTAGCGTTTTTACCCGCTAATACGGTAGCTGCGGCATCTGTGTATTACCTGTGGATATTGGCCAGCACAGa                                             <  2:766377/150‑1 (MQ=255)
                                                                                                        ggTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCGACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATAGTGGCCAGCACAGa                                             >  1:3986988/1‑150 (MQ=255)
                                                                                                         gTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCGACAACGGCCGGCGCCCAACCATCGGTGTGTGCCCACAAAGGAGCGGGTTTACCCCCGAATACGGTAGCGGCGGCAACGGTGGGTTACCGGTGGGTTTTGGGCCGCACAGAc                                            >  1:3768570/1‑150 (MQ=37)
                                                                                                             cGTGTGCAGGCCCATTCGCGGTTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCTTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGc                                        <  2:1039973/150‑1 (MQ=255)
                                                                                                                 tgCAGGCCCATTCGCGGGTAAACCCGCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACAGAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCcac                                    <  2:1190158/150‑1 (MQ=255)
                                                                                                                    aGGCCCTTTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACaaa                                 <  2:5343927/150‑1 (MQ=255)
                                                                                                                    aGGCCCATTCGCGGTTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACaaa                                 <  2:2950551/150‑1 (MQ=255)
                                                                                                                    aGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATCGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACaaa                                 <  2:2845122/150‑1 (MQ=255)
                                                                                                                     ggCCCTTTGGGGGGTAACCCCGCTCCTACAAGGGCCGGCGCCCAACCTTGGTTGTGTGCCCACTTTGGAGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGTTTACCGGTGGATATTGGCCTGCACAGACATGCCCACAAAt                                <  1:4645239/150‑1 (MQ=25)
                                                                                                                      gCCCATTCGCGGGTAAAACCGCGCCGACAACGGCCTGCGCCCAACCCACGTTGTGTGCCCCCATAGGAGCGGGTTTACCCGCGAAAACGGGAGCGGCGGCAACCGTGGATTAACGGTGGATATTGGCCAGCACACACATGCCCACAAATg                               >  1:3377082/1‑150 (MQ=37)
                                                                                                                        ccATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTCCCGGTGGATATTGGCCAGCACAGACATGCCCACCAAGGaa                             >  2:5628107/1‑150 (MQ=255)
                                                                                                                         cATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCCTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGGAGCGGCGGCAACGGGGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGaaa                            >  1:1162047/1‑150 (MQ=255)
                                                                                                                          aTTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGaaaa                           >  2:1903487/1‑150 (MQ=255)
                                                                                                                           ttCGCGGGTAAACCCGCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAAc                          <  2:1671426/150‑1 (MQ=255)
                                                                                                                           ttCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAAc                          <  1:405676/150‑1 (MQ=255)
                                                                                                                             cgcgGGTAAAACCGCTCCAACAACGGCCAGCGCCCAACCATCGTTTTTTGCCAACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGAACAGACATGCCCACAAATGAAAAcgc                        >  2:1323933/1‑150 (MQ=37)
                                                                                                                              gcgGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAAcgcc                       >  2:1753466/1‑150 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAAcgcc                       >  1:1345067/1‑150 (MQ=255)
                                                                                                                                    aaccccGCTCCTACACGGGCGTGCGCCAACCCTTGTTTGTGTGCCAACTTTGTTGGGGTTTTACCCGCTAATACGGTAGCGGCGGCAACGGTGGATTACGGGTGTATATTGGCCAGCACATACATGCCCACAAATGAAAACGCCGCGGcc                 <  1:2874428/147‑1 (MQ=19)
                                                                                                                                     aaCCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGc                                                      <  1:370034/112‑1 (MQ=255)
                                                                                                                                      aCCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCt               >  2:3714902/1‑150 (MQ=255)
                                                                                                                                        ccGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCtgc             >  2:3800606/1‑148 (MQ=255)
                                                                                                                                          gCTCCTACCACGGCCTGCGCCCAACCATGGTTGTGGGCCCACATAGGAGCGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtctc           <  2:3771166/150‑1 (MQ=255)
                                                                                                                                          gCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTGACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtctc           <  2:5512354/150‑1 (MQ=255)
                                                                                                                                          gCTCCTACAACGGCCTGCGCCCAACAATCGTTGTGTGCCCACATAGGAGGGGTTTTACCCGCGAATACGGTAGCGGCGGAAACGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtctc           <  1:2050077/150‑1 (MQ=39)
                                                                                                                                            tCCTACAACGGCCTGCGCCAAACAATCTTTGTGTGCCCAAATAGGAGCGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAg         <  1:438072/150‑1 (MQ=38)
                                                                                                                                             ccTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGCTCAgg        >  2:4203777/1‑150 (MQ=255)
                                                                                                                                               tACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGc      >  1:4685532/1‑150 (MQ=255)
                                                                                                                                               tACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGc      >  2:2640487/1‑150 (MQ=255)
                                                                                                                                                  aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgc   <  2:4935955/150‑1 (MQ=255)
                                                                                                                                                   aCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgcg  >  1:5614018/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCG  >  NC_002947/5562638‑5562934

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