Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,435,079 A→C 6.0% intergenic (+39/‑145) PP_3924 → / → PP_3925 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,435,0790AC 6.0% 795.2 / 26.2 228intergenic (+39/‑145)PP_3924/PP_3925hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (90/124);  new base C (5/9);  total (95/133)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

GGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATG  >  NC_002947/4434930‑4435226
                                                                                                                                                     |                                                                                                                                                    
         cgcggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCgg                                                                                                                                             <  1:5808674/149‑1 (MQ=255)
           cggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGtt                                                                                                                                          <  1:6114834/150‑1 (MQ=255)
              cgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttttt                                                                                                                                       <  1:7315524/150‑1 (MQ=255)
              cgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttttt                                                                                                                                       <  1:4971213/150‑1 (MQ=255)
              cgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGttttt                                                                                                                                       <  2:1056863/150‑1 (MQ=255)
               gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGtttttc                                                                                                                                      <  2:7319373/150‑2 (MQ=255)
                 aaCCGACGGCATGCAGGCCTTCTTGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttttt                                                                                                                                    <  2:6094766/150‑1 (MQ=255)
                  aCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCGTCGGTTTTTTTTg                                                                                                                                   >  2:5636576/1‑150 (MQ=255)
                   ccGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGc                                                                                                                                  >  2:767659/1‑150 (MQ=255)
                    cGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGcc                                                                                                                                 <  2:3455036/150‑1 (MQ=255)
                    cGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGcc                                                                                                                                 <  2:2030501/150‑1 (MQ=255)
                      aCGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCt                                                                                                                               <  2:663249/150‑1 (MQ=255)
                      aCGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCt                                                                                                                               <  1:4790393/150‑1 (MQ=255)
                       cGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTg                                                                                                                              <  1:6039378/150‑1 (MQ=255)
                        ggCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGc                                                                                                                             <  2:5321625/150‑1 (MQ=255)
                         gCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCa                                                                                                                            <  1:6171654/150‑1 (MQ=255)
                             gCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCg                                                                                                                        <  1:1990393/150‑1 (MQ=255)
                                   cTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAg                                                                                                                  <  2:6738521/150‑1 (MQ=255)
                                     tCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTg                                                                                                                >  2:3087346/1‑150 (MQ=255)
                                       aTGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGAt                                                                                                              <  2:5478888/150‑1 (MQ=255)
                                        tGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATc                                                                                                             <  2:1435758/150‑1 (MQ=255)
                                         gggTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATcc                                                                                                            <  1:5086085/150‑1 (MQ=255)
                                           gTGGCACGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccac                                                                                                          <  2:5019874/150‑1 (MQ=255)
                                            tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacc                                                                                                         >  2:5992847/1‑150 (MQ=255)
                                            tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacc                                                                                                         >  2:6273671/1‑150 (MQ=255)
                                            tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacc                                                                                                         >  2:3972962/1‑150 (MQ=255)
                                             ggCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacca                                                                                                        <  1:300360/150‑1 (MQ=255)
                                               cAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGc                                                                                                      >  1:3960192/1‑150 (MQ=255)
                                                aaGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCa                                                                                                     >  1:3955022/1‑150 (MQ=255)
                                                              aGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAt                                                                                       >  2:70326/1‑150 (MQ=255)
                                                               ggCGACATGATGCTGCCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAtt                                                                                      >  2:6332561/1‑150 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAttt                                                                                     >  1:447835/1‑150 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAttt                                                                                     <  1:7434946/150‑1 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCGTCGGTTTTTTTTGCCCTGCATTAGGTGCAGATATCCACCAGCAACACCCTCAACAAttt                                                                                     >  2:4718044/1‑150 (MQ=255)
                                                                  gACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTgg                                                                                   >  2:1441112/1‑150 (MQ=255)
                                                                   aCATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGc                                                                                  >  1:6390074/1‑150 (MQ=255)
                                                                     atgatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGcac                                                                                >  1:1303679/1‑150 (MQ=255)
                                                                       gattCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGAc                                                                                                                                                     >  1:407510/1‑80 (MQ=255)
                                                                       gattCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGAc                                                                                                                                                     <  2:407510/80‑1 (MQ=255)
                                                                         tgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATgg                                                                            <  1:2099732/150‑1 (MQ=255)
                                                                               cGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGACAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGAt                                                                      <  2:4335808/150‑1 (MQ=255)
                                                                                 aTGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTg                                                                    >  2:3542789/1‑150 (MQ=255)
                                                                                   gAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGtt                                                                  <  1:3125056/150‑1 (MQ=255)
                                                                                   gAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGtt                                                                  >  1:3408245/1‑150 (MQ=255)
                                                                                    aaaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 <  2:3219854/150‑1 (MQ=255)
                                                                                    aaaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 <  2:3112102/150‑1 (MQ=255)
                                                                                    aaaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 <  2:6351774/150‑1 (MQ=255)
                                                                                      aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCa                                                               >  1:3622640/1‑150 (MQ=255)
                                                                                      aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGAATTCGGTGCAGTGATCAACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCa                                                               >  2:6226902/1‑150 (MQ=255)
                                                                                       cTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCaa                                                              <  1:2698700/150‑1 (MQ=255)
                                                                                        tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAAc                                                             >  1:4362118/1‑150 (MQ=255)
                                                                                             gAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTgg                                                        <  1:1386610/150‑1 (MQ=255)
                                                                                             gAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTgg                                                        <  2:4400489/150‑1 (MQ=255)
                                                                                                 tGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGa                                                    >  1:5633280/1‑150 (MQ=255)
                                                                                                 tGTTCCAGTCCTGAGGGCTGGGAATACCTGATTATCGAAAAAACGAAGCATTAAGGATTAGGTTTTTTTTTTCAATGCATTCGGTGCATTTATCAACCAGCAAAAACCTCACCTATTTGGCACATGGGACGATTTTTCCAAAACTGGCgg                                                     >  2:1587076/1‑150 (MQ=255)
                                                                                                    tCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTgg                                                        >  2:4138884/1‑143 (MQ=255)
                                                                                                    tCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCa                                                 <  2:4213115/150‑1 (MQ=255)
                                                                                                      ccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAtt                                               >  1:5260279/1‑150 (MQ=255)
                                                                                                      ccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAtt                                               <  2:1923356/150‑1 (MQ=255)
                                                                                                      ccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAtt                                               <  1:4946618/150‑1 (MQ=255)
                                                                                                       cGTCCTTAGTGCTTGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCTGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                              <  2:4118101/150‑1 (MQ=255)
                                                                                                       cGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                              <  2:5782069/150‑1 (MQ=255)
                                                                                                       cGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                              >  1:703082/1‑150 (MQ=255)
                                                                                                       cGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                              <  2:2898251/150‑1 (MQ=255)
                                                                                                       cGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                              <  2:2031991/150‑1 (MQ=255)
                                                                                                        gTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAg                                             <  1:6805209/150‑1 (MQ=255)
                                                                                                        gTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGATATTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAg                                             <  1:7213760/150‑1 (MQ=255)
                                                                                                            tGAGGGCTGGGAATACCGGATGATCGAAAATCCGAAGCCTGACGGCTTCGGTTTTGTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAACTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGAttc                                         >  1:1899465/1‑149 (MQ=255)
                                                                                                            tGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTa                                         >  2:6595266/1‑150 (MQ=255)
                                                                                                                ggCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCa                                     >  1:4135084/1‑150 (MQ=255)
                                                                                                                ggCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCa                                     >  2:4525404/1‑150 (MQ=255)
                                                                                                                 gCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCaa                                    <  1:7113569/150‑1 (MQ=255)
                                                                                                                 gCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCaa                                    <  1:974710/150‑1 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  <  1:6454546/150‑1 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  >  1:3517435/1‑150 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  <  1:5139803/150‑1 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  <  1:2591009/150‑1 (MQ=255)
                                                                                                                      gAACACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTa                               <  2:604508/150‑1 (MQ=255)
                                                                                                                         tACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACtt                            <  2:3067952/150‑1 (MQ=255)
                                                                                                                          aCCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 <  2:1744592/112‑1 (MQ=255)
                                                                                                                          aCCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 >  1:1744592/1‑112 (MQ=255)
                                                                                                                            cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCt                         >  2:2746112/1‑150 (MQ=255)
                                                                                                                              gatgatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTAt                       <  2:5260279/150‑1 (MQ=255)
                                                                                                                              gatgatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTAt                       >  2:5668361/1‑150 (MQ=255)
                                                                                                                               atgatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGtt                                                                  >  1:2364881/1‑106 (MQ=255)
                                                                                                                               atgatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGtt                                                                  <  2:2364881/106‑1 (MQ=255)
                                                                                                                                 gatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga                    <  1:425458/150‑1 (MQ=255)
                                                                                                                                 gatCGAAAAACCGAAGCCTGACGGCGTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga                    >  1:184243/1‑150 (MQ=255)
                                                                                                                                  atCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGCCGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCgag                   >  1:5867924/1‑150 (MQ=255)
                                                                                                                                  atCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAGTAGATTACCCAATAGTACTTGCTATCgag                   >  1:3475943/1‑150 (MQ=255)
                                                                                                                                   tCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCgaga                  >  2:1212939/1‑150 (MQ=255)
                                                                                                                                    cGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGaa                 <  2:21928/150‑1 (MQ=255)
                                                                                                                                     gAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAAt                >  1:3890570/1‑150 (MQ=255)
                                                                                                                                     gAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAAt                >  1:5208500/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAAt                >  1:5450556/2‑150 (MQ=255)
                                                                                                                                     aaaaaaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCAACAGCATCACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCAGGCGGACCCTTAGATTACCCAATAGTACTTGCTATCGCGAAt                >  1:1209909/2‑150 (MQ=255)
                                                                                                                                          aCCGAAGCCTGACGGCTTCGGTTTTTATTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCAGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGa           >  1:5631782/1‑150 (MQ=255)
                                                                                                                                           ccGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGaa          <  1:3845949/150‑1 (MQ=255)
                                                                                                                                             gAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGa        >  1:2553242/1‑150 (MQ=255)
                                                                                                                                                gCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGc     <  2:1793400/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATG  >  NC_002947/4434930‑4435226

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