Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 26,837:1 +A 15.4% intergenic (+33/+152) PP_0021 → / ← PP_0022 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294726,8371.A15.4% 436.7 / 89.2 130intergenic (+33/+152)PP_0021/PP_0022hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (60/50);  new base A (7/13);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 133 total aligned reads displayed.

TCGGCTTGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGG  >  NC_002947/26688‑26984
                                                                                                                                                      |                                                                                                                                                   
 cgGCTTGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTa                                                                                                                                                    <  2:3147179/150‑1 (MQ=255)
     ttGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACca                                                                                                                                                >  1:1266633/1‑150 (MQ=255)
     ttGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACca                                                                                                                                                >  1:1960023/1‑150 (MQ=255)
     ttGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTAACTGTAACca                                                                                                                                                >  1:1441177/1‑150 (MQ=255)
      tGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACcac                                                                                                                                               >  2:2343837/1‑150 (MQ=255)
         ttCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGt                                                                                                                                            >  2:5233837/1‑150 (MQ=255)
         ttCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGt                                                                                                                                            >  2:169823/1‑150 (MQ=255)
               tCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTa                                                                                                                                      <  1:127011/150‑1 (MQ=255)
               tCCCTTATCCGGAGGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGACGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGtta                                                                                                                                       <  1:5043503/150‑2 (MQ=39)
                 ccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACa                                                                                                                                    <  2:603588/150‑1 (MQ=255)
                 ccTTATCCGGAGGCCGGTTTATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGTGGAACTGGGCGATCTGCTGGAGCTGAATCAGGCCTGAGGAGGCATTTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACa                                                                                                                                    <  2:4190967/150‑1 (MQ=255)
                   ttctCCGGAGGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGt                                                                                                                                  <  1:2930792/147‑1 (MQ=255)
                   ttATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGt                                                                                                                                  <  1:5026278/150‑1 (MQ=255)
                   ttATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGt                                                                                                                                  <  1:4185793/150‑1 (MQ=255)
                   ttATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGt                                                                                                                                  <  1:2930852/150‑1 (MQ=255)
                    tATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTg                                                                                                                                 <  1:4549629/150‑1 (MQ=255)
                      tCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGaa                                                                                                                               <  2:4900470/150‑1 (MQ=255)
                         ggATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTTaaaaa                                                                                                                            >  1:609166/1‑150 (MQ=255)
                           aTGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCCGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGaaaaaaa                                                                                                                          <  2:4689264/150‑1 (MQ=38)
                             gCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGt                                                                                                                        <  1:3473948/150‑1 (MQ=255)
                             gCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACATTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGt                                                                                                                        <  2:748196/150‑1 (MQ=255)
                              ccGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTg                                                                                                                       <  2:1259381/150‑1 (MQ=255)
                              ccGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTg                                                                                                                       <  1:4305109/150‑1 (MQ=255)
                              ccGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGAGGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTg                                                                                                                       <  2:2068495/150‑1 (MQ=37)
                              ccGGTTGATCATCCTGATCGCAGAATGGAGAGAGGGAAGAGGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGCGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTg                                                                                                                       <  2:4805685/150‑1 (MQ=25)
                               cGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTgg                                                                                                                      >  2:1414592/1‑150 (MQ=255)
                                ggTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGt                                                                                                                     <  2:565827/150‑1 (MQ=255)
                                ggTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGt                                                                                                                     <  2:1767736/150‑1 (MQ=255)
                                     atcatcCTTATCGCATAATGGATAGTGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGTCCTTTTGTAAACCACTTTTTTTCACTGTAAGCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACa                                                                                                                <  2:4317003/150‑1 (MQ=38)
                                     atcatcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACa                                                                                                                <  2:1562303/150‑1 (MQ=255)
                                      tcatcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACaa                                                                                                               <  2:2723308/150‑1 (MQ=255)
                                        atcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAACAAAGTGGTTTACAACg                                                                                                             >  2:4306567/1‑150 (MQ=255)
                                        atcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACg                                                                                                             >  2:2716266/1‑150 (MQ=255)
                                          ccTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTATTTTCACTGTAACCACAGTGGATTACAGTGAAAAAAAGTGGTTTACAACGTa                                                                                                           <  2:2364030/150‑1 (MQ=255)
                                           cTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAc                                                                                                                                                  >  2:1323543/1‑110 (MQ=255)
                                           cTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAc                                                                                                                                                  <  1:1323543/110‑1 (MQ=255)
                                             tATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCt                                                                                                        >  2:3832249/1‑150 (MQ=255)
                                                cGCAGAATGGATAGAGTGAAGATGGCCCGGGCGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTGTACAGTGAAAAACAGTGGTTTACAACGTAGCTACg                                                                                                     >  2:3154501/1‑150 (MQ=255)
                                                 gCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGa                                                                                                    >  2:2889245/1‑150 (MQ=255)
                                                  cAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGac                                                                                                   <  1:2605543/150‑1 (MQ=255)
                                                     aaTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                          <  2:1942034/124‑1 (MQ=255)
                                                     aaTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaaaaaaa                                                                                                                          >  1:1942034/1‑124 (MQ=255)
                                                              agTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGc                                                                                       <  2:255925/150‑1 (MQ=255)
                                                                tGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCa                                                                                     <  1:2627675/150‑1 (MQ=255)
                                                                        gCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAg                                                                             <  1:285898/150‑1 (MQ=255)
                                                                         cccGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAgg                                                                            >  2:2432263/1‑150 (MQ=255)
                                                                          ccGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGc                                                                           >  2:4004063/1‑150 (MQ=255)
                                                                                    cTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctc                                                                 <  2:1535539/150‑1 (MQ=255)
                                                                                    cTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctc                                                                 <  1:4522934/150‑1 (MQ=255)
                                                                                    cTGGGCGATCTGCGGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctc                                                                 <  2:4166314/150‑1 (MQ=255)
                                                                                     tGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctcc                                                                >  1:4910387/1‑150 (MQ=255)
                                                                                         cGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAt                                                            >  2:3718863/1‑150 (MQ=255)
                                                                                          gATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                           <  2:4704585/150‑1 (MQ=255)
                                                                                          gATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                           <  1:2391229/150‑1 (MQ=255)
                                                                                           ctCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTGACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTg                                                          <  2:4838133/149‑1 (MQ=255)
                                                                                            tCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGa                                                         <  1:2343837/150‑1 (MQ=255)
                                                                                               gctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTc                                                      <  1:3017787/150‑1 (MQ=255)
                                                                                               gctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTc                                                      <  1:2595259/150‑1 (MQ=255)
                                                                                                 tgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGt                                                    >  1:3793326/1‑150 (MQ=255)
                                                                                                     gCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAg                                                >  1:3500617/1‑150 (MQ=255)
                                                                                                           tCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGc                                          >  1:3458155/1‑150 (MQ=255)
                                                                                                           tCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGc                                          >  2:4722916/1‑150 (MQ=255)
                                                                                                             aGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGt                                        >  1:1861106/1‑150 (MQ=255)
                                                                                                             aGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGt                                        >  2:1046292/1‑150 (MQ=255)
                                                                                                              ggCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                       <  2:4153037/150‑1 (MQ=255)
                                                                                                              ggCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                       <  2:3793326/150‑1 (MQ=255)
                                                                                                              ggCCTGAGGAGGCCTGTTGTAAACCACTTTTCTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                       <  2:3119272/150‑1 (MQ=255)
                                                                                                                   gaggagGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAaca                                  <  2:5316167/150‑1 (MQ=255)
                                                                                                                   gaggagGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAaca                                  <  2:734561/150‑1 (MQ=255)
                                                                                                                      gaggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCt                               >  1:2257779/1‑150 (MQ=255)
                                                                                                                        ggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                             <  2:877975/150‑1 (MQ=255)
                                                                                                                        ggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                             <  2:1960023/150‑1 (MQ=255)
                                                                                                                        ggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                             <  2:1960238/150‑1 (MQ=255)
                                                                                                                        ggCCTGTTGTAAACCACTTTGTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                             <  1:2916508/150‑1 (MQ=255)
                                                                                                                         gCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtg                            >  2:2596109/1‑150 (MQ=255)
                                                                                                                         gCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtg                            <  1:2889245/150‑1 (MQ=255)
                                                                                                                         gCCTGTTGTAAACCACTTTATTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtg                            <  2:4051188/150‑1 (MQ=255)
                                                                                                                            tgttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAGCGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgctg                         >  2:1119182/1‑150 (MQ=255)
                                                                                                                            tgttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCGGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgctg                         >  1:909615/1‑150 (MQ=255)
                                                                                                                                   aaCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCgggg                  >  1:5427648/1‑150 (MQ=255)
                                                                                                                                     ccACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGcc                >  2:1245674/1‑150 (MQ=255)
                                                                                                                                     ccACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGcc                >  2:2720396/1‑150 (MQ=255)
                                                                                                                                       aCTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCtt              >  1:2229567/1‑150 (MQ=255)
                                                                                                                                         tttttttCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGc            >  1:230577/1‑150 (MQ=255)
                                                                                                                                         tttttttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGc            >  2:2935600/1‑150 (MQ=255)
                                                                                                                                         tttttttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGc            >  2:2935809/1‑150 (MQ=255)
                                                                                                                                              ttCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACATAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgc       >  2:1647156/1‑150 (MQ=38)
                                                                                                                                              ttCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgc       >  1:5435572/1‑150 (MQ=255)
                                                                                                                                               tCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgcc      >  1:4364611/1‑150 (MQ=255)
                                                                                                                                               tCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgcc      >  1:5395354/1‑150 (MQ=255)
                                                                                                                                                cACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAg                                                                             <  2:4014439/78‑1 (MQ=18)
                                                                                                                                                cACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAg                                                                             >  1:4014439/1‑78 (MQ=18)
                                                                                                                                                cACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccg     >  1:419520/1‑150 (MQ=255)
                                                                                                                                                 aCTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                             >  1:5090194/1‑125 (MQ=255)
                                                                                                                                                 aCTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAt                             <  2:5090194/125‑1 (MQ=255)
                                                                                                                                                 aCTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccgc    >  1:552787/1‑150 (MQ=255)
                                                                                                                                                 aCTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccgc    >  1:1967895/1‑150 (MQ=255)
                                                                                                                                                   tGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCgg  >  2:4629787/1‑150 (MQ=255)
                                                                                                                                                   tGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCgg  <  2:609166/150‑1 (MQ=255)
                                                                                                                                                   tGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCgg  >  1:3981114/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCGGCTTGGTTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGG  >  NC_002947/26688‑26984

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: