Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,353,312 A→C 9.7% intergenic (+21/+55) PP_4706 → / ← PP_4707 hypothetical protein/OsmY‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,353,3120AC 9.7% 792.6 / 58.5 256intergenic (+21/+55)PP_4706/PP_4707hypothetical protein/OsmY‑like protein
Reads supporting (aligned to +/- strand):  ref base A (119/112);  new base C (9/16);  total (128/128)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 256 total aligned reads displayed.

TGATGCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAAC  >  NC_002947/5353165‑5353459
                                                                                                                                                   |                                                                                                                                                   
tgatgCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACg                                                                                                                                                   <  1:8989717/150‑1 (MQ=255)
   tgCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGc                                                                                                                                                <  1:1186867/150‑1 (MQ=255)
     cggcgTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCa                                                                                                                                              <  2:7431627/150‑1 (MQ=255)
       gcgTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGACTGCTTCGAAAGCGATGTGATCGCGCGCCACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAgg                                                                                                                                            >  1:3590243/1‑150 (MQ=255)
        cgTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGCCAggg                                                                                                                                           <  1:7188693/150‑1 (MQ=255)
        cgTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAggg                                                                                                                                           <  1:1060452/150‑1 (MQ=255)
         gTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGc                                                                                                                                          <  2:1334393/150‑1 (MQ=255)
           aTGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGCCAGGGCtt                                                                                                                                        <  1:6483729/150‑1 (MQ=255)
           aTGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGCCAGGGCtt                                                                                                                                        <  2:1814935/150‑1 (MQ=255)
              aCGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGa                                                                                                                                     >  1:3312633/1‑150 (MQ=255)
               cgcgCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGCCAGGGCtttttc                                                                                                                                    <  2:4132419/150‑5 (MQ=255)
                   cATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGaaa                                                                                                                                <  2:6545816/150‑1 (MQ=255)
                   cATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGaaa                                                                                                                                <  2:497676/150‑1 (MQ=255)
                       cccAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTgg                                                                                                                            >  1:2262243/1‑150 (MQ=255)
                          aGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCa                                                                                                                         >  2:7887359/1‑150 (MQ=255)
                             cTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAAt                                                                                                                      <  1:296230/150‑1 (MQ=255)
                               gcgcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATcc                                                                                                                    >  1:3328950/1‑150 (MQ=255)
                                cgcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCa                                                                                                                   >  1:8041050/1‑150 (MQ=255)
                                 gcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAt                                                                                                                  >  2:2283157/1‑150 (MQ=255)
                                 gcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAt                                                                                                                  >  2:6280942/1‑150 (MQ=255)
                                   aTTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAttt                                                                                                                <  1:3975563/150‑1 (MQ=255)
                                   aTTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAttt                                                                                                                <  1:7856446/150‑1 (MQ=255)
                                     tttCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCTGTGTATTGAGCGCAAGAAAAAGCCCTGGCACTGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTgg                                                                                                              <  1:6229259/150‑1 (MQ=255)
                                     tttCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTgg                                                                                                              <  1:5624348/150‑1 (MQ=255)
                                      ttCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGCGGCCAAATGGATTTGGc                                                                                                             <  1:3638640/150‑1 (MQ=255)
                                        cgcgACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCg                                                                                                           <  2:1120206/150‑1 (MQ=255)
                                              aGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGa                                                                                                     >  2:5084767/1‑150 (MQ=255)
                                                ccGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGACAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAACTCCATTTGGCCGACGGGACa                                                                                                   <  2:4840193/150‑1 (MQ=255)
                                                   tCCATTGGTGGAGCCACTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGAATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTc                                                                                                >  2:4531253/1‑150 (MQ=255)
                                                    ccATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTcc                                                                                               <  1:5841124/150‑1 (MQ=255)
                                                      aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGACAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCtt                                                                                             <  1:6287068/150‑1 (MQ=255)
                                                      aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCtt                                                                                             >  2:129448/1‑150 (MQ=255)
                                                      aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCtt                                                                                             >  2:8765529/1‑150 (MQ=255)
                                                       ttGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTa                                                                                            >  1:8342273/1‑150 (MQ=255)
                                                          gtggAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTc                                                                                         <  1:6710845/150‑1 (MQ=255)
                                                          gtggAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTc                                                                                         <  2:4956497/150‑1 (MQ=255)
                                                          gtggAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTc                                                                                         <  1:2095830/150‑1 (MQ=255)
                                                               gCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGcgt                                                                                    >  2:5448739/1‑148 (MQ=255)
                                                                  tctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCa                                                                                 <  2:2084389/150‑1 (MQ=255)
                                                                  tctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCa                                                                                 >  2:1556826/1‑150 (MQ=255)
                                                                   ctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCtt                                                                                                                                        <  1:428387/94‑1 (MQ=255)
                                                                   ctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCtt                                                                                                                                        >  2:428387/1‑94 (MQ=255)
                                                                    tGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCACGAAAATGCCCTGGCACGGGCCGGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAgg                                                                               >  2:256556/1‑150 (MQ=255)
                                                                    tGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAgg                                                                               >  2:7469507/1‑150 (MQ=255)
                                                                       tCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTa                                                                                            >  2:8364004/1‑134 (MQ=255)
                                                                         cTTGAGCAGGAATGCTTCGACAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGACAGTGGCCAAAGCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTc                                                                          <  2:3093475/150‑1 (MQ=255)
                                                                                cggAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACAtt                                                                   <  2:7107328/149‑1 (MQ=255)
                                                                                aGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACAtt                                                                   >  1:1901707/1‑150 (MQ=255)
                                                                                aGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACAtt                                                                   <  1:4825089/150‑1 (MQ=255)
                                                                                 ggAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTa                                                                  >  1:3690544/1‑150 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGCCGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           >  1:885654/1‑150 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           >  1:4929351/1‑150 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           <  1:8738059/150‑1 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           >  2:6990110/1‑150 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           >  1:2458014/1‑150 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           >  1:1439114/1‑150 (MQ=255)
                                                                                        ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                           >  1:8479792/1‑150 (MQ=255)
                                                                                                   tgtgATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGa                                                <  2:5388934/150‑1 (MQ=255)
                                                                                                   tgtgATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGa                                                <  2:8913577/150‑1 (MQ=255)
                                                                                                   ggtgATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAACAGCCCTGGCACGGGCCAGGGCTTCGATGACAGTGGCCAACTCCATTTGGCCGACGGGACAGTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGa                                                <  2:8143338/149‑1 (MQ=255)
                                                                                                    gtgATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGaa                                               <  2:6475473/150‑2 (MQ=255)
                                                                                                    gtgATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGAt                                               <  2:1858979/150‑1 (MQ=255)
                                                                                                     tgATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCTGACATTACGCTCTTGACGCCTTTGATg                                              >  2:9093299/1‑150 (MQ=255)
                                                                                                      gATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTGTGATGt                                             >  2:5244013/1‑150 (MQ=255)
                                                                                                       aTCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCTGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGtt                                            <  1:1316670/150‑1 (MQ=255)
                                                                                                       aTCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGtt                                            <  2:3201040/150‑1 (MQ=255)
                                                                                                         cgcgcgCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCt                                          <  1:5382643/150‑1 (MQ=255)
                                                                                                           cgcgcAACCCGTTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTg                                        <  2:6875314/150‑1 (MQ=255)
                                                                                                           cgcgcAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTg                                        <  1:2068488/150‑1 (MQ=255)
                                                                                                            gcgcAACCCTGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGCCAGGGCtt                                                                                                                                        >  1:696922/1‑53 (MQ=25)
                                                                                                               cAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTaa                                    <  2:3986867/150‑1 (MQ=255)
                                                                                                                aaCCCGGTGTATTGAGCGCAAGAAATAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAAt                                   >  2:630530/1‑150 (MQ=255)
                                                                                                                aaCCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACTGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATACCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAAt                                   >  2:7307883/1‑150 (MQ=255)
                                                                                                                  cccGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATgc                                 >  2:3675336/1‑150 (MQ=255)
                                                                                                                  cccGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATgc                                 >  1:2862149/1‑150 (MQ=255)
                                                                                                                      gtgtATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCAc                             <  1:6138504/150‑1 (MQ=255)
                                                                                                                           ttGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCtt                        <  1:5910656/150‑1 (MQ=255)
                                                                                                                              aGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTg                     >  1:7025360/1‑150 (MQ=255)
                                                                                                                                 gcAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCtt                  >  1:8436382/1‑150 (MQ=255)
                                                                                                                                   aaGACAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAACTCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                <  1:7836867/150‑1 (MQ=255)
                                                                                                                                   aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                <  1:1945402/150‑1 (MQ=255)
                                                                                                                                   aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                >  2:1645576/1‑150 (MQ=255)
                                                                                                                                   aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                >  2:8672542/1‑150 (MQ=255)
                                                                                                                                   aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                >  1:1567127/1‑150 (MQ=255)
                                                                                                                                    aGACAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTCTTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCt               <  2:8539666/150‑1 (MQ=255)
                                                                                                                                    aGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCt               <  2:2756508/150‑1 (MQ=255)
                                                                                                                                     gAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCtt                                         <  1:6703892/123‑1 (MQ=255)
                                                                                                                                     atcgaaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTg              >  1:1253453/5‑150 (MQ=255)
                                                                                                                                      tcgaaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGc             >  1:4421308/4‑150 (MQ=255)
                                                                                                                                      aaaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGc             <  1:2638360/150‑1 (MQ=255)
                                                                                                                                       aaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCg            >  2:7839338/1‑150 (MQ=255)
                                                                                                                                        gaaGCCCTGGCCCGTGCCAGGGCTTTGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGa           >  2:5815471/2‑150 (MQ=255)
                                                                                                                                           gCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGTCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTc        <  2:8400328/150‑1 (MQ=255)
                                                                                                                                           gCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTc        >  2:6872291/1‑150 (MQ=255)
                                                                                                                                           gCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTc        <  2:4706039/150‑1 (MQ=255)
                                                                                                                                            cccTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCg       <  2:569890/150‑1 (MQ=255)
                                                                                                                                            cccTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCg       <  1:5478057/150‑1 (MQ=255)
                                                                                                                                            cccTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCg       <  1:7533530/150‑1 (MQ=255)
                                                                                                                                              cTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGt     >  2:3487467/1‑150 (MQ=255)
                                                                                                                                                 gCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAAc  <  2:20707/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGATGCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAAC  >  NC_002947/5353165‑5353459

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