Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,475,003 C→A 9.2% intergenic (‑776/+25) algD ← / ← PP_1289 GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,475,0030CA 9.2% 407.5 / 24.7 130intergenic (‑776/+25)algD/PP_1289GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (56/62);  new base A (2/10);  total (58/72)
Fisher's exact test for biased strand distribution p-value = 6.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 133 total aligned reads displayed.

TCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGATGA  >  NC_002947/1474854‑1475152
                                                                                                                                                     |                                                                                                                                                     
tCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTc                                                                                                                                                       >  2:3947499/1‑150 (MQ=255)
  aaaaaaGGTATGAACGGTCTTCTGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAg                                                                                                                                                     <  2:4564960/150‑1 (MQ=255)
  aaaaaaGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAg                                                                                                                                                     >  1:3312798/1‑150 (MQ=255)
  aaaaaaGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAg                                                                                                                                                     >  1:2378027/1‑150 (MQ=255)
    aaaaGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGcc                                                                                                                                                   >  1:1362925/1‑150 (MQ=255)
         gTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGaa                                                                                                                                              <  2:633911/150‑1 (MQ=255)
           aTGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAaga                                                                                                                                            <  2:4284365/150‑1 (MQ=255)
           aTGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAaga                                                                                                                                            <  1:3639802/150‑1 (MQ=255)
           aTGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAaga                                                                                                                                            <  1:3829876/150‑1 (MQ=255)
            tGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAagag                                                                                                                                           <  1:4837107/150‑1 (MQ=255)
            tGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAagag                                                                                                                                           <  2:4234589/150‑1 (MQ=255)
            tGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTGTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAagag                                                                                                                                           <  2:4762171/150‑1 (MQ=255)
             gAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAgg                                                                                                                                          >  1:274404/1‑150 (MQ=255)
              accGGACTTATGGCGTTTGGCGATAAAGCGAGGTGCGACGCATTTTTGGCGACGGATTCGAGGCGTGAAGGCGGAGGGCCGCGTGCCCAGTTGCAGGCCAGAGCCGGATCACGCCGGCTGACCTCTACGTGGCTGCAGCCGCGACGAggg                                                                                                                                         >  2:30094/3‑150 (MQ=255)
              aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAggg                                                                                                                                         <  2:3989663/150‑1 (MQ=255)
               aCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGc                                                                                                                                        >  1:1225144/1‑150 (MQ=255)
               aCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGc                                                                                                                                        <  1:4755378/150‑1 (MQ=255)
                cGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGCCCTCTTCGCGGCTGAAGCCGCGAAGAggtca                                                                                                                                       <  2:2888250/150‑4 (MQ=255)
                 ggTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGt                                                                                                                                      >  2:4266626/1‑150 (MQ=255)
                  gTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTg                                                                                                                                     >  2:3656993/1‑150 (MQ=255)
                        tGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGc                                                                                                                               >  1:1295495/1‑150 (MQ=255)
                         ggCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCt                                                                                                                              <  2:4402613/150‑1 (MQ=255)
                            gTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCAc                                                                                                                           >  1:2958548/1‑150 (MQ=255)
                                  cGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCcgcg                                                                                                                     <  1:3548844/150‑1 (MQ=255)
                                  cGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCcgcg                                                                                                                     <  2:3269149/150‑1 (MQ=255)
                                  cGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCcgcg                                                                                                                     <  1:3549545/150‑1 (MQ=255)
                                   gATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCgg                                                                                                                    >  2:1791287/1‑150 (MQ=255)
                                    aTAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCggg                                                                                                                   <  2:3485598/150‑1 (MQ=255)
                                    aTAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCggg                                                                                                                   <  2:2487959/150‑1 (MQ=255)
                                     tAAAGGGAGGTGTGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCgggg                                                                                                                  <  2:288105/150‑1 (MQ=255)
                                      aaaGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGt                                                                                                                 <  2:1225144/150‑1 (MQ=255)
                                          ggAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGAAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATc                                                                                                             >  1:1295977/1‑150 (MQ=255)
                                              gTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCgcgc                                                                                                         <  1:4897443/150‑1 (MQ=255)
                                               tGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCa                                                                                                        <  1:3928471/150‑1 (MQ=255)
                                                 cGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAgg                                                                                                      <  2:3459173/150‑1 (MQ=255)
                                                   aCGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGaa                                                                                                    >  2:4975125/1‑150 (MQ=255)
                                                   aCGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGaa                                                                                                    >  1:1730680/1‑150 (MQ=255)
                                                   aCGGATTTTTGGCGTCGGATTCGTGGCATGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGag                                                                                                    <  2:3201016/150‑2 (MQ=255)
                                                     ggATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAaca                                                                                                  >  2:3843958/1‑150 (MQ=255)
                                                       aTTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACAcc                                                                                                <  1:4169311/150‑1 (MQ=255)
                                                        tttttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCa                                                                                               >  1:4197537/1‑150 (MQ=255)
                                                        tttttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCa                                                                                               >  1:4197457/1‑150 (MQ=255)
                                                         ttttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCaa                                                                                              <  2:4364003/150‑1 (MQ=255)
                                                          tttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCaaa                                                                                             >  2:920674/1‑150 (MQ=255)
                                                          tttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCaaa                                                                                             <  1:807646/150‑1 (MQ=255)
                                                          tttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCaaa                                                                                             <  1:1854072/150‑1 (MQ=255)
                                                                gTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCt                                                                                       <  2:2732758/150‑1 (MQ=255)
                                                                gTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCt                                                                                       >  1:3079497/1‑150 (MQ=255)
                                                                ctctgATTCGTGGCGTGATGGCGGATTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCt                                                                                       <  2:3410999/146‑1 (MQ=255)
                                                                 tCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTg                                                                                      <  2:2520363/150‑1 (MQ=255)
                                                                 tCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTg                                                                                      >  1:3705803/1‑150 (MQ=255)
                                                                 tCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTg                                                                                      >  2:2969680/1‑150 (MQ=255)
                                                                 tCGGAATCGTGGCGAGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCTGCCGCGAAGTGGGCGTGACAGGCTCACTCCGCGGGGTGTACGCGCAGGAACACCAAATGGTCTg                                                                                      >  2:160809/1‑150 (MQ=255)
                                                                  cGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTgg                                                                                     <  2:1874245/150‑1 (MQ=255)
                                                                   ggATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGc                                                                                    >  1:196209/1‑150 (MQ=255)
                                                                        cGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGa                                                                               <  1:4351260/150‑1 (MQ=255)
                                                                             cGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCt                                                                          <  1:4266626/150‑1 (MQ=255)
                                                                               tGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCACATGGTCTGGCTTGGACTGCTCg                                                                        >  2:3686878/1‑150 (MQ=255)
                                                                                gATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCgg                                                                       >  2:4403717/1‑150 (MQ=255)
                                                                                  tGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCa                                                                     <  1:4545893/150‑1 (MQ=255)
                                                                                      ggTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgt                                                                 >  1:5007727/1‑150 (MQ=255)
                                                                                        tttGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgtag                                                               <  1:3378099/150‑1 (MQ=255)
                                                                                             cctctGCCCTGCAACAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgtagtagta                                                          >  1:429214/5‑150 (MQ=255)
                                                                                             cctctGCCCTGCAACAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgtagtagta                                                          >  1:5087476/5‑150 (MQ=255)
                                                                                              gCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAacac                                                                                                 >  1:2596022/1‑110 (MQ=255)
                                                                                              gCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAacac                                                                                                 <  2:2596022/110‑1 (MQ=255)
                                                                                                gTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTAccc                                                       <  1:4524029/150‑1 (MQ=255)
                                                                                                 tGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTAcccc                                                      >  2:4631319/1‑150 (MQ=255)
                                                                                                    ccTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGa                                                   <  1:3214699/150‑1 (MQ=255)
                                                                                                    ccTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGa                                                   <  1:362472/150‑1 (MQ=255)
                                                                                                    ccTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGa                                                   <  2:655412/150‑1 (MQ=255)
                                                                                                        ttGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATcgcg                                                                                                          <  2:4593410/91‑1 (MQ=255)
                                                                                                        ttGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATcgcg                                                                                                          >  1:4593410/1‑91 (MQ=255)
                                                                                                              ggCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACt                                         >  2:1860872/1‑150 (MQ=255)
                                                                                                              ggCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACt                                         >  2:1860449/1‑150 (MQ=255)
                                                                                                                cAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGc                                       <  2:1232058/150‑1 (MQ=255)
                                                                                                                 agagGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCt                                      <  1:4770386/150‑1 (MQ=255)
                                                                                                                  gagGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCtt                                     <  2:2794831/150‑1 (MQ=255)
                                                                                                                    ggCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGa                                   <  2:1248251/150‑1 (MQ=255)
                                                                                                                     gCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAg                                  >  2:3614104/1‑150 (MQ=255)
                                                                                                                       ggATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGTCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCt                                <  2:3030002/150‑1 (MQ=255)
                                                                                                                       ggATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCt                                <  1:2879956/150‑1 (MQ=255)
                                                                                                                        gATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTg                               >  2:2022870/1‑150 (MQ=255)
                                                                                                                                        ctctTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATg               <  1:3947499/150‑1 (MQ=255)
                                                                                                                                        ctctTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATg               <  2:196209/150‑1 (MQ=255)
                                                                                                                                         tctTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgc              >  1:2193630/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgc              >  1:2193393/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCAGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgc              >  2:4745906/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGCTTCTGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgcgc            >  1:4685417/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgcgc            >  1:4507998/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgcgc            >  2:1582627/1‑150 (MQ=255)
                                                                                                                                             cgcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctg          <  2:3079497/150‑1 (MQ=255)
                                                                                                                                              gcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctgc         >  2:220166/1‑150 (MQ=255)
                                                                                                                                              gcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGCAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctgc         >  1:2181859/1‑150 (MQ=255)
                                                                                                                                               cgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctgct        <  1:3567802/150‑1 (MQ=255)
                                                                                                                                               cgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctgct        <  2:274404/150‑1 (MQ=255)
                                                                                                                                                   ttCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGAt    <  1:920674/150‑1 (MQ=255)
                                                                                                                                                   ttCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGAt    <  2:1981250/150‑1 (MQ=255)
                                                                                                                                                     cAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGATGa  >  1:643686/1‑150 (MQ=255)
                                                                                                                                                     cAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGATGa  >  1:4423451/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGATGA  >  NC_002947/1474854‑1475152

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