Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 479,496 A→G 5.4% intergenic (‑25/+31) cca ← / ← ycgB multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947479,4960AG 5.4% 1253.4 / 35.7 386intergenic (‑25/+31)cca/ycgBmultifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase/type IV piliation protein
Reads supporting (aligned to +/- strand):  ref base A (167/198);  new base G (10/11);  total (177/209)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 400 total aligned reads displayed.

ACCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGT  >  NC_002947/479347‑479645
                                                                                                                                                     |                                                                                                                                                     
aCCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACa                                                                                                                                                       <  1:7352364/150‑1 (MQ=255)
aCCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACa                                                                                                                                                       <  2:2716887/150‑1 (MQ=255)
 ccGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAg                                                                                                                                                      >  2:164445/1‑150 (MQ=255)
  cGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAgg                                                                                                                                                     <  2:5337131/150‑1 (MQ=255)
  cGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAgg                                                                                                                                                     <  1:4572012/150‑1 (MQ=255)
       gATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAAc                                                                                                                                                <  1:6852967/150‑1 (MQ=255)
       gATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACCGCCGACGGGTTatc                                                                                                                                                <  1:1803977/150‑3 (MQ=255)
          aGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCg                                                                                                                                             <  1:5427658/150‑1 (MQ=255)
            ccATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTc                                                                                                                                           >  1:2813381/1‑150 (MQ=255)
             cATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTcg                                                                                                                                          >  1:747936/1‑150 (MQ=255)
             cATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTcg                                                                                                                                          >  2:7750490/1‑150 (MQ=255)
              aTTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTcgg                                                                                                                                         >  1:2405626/1‑150 (MQ=255)
               tttGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTcggc                                                                                                                                        >  1:5283922/1‑150 (MQ=255)
                ttGCTAGCGTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACGGGTTATCCTGTcggca                                                                                                                                       <  2:3554279/150‑2 (MQ=255)
                   cTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTc                                                                                                                                    >  1:6178770/1‑150 (MQ=255)
                    tAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCa                                                                                                                                   >  2:7567836/1‑150 (MQ=255)
                        aTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATGACCGCCGACGGGTTATCCTGTCGGCAGTTATgga                                                                                                                               <  2:3068624/150‑2 (MQ=38)
                        aTTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATGACCGCCGACGGGTTATCCTGTCGGCAGTTATgga                                                                                                                               <  2:7406160/150‑2 (MQ=38)
                         tttCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATgggg                                                                                                                              >  2:5716620/1‑150 (MQ=255)
                           tcttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGGCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGtt                                                                                                                            <  2:6985604/150‑1 (MQ=255)
                             ttcGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGa                                                                                                                                                    <  2:7422831/124‑1 (MQ=255)
                               cGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCAt                                                                                                                        <  1:6681720/150‑1 (MQ=255)
                                gACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATc                                                                                                                       <  1:2246478/150‑1 (MQ=255)
                                   ggTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                                    <  2:8281706/150‑1 (MQ=255)
                                   ggTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                                    <  1:7175690/150‑1 (MQ=255)
                                   ggTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCaga                                                                                                                    <  1:7838299/150‑1 (MQ=255)
                                    gTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCagag                                                                                                                   <  1:6207011/150‑1 (MQ=255)
                                     tGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAgg                                                                                                                  <  2:427315/150‑1 (MQ=255)
                                         gcCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                              <  1:447735/150‑1 (MQ=255)
                                         gcCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAggtggt                                                                                                              <  1:521934/150‑1 (MQ=255)
                                           ccaccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGtggtgg                                                                                                            <  2:2409622/150‑1 (MQ=255)
                                           ccaccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGtggtgg                                                                                                            >  1:5096989/1‑150 (MQ=255)
                                           ccaccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGtggtgg                                                                                                            >  2:5923514/1‑150 (MQ=255)
                                            caccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCCTCAGAGGTGGTGGa                                                                                                           >  1:2616239/1‑150 (MQ=255)
                                            caccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGa                                                                                                           >  2:7987707/1‑150 (MQ=255)
                                            caccaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGa                                                                                                           <  2:1029296/150‑1 (MQ=255)
                                             accaccaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGAt                                                                                                          >  2:5528033/1‑150 (MQ=255)
                                                accaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAAc                                                                                                       >  2:1288101/1‑150 (MQ=255)
                                                accaGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAAc                                                                                                       >  2:3328314/1‑150 (MQ=255)
                                                 ccagccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                                      <  2:2064857/150‑1 (MQ=255)
                                                 ccagccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                                      <  1:1630182/150‑1 (MQ=255)
                                                 ccagccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACa                                                                                                      <  1:1393668/150‑1 (MQ=255)
                                                   agccagTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGc                                                                                                    <  1:2015157/150‑1 (MQ=255)
                                                        gTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGt                                                                                               >  1:2069525/1‑150 (MQ=255)
                                                         tCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTc                                                                                              >  2:4833013/1‑150 (MQ=255)
                                                            atatCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCAt                                                                                           <  2:6579753/150‑1 (MQ=255)
                                                              atCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGc                                                                                         <  1:6487376/150‑1 (MQ=255)
                                                               tCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCg                                                                                        >  2:8142407/1‑150 (MQ=255)
                                                                cGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCgg                                                                                       <  1:3974079/150‑1 (MQ=255)
                                                                 gCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGc                                                                                      <  1:3676222/150‑1 (MQ=255)
                                                                   tGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCa                                                                                    <  2:5154529/150‑1 (MQ=255)
                                                                   tGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCa                                                                                    >  1:3125848/1‑150 (MQ=255)
                                                                     aCAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATa                                                                                  <  1:5770258/150‑1 (MQ=255)
                                                                       aGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGt                                                                                >  1:4732543/1‑150 (MQ=255)
                                                                           cggcCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGgc                                                                            <  2:4011682/150‑1 (MQ=255)
                                                                             gcCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGgcgc                                                                          <  1:8597676/150‑1 (MQ=255)
                                                                               cGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCtt                                                                        <  1:8246352/150‑1 (MQ=255)
                                                                                gAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCttt                                                                       <  1:3868352/150‑1 (MQ=255)
                                                                                    aaaCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCt                                                                   >  1:5611171/1‑150 (MQ=255)
                                                                                      aCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTAt                                                                 <  2:7889581/150‑1 (MQ=255)
                                                                                      aCGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTAt                                                                 <  1:2750119/150‑1 (MQ=255)
                                                                                          tCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTc                                                             <  1:4409596/150‑1 (MQ=255)
                                                                                           cgcgCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTcg                                                            <  2:857713/150‑1 (MQ=255)
                                                                                                 cGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCcgcgc                                                      <  1:5124119/150‑1 (MQ=255)
                                                                                                    cgccgcCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCgcggc                                                   <  2:1735210/150‑1 (MQ=255)
                                                                                                     gccgccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCgcggcg                                                  <  2:5254163/150‑1 (MQ=255)
                                                                                                       cgccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCggcggc                                                >  2:1417173/1‑150 (MQ=255)
                                                                                                        gccAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCa                                               >  1:3294462/1‑150 (MQ=255)
                                                                                                           aaCTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCAtgt                                            >  2:3815319/1‑150 (MQ=255)
                                                                                                           aaCTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCAtgt                                            <  1:8425525/150‑1 (MQ=255)
                                                                                                            aCTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCAtgtg                                           <  1:1052816/150‑1 (MQ=255)
                                                                                                              tttGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                                         >  2:879801/1‑150 (MQ=255)
                                                                                                              tttGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATgtggt                                         <  1:1896077/150‑1 (MQ=255)
                                                                                                                tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGTGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                                       <  1:48888/150‑1 (MQ=255)
                                                                                                                tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                                       <  2:935483/150‑1 (MQ=255)
                                                                                                                tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                                       <  1:1171383/150‑1 (MQ=255)
                                                                                                                tGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGa                                       <  1:7750490/150‑1 (MQ=255)
                                                                                                                 gaaGATGTTCATGTGCAGCTATCCATAACTGCCGACAGGATAACCAGTCGGCGGTCATGGGGTTGCATAAGAGGAGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAg                                      >  2:834557/3‑150 (MQ=255)
                                                                                                                   aGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTa                                    <  1:6936798/150‑1 (MQ=255)
                                                                                                                      tgtgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtc                                 <  1:7983819/150‑1 (MQ=255)
                                                                                                                       gtgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtca                                <  2:1808761/150‑1 (MQ=255)
                                                                                                                       gtgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtca                                >  2:6377875/1‑150 (MQ=255)
                                                                                                                        tgcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtcat                               >  2:2642356/1‑150 (MQ=255)
                                                                                                                         gcatgtgcaGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACtcatc                              >  2:494197/1‑150 (MQ=255)
                                                                                                                                 aGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTc                      <  2:2154351/150‑1 (MQ=255)
                                                                                                                                    tctcCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGcc                   <  2:4671874/150‑1 (MQ=255)
                                                                                                                                      tcCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCtt                 >  1:4127348/1‑150 (MQ=255)
                                                                                                                                      tcCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCtt                 >  1:3742480/1‑150 (MQ=255)
                                                                                                                                       ccATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTg                >  1:8812005/1‑150 (MQ=255)
                                                                                                                                        cATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGc               >  2:7682795/1‑150 (MQ=255)
                                                                                                                                         atgaccgCCGACGGGTTATCCTGTCGGCAGTTATGGGGTTGCCTCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCa              >  2:7276758/7‑150 (MQ=38)
                                                                                                                                          tAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCAc             <  2:3457029/150‑1 (MQ=255)
                                                                                                                                             ccgCCGACGGGTTATCCTGTCGGCAGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGt          >  2:2839670/3‑150 (MQ=255)
                                                                                                                                              tGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc         <  1:4206001/150‑1 (MQ=255)
                                                                                                                                              tGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGtc         <  2:6855978/150‑1 (MQ=255)
                                                                                                                                                 cGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTcc      <  2:4374486/150‑1 (MQ=255)
                                                                                                                                                   aCAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTTTGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAg    <  1:8122839/150‑1 (MQ=255)
                                                                                                                                                     aGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGt  <  2:7323604/150‑1 (MQ=255)
                                                                                                                                                     aGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGt  <  2:1590268/150‑1 (MQ=255)
                                                                                                                                                     aGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGt  <  2:6130129/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGCCGATAGCCATTTGCTAGCATTTCTTCGACGGTGGCGCCCACCACCAGCCAGTCGATATCGCTGACAGGGCGGCCGAGCAAACGGTCGCGCACGGCGCCGCCAACTTTGTAGATGTGCATGTGCAGCTCTCCATAACTGCCGACAGGATAACCCGTCGGCGGTCATGGGGTTGCATCAGAGGTGGTGGATAACAGCCAGGTCCATGCGGCCATAGTCGGCGCTTTCGCTATGCTCGCCGCGCGGCGGCATGTGGTGAGTACTCATCACTTGCTCGCCTTGCACGGTCTCCAGGT  >  NC_002947/479347‑479645

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