Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 188,613 T→C 5.5% L241P (CTG→CCG)  PP_0165 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947188,6130TC 5.5% 1061.9 / 29.8 330L241P (CTG→CCG) PP_0165GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (159/151);  new base C (6/12);  total (165/165)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 332 total aligned reads displayed.

ATTCCTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACC  >  NC_002947/188464‑188762
                                                                                                                                                     |                                                                                                                                                     
  tCCTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGa                                                                                                                                                     <  2:8657752/150‑1 (MQ=255)
    cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                   >  1:7148512/1‑150 (MQ=255)
    cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                   >  1:6306403/1‑150 (MQ=255)
    cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                   >  1:1670895/1‑150 (MQ=255)
    cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                   >  2:8199109/1‑150 (MQ=255)
    cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                   <  1:511495/150‑1 (MQ=255)
     tCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCg                                                                                                                                                  >  1:3441965/1‑150 (MQ=255)
       aGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCggg                                                                                                                                                >  2:6753217/1‑150 (MQ=255)
        gCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGc                                                                                                                                               <  2:7220815/150‑1 (MQ=255)
         ccTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCt                                                                                                                                              <  2:856530/150‑1 (MQ=255)
            gCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTc                                                                                                                                           >  2:3321805/1‑150 (MQ=255)
            gCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTc                                                                                                                                           <  2:5192689/150‑1 (MQ=255)
             ccAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCg                                                                                                                                          >  1:5995637/1‑150 (MQ=255)
              cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTcct                                                                                                                                         <  1:3778284/150‑3 (MQ=255)
              cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcct                                                                                                                                         <  1:195353/150‑3 (MQ=255)
              cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcct                                                                                                                                         <  1:4777948/150‑3 (MQ=255)
              cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcct                                                                                                                                         <  1:6858244/150‑3 (MQ=255)
               agagTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGaa                                                                                                                                        >  1:318197/1‑150 (MQ=255)
                gagTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAc                                                                                                                                       <  1:8300291/150‑1 (MQ=255)
                gagTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAc                                                                                                                                       <  2:5421178/150‑1 (MQ=255)
                  gTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcctggta                                                                                                                                     <  2:2693662/150‑7 (MQ=255)
                   ttGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAccgc                                                                                                                                    <  2:5276426/150‑1 (MQ=255)
                    tGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAccgcc                                                                                                                                   <  2:1120543/150‑1 (MQ=255)
                     gCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAccgccg                                                                                                                                  >  2:4154242/1‑150 (MQ=255)
                        acgcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTa                                                                                                                               >  1:3707902/1‑150 (MQ=255)
                              gcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGa                                                                                                                         <  1:5359349/150‑1 (MQ=255)
                                cGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGaaa                                                                                                                       <  1:6077032/150‑1 (MQ=255)
                                      tACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGc                                                                                                                 >  2:3206156/1‑150 (MQ=255)
                                       aCGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCt                                                                                                                >  1:7194894/1‑150 (MQ=255)
                                        cgccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTg                                                                                                               <  2:73199/150‑1 (MQ=255)
                                        cgccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTg                                                                                                               >  1:8769461/1‑150 (MQ=255)
                                         gccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGa                                                                                                              >  2:1661582/1‑150 (MQ=255)
                                         gccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGa                                                                                                              <  1:3330131/150‑1 (MQ=255)
                                         gccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGa                                                                                                              >  1:5208402/1‑150 (MQ=255)
                                              gTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAAt                                                                                                         >  2:2969233/1‑150 (MQ=255)
                                                 gTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTTGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCgg                                                                                                      <  1:4820644/150‑1 (MQ=255)
                                                    cAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTa                                                                                                   >  1:5029054/1‑150 (MQ=255)
                                                     aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTaa                                                                                                  >  2:3962659/1‑150 (MQ=255)
                                                     aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTaa                                                                                                  <  2:6396021/150‑1 (MQ=255)
                                                     aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTaa                                                                                                  >  1:1317230/1‑150 (MQ=255)
                                                       gCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGc                                                                                                >  2:7175434/1‑150 (MQ=255)
                                                          aTGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAAc                                                                                             <  1:6649821/150‑1 (MQ=255)
                                                             aaCCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTg                                                                                          >  1:1441224/1‑150 (MQ=255)
                                                             aaCCAGATGGTCGAGAAACTCAAAGCCCTGATCACCGAGCAGGCCGAGCGCAGTGAGTGGCTGCGAGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTg                                                                                          >  1:6438278/1‑150 (MQ=255)
                                                              aCCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTgg                                                                                         >  1:1085722/1‑150 (MQ=255)
                                                                cAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGaa                                                                                       <  1:6561967/150‑1 (MQ=255)
                                                                   aTGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAc                                                                                    <  2:8410020/150‑1 (MQ=255)
                                                                   aTGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAc                                                                                    >  1:116505/1‑150 (MQ=255)
                                                                     ggTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgc                                                                                  >  1:4528868/1‑150 (MQ=255)
                                                                     ggTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgc                                                                                  <  2:1988480/150‑1 (MQ=255)
                                                                       tCGAGAAACTCAAAGCCCTGTTCACCGAGCATGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTAACAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATTCAGCTGAACAATCGGGTAAGCAACCTGTAAGAcgcgc                                                                                >  1:4853047/1‑150 (MQ=255)
                                                                       tCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgc                                                                                <  2:3492939/150‑1 (MQ=255)
                                                                         gagaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcgc                                                                              >  1:2567562/1‑150 (MQ=255)
                                                                         gagaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                          >  1:1249980/1‑89 (MQ=255)
                                                                          agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCTCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcgcg                                                                             <  1:4462537/150‑1 (MQ=255)
                                                                            aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGcc                                                                           <  2:7748561/150‑1 (MQ=255)
                                                                            aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGcc                                                                           <  1:6984098/150‑1 (MQ=255)
                                                                             aaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCg                                                                          <  1:2870345/150‑1 (MQ=255)
                                                                               cTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGt                                                                        >  2:694322/1‑150 (MQ=255)
                                                                                    aGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACACTCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                   <  2:3611911/150‑1 (MQ=255)
                                                                                         tGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                              <  1:7257183/150‑1 (MQ=255)
                                                                                         tGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                              >  2:3154341/1‑150 (MQ=255)
                                                                                          gTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGtt                                                             >  2:1948122/1‑150 (MQ=255)
                                                                                             cACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgct                                                          <  1:7625777/150‑1 (MQ=255)
                                                                                              aCCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgctg                                                         >  1:8253624/1‑150 (MQ=255)
                                                                                               ccGAGCAGGTCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgctgc                                                        <  1:7378395/150‑1 (MQ=255)
                                                                                                cGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCg                                                       >  1:5849426/1‑150 (MQ=255)
                                                                                                 gAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCgg                                                      <  1:4028890/150‑1 (MQ=255)
                                                                                                         cGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAggg                                              <  2:6306403/150‑1 (MQ=255)
                                                                                                          gAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGCACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGc                                             <  2:44710/150‑1 (MQ=255)
                                                                                                                  tGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCctggctgg                                     >  1:5410973/1‑150 (MQ=255)
                                                                                                                  tGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCctggctgg                                     >  2:5995479/1‑150 (MQ=255)
                                                                                                                   ggcctgCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCtggctggc                                    <  2:5013882/145‑1 (MQ=255)
                                                                                                                   gAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCtggctggc                                    <  2:5385955/150‑1 (MQ=255)
                                                                                                                    aGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGcc                                   <  2:7378955/150‑1 (MQ=255)
                                                                                                                    aGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGcc                                   <  2:144591/150‑1 (MQ=255)
                                                                                                                         cTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAc                              <  1:8640254/150‑1 (MQ=255)
                                                                                                                          tgcgtgcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcg                             >  1:4195403/1‑150 (MQ=255)
                                                                                                                            cgtgcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgcc                           >  2:4323747/1‑150 (MQ=255)
                                                                                                                             gtgcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                          <  2:5778818/150‑1 (MQ=255)
                                                                                                                                  gAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCtt                     >  1:5586401/1‑150 (MQ=255)
                                                                                                                                       cTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcgg                >  2:8617555/1‑150 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACTGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               <  1:3089106/150‑1 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               >  1:2323425/1‑150 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               >  2:6057219/1‑150 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               >  2:6235222/1‑150 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               >  1:2394391/1‑150 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               <  2:7582777/150‑1 (MQ=255)
                                                                                                                                        tACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc               >  2:7015354/1‑150 (MQ=255)
                                                                                                                                         aCCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAg                                                <  1:5720936/116‑1 (MQ=255)
                                                                                                                                          ccAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCa             >  2:1316395/1‑150 (MQ=255)
                                                                                                                                          ccAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCa             <  2:2342397/150‑1 (MQ=255)
                                                                                                                                           ttcgACAGCCCGGTCAGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAg            >  1:2739830/4‑150 (MQ=255)
                                                                                                                                              gACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGt         >  1:7779321/1‑150 (MQ=255)
                                                                                                                                              gACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCAAggcgggcgttca         <  2:668860/150‑8 (MQ=255)
                                                                                                                                                   ccTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGa    >  1:7479215/1‑150 (MQ=255)
                                                                                                                                                   ccTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGa    >  2:4671896/1‑150 (MQ=255)
                                                                                                                                                   ccTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGa    >  1:8321090/1‑150 (MQ=255)
                                                                                                                                                   ccTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGa    >  2:3443377/1‑150 (MQ=255)
                                                                                                                                                    cTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc   <  1:923683/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCCTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACC  >  NC_002947/188464‑188762

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: