Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,056,922 G→C 6.4% A466G (GCT→GGT) ‡ PP_0910 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,056,9220GC 6.4% 1163.7 / 50.0 345A466G (GCT→GGT) ‡PP_0910hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (167/156);  new base C (8/14);  total (175/170)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 355 total aligned reads displayed.

GGCGATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGCCC  >  NC_002947/1056773‑1057071
                                                                                                                                                     |                                                                                                                                                     
ggCGATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGAGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAg                                                                                                                                                       >  2:5209800/1‑150 (MQ=255)
  cGATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAgcg                                                                                                                                                     <  2:235868/150‑1 (MQ=255)
      ggcggcGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCgct                                                                                                                                                 <  1:4650340/150‑1 (MQ=255)
          gcggAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAGCGCGCTGCGTTGCTc                                                                                                                                             <  2:6181591/150‑1 (MQ=255)
          gcggAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAGCGCGCTGCGTTGCTc                                                                                                                                             <  1:3378751/150‑1 (MQ=255)
          gcggAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAGCGCGCTGCGTTGCTc                                                                                                                                             <  1:78220/150‑1 (MQ=255)
          gcggAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAGCGCGCTGCGTTGCTc                                                                                                                                             <  1:543837/150‑1 (MQ=255)
          gcggAGGTCACTGCTTGCGGGTCGTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAGCGCGCTGCGTTGCTc                                                                                                                                             <  2:4408102/150‑1 (MQ=255)
              aGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCa                                                                                                                                         <  1:6041092/150‑1 (MQ=255)
               ggTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCAt                                                                                                                                        >  2:5530169/1‑150 (MQ=255)
                gTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAGCGCGCTGCGTTGCTCgcgttc                                                                                                                                       <  1:3264791/150‑6 (MQ=255)
                  cACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGa                                                                                                                                     >  2:681198/1‑150 (MQ=255)
                  cACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGa                                                                                                                                     >  2:785218/1‑150 (MQ=255)
                    cTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCAtgttcc                                                                                                                                   <  2:1525518/150‑6 (MQ=255)
                    cTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGaaa                                                                                                                                   <  1:1733652/150‑1 (MQ=255)
                    cTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGaaa                                                                                                                                   <  1:2771298/150‑1 (MQ=255)
                     tGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAg                                                                                                                                  >  2:3417698/1‑150 (MQ=255)
                      gCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGt                                                                                                                                 <  1:2477488/150‑1 (MQ=255)
                       cTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTg                                                                                                                                >  1:2081305/1‑150 (MQ=255)
                        ttGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACTCAGCGCGCTGCTCGGCATCGAAAGTGa                                                                                                                               <  1:6624856/150‑1 (MQ=255)
                        ttGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGa                                                                                                                               >  1:1074163/1‑150 (MQ=255)
                        ttGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGa                                                                                                                               >  2:2339405/1‑150 (MQ=255)
                        ttGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGa                                                                                                                               <  1:3162400/150‑1 (MQ=255)
                        ttGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGa                                                                                                                               <  1:2240537/150‑1 (MQ=255)
                        ttGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGa                                                                                                                               <  1:2565039/150‑1 (MQ=255)
                          gCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGt                                                                                                                             <  1:1169114/150‑1 (MQ=255)
                            gggTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTAc                                                                                                                           <  1:280955/150‑1 (MQ=255)
                                 ttGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAAt                                                                                                                      >  1:4235480/1‑150 (MQ=255)
                                  tgtgGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAAtt                                                                                                                     <  2:6284647/150‑1 (MQ=255)
                                    tgGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGc                                                                                                                   >  2:5305928/1‑150 (MQ=255)
                                    tgGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGc                                                                                                                   >  2:3999277/1‑150 (MQ=255)
                                     ggCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCg                                                                                                                  >  2:1119729/1‑150 (MQ=255)
                                     ggCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTGTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCg                                                                                                                  >  2:4035920/1‑150 (MQ=255)
                                        gAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTgc                                                                                                               >  2:5835379/1‑150 (MQ=255)
                                        gAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTgc                                                                                                               >  1:5581181/1‑150 (MQ=255)
                                          gCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGATTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTgccg                                                                                                             <  2:5129963/150‑1 (MQ=255)
                                                aCCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTgccgccaacg                                                                                                       >  1:1858545/1‑146 (MQ=255)
                                                  cGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCt                                                                                                     >  1:6004600/1‑150 (MQ=255)
                                                    aGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTc                                                                                                   <  1:2888150/150‑1 (MQ=255)
                                                    aGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTc                                                                                                   <  2:2915262/150‑1 (MQ=255)
                                                     gggCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCg                                                                                                  <  1:6760529/150‑1 (MQ=255)
                                                      ggCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGc                                                                                                 <  1:4005073/150‑1 (MQ=255)
                                                      ggCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGc                                                                                                 <  2:1861939/150‑1 (MQ=255)
                                                              cggcgAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGc                                                                                         <  1:2669406/150‑1 (MQ=255)
                                                                gcgAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGACCGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCca                                                                                       >  2:3203653/1‑150 (MQ=255)
                                                                  gAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCcaca                                                                                     >  2:4748141/1‑150 (MQ=255)
                                                                      cTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAc                                                                                 >  2:1222927/1‑150 (MQ=255)
                                                                      cTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCAc                                                                                 <  1:3516492/150‑1 (MQ=255)
                                                                       tGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACg                                                                                >  1:4758372/1‑150 (MQ=255)
                                                                       tGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGATTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACg                                                                                >  2:1056774/1‑150 (MQ=255)
                                                                         gTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGc                                                                              <  2:5014720/150‑1 (MQ=255)
                                                                             acaacaGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCa                                                                          <  2:422325/150‑1 (MQ=255)
                                                                             acaacaGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCa                                                                          >  1:4592224/1‑150 (MQ=255)
                                                                               aacaGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAg                                                                        >  1:5203781/1‑150 (MQ=255)
                                                                                 caGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCa                                                                      <  1:3141348/150‑1 (MQ=255)
                                                                                  aGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAg                                                                     <  1:6152877/150‑1 (MQ=255)
                                                                                  aGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAg                                                                     >  2:4621170/1‑150 (MQ=255)
                                                                                   gTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAgg                                                                    <  2:592262/150‑1 (MQ=255)
                                                                                   gTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAgg                                                                    >  2:1735006/1‑150 (MQ=255)
                                                                                       gggTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGc                                                                <  2:6002175/150‑1 (MQ=255)
                                                                                        ggTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCg                                                               >  1:270223/1‑150 (MQ=255)
                                                                                            ggAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGgcg                                                           >  1:4706060/1‑150 (MQ=255)
                                                                                              aaCATGCCCGGCGCCTGTGCGCGTTTACAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGcgcg                                                         <  2:5782429/150‑1 (MQ=255)
                                                                                               aCATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCgg                                                        <  1:834953/150‑1 (MQ=255)
                                                                                               aCATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCgg                                                        <  1:5059994/150‑1 (MQ=255)
                                                                                                 aTGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGcc                                                      <  2:4701949/150‑1 (MQ=255)
                                                                                                 aTGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGcc                                                      <  1:1692430/150‑1 (MQ=255)
                                                                                                 aTGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGcc                                                      <  1:5871007/150‑1 (MQ=255)
                                                                                                 aTGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGcc                                                      <  2:6495541/150‑1 (MQ=255)
                                                                                                    cccGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTc                                                   >  2:5492573/1‑150 (MQ=255)
                                                                                                      cGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCAt                                                 >  1:1434822/1‑150 (MQ=255)
                                                                                                        gCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGc                                               <  2:6055758/150‑1 (MQ=255)
                                                                                                        gCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGc                                               <  1:6837104/150‑1 (MQ=255)
                                                                                                         cTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCa                                              <  2:6561859/150‑1 (MQ=255)
                                                                                                          tCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAg                                             >  2:1980288/1‑150 (MQ=255)
                                                                                                          tCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAg                                             >  2:2090586/1‑150 (MQ=255)
                                                                                                            cTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGa                                           <  2:2756156/150‑1 (MQ=255)
                                                                                                             tttGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGaa                                          <  1:5942528/150‑1 (MQ=255)
                                                                                                                 cgcgTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGAAGGAACTgg                                      >  1:4054856/1‑150 (MQ=255)
                                                                                                                  gcgTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGc                                     >  1:4997225/1‑150 (MQ=255)
                                                                                                                         aaGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGcc                              >  2:365339/1‑150 (MQ=255)
                                                                                                                           gTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCAc                            >  1:3469077/1‑150 (MQ=255)
                                                                                                                            tGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACt                           >  1:1941552/1‑150 (MQ=255)
                                                                                                                            tGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACt                           <  1:178880/150‑1 (MQ=255)
                                                                                                                            tGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACt                           <  1:3591732/150‑1 (MQ=255)
                                                                                                                            tGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACt                           >  1:488127/1‑150 (MQ=255)
                                                                                                                               gTTCAGGAACGCGGGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGAc                        <  1:1119749/150‑1 (MQ=255)
                                                                                                                                   aGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCACGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGtt                    <  1:2127543/150‑1 (MQ=255)
                                                                                                                                     gatgccgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTc                  >  1:3724295/6‑150 (MQ=255)
                                                                                                                                     gAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGgcg                                                           >  2:6804177/1‑109 (MQ=255)
                                                                                                                                      aaCGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCg                 >  2:2130900/1‑150 (MQ=255)
                                                                                                                                        cgcgAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGACCTGGCTATGGCCACTGACCGTTTCGCt               >  1:6563201/1‑150 (MQ=255)
                                                                                                                                         gcgAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTc              <  1:1009841/150‑1 (MQ=255)
                                                                                                                                         gcgAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTc              >  2:2627335/1‑150 (MQ=255)
                                                                                                                                          cgAGCAGCGCGCTGCGTTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCa             >  1:4106451/1‑150 (MQ=255)
                                                                                                                                          cgAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCa             <  1:5305953/150‑1 (MQ=255)
                                                                                                                                          cgAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCa             <  1:3239218/150‑1 (MQ=255)
                                                                                                                                             gCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGc          <  1:2199934/150‑1 (MQ=255)
                                                                                                                                                    aGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGcc   >  2:3601401/1‑150 (MQ=255)
                                                                                                                                                     gcgcgcTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGccc  <  1:1467831/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGATGGCGGCGGAGGTCACTGCTTGCGGGTCTTGTGGCGAGCAGGCACCGAGGGCGAGTGCGGCGAGGCTGGTGAACAACAGTTTGGGTCGGAACATGCCCGGCTCCTTTGCGCGTTTAAAGTGAGTTCAGGAACGCGAGCAACGCAGCGCGCTGCTCGGCATCGAAAGTGAGTACGAAATTGCGTGCCGCCTCGGCTTCGCCACCGTGCCACAGCACGGCCTCAAGCAGGTTGCGGGCGCGGCCGTCATGCAGGAACTGGCTATGGCCACTGACCGTTTCGCTCAGGCCGATGCCC  >  NC_002947/1056773‑1057071

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: ‑