Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 233,832 A→T 9.2% intergenic (+20/+43) PP_0180 → / ← PP_0181 cytochrome c family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947233,8320AT 9.2% 338.4 / 20.1 128intergenic (+20/+43)PP_0180/PP_0181cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (48/64);  new base T (4/8);  total (56/72)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 142 total aligned reads displayed.

TGAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCC  >  NC_002947/233683‑233981
                                                                                                                                                     |                                                                                                                                                     
tGAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCca                                                                                                                                                       >  2:855609/1‑150 (MQ=255)
tGAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCca                                                                                                                                                       <  2:3537018/150‑1 (MQ=255)
 gAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcac                                                                                                                                                      <  2:5394226/150‑1 (MQ=255)
  aGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcaca                                                                                                                                                     <  2:5103430/150‑1 (MQ=255)
   gTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAg                                                                                                                                                    <  1:3729375/150‑1 (MQ=255)
   gTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAg                                                                                                                                                    >  2:2488361/1‑150 (MQ=255)
   gTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAg                                                                                                                                                    <  1:3729339/150‑1 (MQ=255)
   gTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAg                                                                                                                                                    <  2:5299556/150‑1 (MQ=255)
    ttCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcccacc                                                                                                                                                   >  1:193491/1‑148 (MQ=255)
    ttCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                   >  1:3543638/1‑150 (MQ=255)
     tCGACTGGCTGGGGAGTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAggg                                                                                                                                                  <  2:4668454/150‑1 (MQ=255)
          tggctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACcgc                                                                                                                                             <  2:212058/150‑1 (MQ=255)
             ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGTCTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  1:678362/150‑1 (MQ=255)
             ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  2:1434823/150‑1 (MQ=255)
             ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  1:3858602/150‑1 (MQ=255)
             ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  2:2113760/150‑1 (MQ=255)
             ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  1:1972453/150‑1 (MQ=255)
               ggggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCCCAGGGACCCCGCaac                                                                                                                                        >  1:2772559/1‑149 (MQ=255)
                 ggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGt                                                                                                                                      <  1:4746625/150‑1 (MQ=255)
                  gATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtc                                                                                                                                     <  1:694120/150‑1 (MQ=255)
                   aTTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtct                                                                                                                                    <  2:5046427/150‑1 (MQ=255)
                   aTTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtct                                                                                                                                    <  2:3796328/150‑1 (MQ=255)
                     tCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGGGCGGTGCCACAGGGACCCCCCAAAGtctaa                                                                                                                                  >  1:106993/1‑148 (MQ=255)
                      cATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAg                                                                                                                                 <  1:4515843/150‑1 (MQ=255)
                       aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGa                                                                                                                                >  2:4991333/1‑150 (MQ=255)
                       aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGa                                                                                                                                <  2:5383798/150‑1 (MQ=255)
                       aTGCCGAGGCCTACTCGCGGTCGGCGCACGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGa                                                                                                                                <  2:4347848/150‑1 (MQ=255)
                           cGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACTGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCtt                                                                                                                            <  1:812876/150‑1 (MQ=255)
                              tGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGATCTGAGACCTTgcg                                                                                                                         <  1:3413923/150‑1 (MQ=255)
                                  tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAg                                                                                                                     <  1:5313167/150‑1 (MQ=255)
                                   aCTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGc                                                                                                                    <  2:674092/150‑1 (MQ=255)
                                      cGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTg                                                                                                                 <  1:2388378/150‑1 (MQ=255)
                                       gCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATAATGGGCCAGCCGGa                                                                                                                >  2:8249/1‑150 (MQ=255)
                                        cTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAt                                                                                                               <  1:2715815/150‑1 (MQ=255)
                                             ggCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCa                                                                                                          <  1:2177309/150‑1 (MQ=255)
                                             ggCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCa                                                                                                          <  1:2177269/150‑1 (MQ=255)
                                                gcAAGCGTTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggc                                                                                                       <  2:1433011/150‑1 (MQ=255)
                                                  aaGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGTGACCGCGCAAGGTCTAAGTACTTTCGCAACCCTGATGCACAggcag                                                                                                     >  2:4632029/1‑148 (MQ=255)
                                                    gCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAGCATCAggcttgc                                                                                                   <  2:2903299/150‑1 (MQ=255)
                                                    gCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAGCATCAggcttgc                                                                                                   <  2:2903399/150‑1 (MQ=255)
                                                       gTGGCCTTGCTGGCCATCGTGTTCCTTTCCGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCtgt                                                                                                <  2:2287627/150‑1 (MQ=255)
                                                            cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGTTGCTCAGGCGGGCGGTGgcgc                                                                                           >  2:432637/1‑150 (MQ=255)
                                                               gCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                        <  2:1864709/150‑1 (MQ=255)
                                                               gCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                        <  1:4900674/150‑1 (MQ=255)
                                                                cTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggc                                                                                                       >  2:4985121/1‑134 (MQ=255)
                                                                 tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                      <  2:1585221/150‑1 (MQ=255)
                                                                 tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                      <  2:2391972/150‑1 (MQ=255)
                                                                 tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                      <  1:983620/150‑1 (MQ=255)
                                                                 tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGGGCCCCAGGGCCCCCCCCACGTGACAGATAATTCGGAAGCCGGATTCTCAGGCTGGCAGTGGCGccgggt                                                                                      >  1:4235682/1‑146 (MQ=255)
                                                                   gCCGTCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAg                                                                                                                                                    <  1:3517427/86‑1 (MQ=255)
                                                                        cGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCCCGCCAGGTCTCAGATCGTGCGCAAGCCAGATGCTCAGGCGGGCGGGGGCGCCGTgagacagtt                                                                               >  1:4984023/1‑145 (MQ=255)
                                                                         gtgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTc                                                                              <  1:4632029/150‑1 (MQ=255)
                                                                           gtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTg                                                                            <  2:4254317/150‑1 (MQ=255)
                                                                             tCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGGCCGCGCCAGGTCTCAGATCGTGCGCAAGCCTGATGCTCAGGCTGGCGGGGGCGCCGTGAGCAAGTTCTGtg                                                                          >  1:3547374/1‑149 (MQ=255)
                                                                              ccTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGGGCCACACGGACCCCCCCAGGTGTCAGAGCTTTCGCCAGCCCGATGCTCAGGCGGGCGGTGGCGCCGTGTGAATGTTgtgttg                                                                         >  1:3807674/1‑148 (MQ=255)
                                                                               cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTACCACAGGGACCGCGCAAGGTCTAAGATCTTGCGCAAGCCAGATGCTCAGGCTGGCGGTGGCGCCTTGAGAAAGTTCTGTAGc                                                                        >  2:1744339/1‑150 (MQ=255)
                                                                                tttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGCCCGCGCCAGGTCTCAGATCTTTCGCAACCCAGATGCGCAGGCGGGCGGTGGCGCCCTGAGAATGTTATTTAGCa                                                                       >  2:630726/1‑150 (MQ=255)
                                                                                    cGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                   <  2:1856748/150‑1 (MQ=255)
                                                                                    cGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                   <  2:5062512/150‑1 (MQ=255)
                                                                                     ggCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTTCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCATGAGAATGTTCTGTAgatcggaa                                                                  >  2:2639267/1‑143 (MQ=255)
                                                                                      gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTGCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCAGATTCTCAGGCTGGCTGTGGCGCCTTGAGCATGTTCTGTAGCAATTGGa                                                                 >  2:3303545/1‑150 (MQ=255)
                                                                                        cGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGaca                                                               <  2:1225661/150‑1 (MQ=255)
                                                                                           tCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCCAGCCTGagcatc                                                                                                           >  2:341557/1‑98 (MQ=255)
                                                                                             gTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTAACACATCTTGCGCAAGCCTGATTCTCAGGCAGGCTGTGGCGCCTTGAGAAAGTTCTGTAGCAATTGGACACTTTc                                                          >  2:4280096/1‑150 (MQ=255)
                                                                                                   ggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGc                                                    <  2:5282909/150‑1 (MQ=255)
                                                                                                    gcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCg                                                   >  2:1780596/1‑150 (MQ=255)
                                                                                                      cGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTGCCAAATGGACCGCGCAAGGTCTGAGATGTTGCGCAAGCCTTATGCTCAGGCTGGCTGTGGCGCCTTGGTAATGTTCTGTATCAATTTGACACTTTAAACGGCGta                                                 >  2:2237613/1‑149 (MQ=255)
                                                                                                      cGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTc                                                 <  1:3117540/150‑1 (MQ=255)
                                                                                                       aaccctagcAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCa                                                <  2:4149346/143‑1 (MQ=255)
                                                                                                        aaaaGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGGCCCACAGGGACCCCGCACGGTCTCAGATCTTGCGCAAGCCGGATGCGCAGGCGGGGTGTGGCGCCGTGAGGAGGGTCGGTAGCAAGTGGACACATTCAACGGCggcca                                               >  2:2735714/1‑146 (MQ=255)
                                                                                                        aaaaGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCCCGGTCCCACAGGGACCGCGCAAGGTCGCAGATCGTGCGCAAGCCAGATGCGCCCGCGGGCGGGGGCGCCGTGCGAAGGTTCGGTAGCCAGTGGGCACGTTCAACGGCGGCaa                                               >  2:3324565/1‑149 (MQ=255)
                                                                                                              gcAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCa                                         <  2:2616642/150‑1 (MQ=255)
                                                                                                                   gctgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGTGAATGTTCTGTAGCAATTGGACACTTTCCACGGCGTCATGACCCAGGCTg                                    >  2:4014478/1‑150 (MQ=255)
                                                                                                                   gctgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTg                                    <  1:3801321/150‑1 (MQ=255)
                                                                                                                   gctgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTg                                    <  2:4346308/150‑1 (MQ=255)
                                                                                                                       cGGCCAACTGACTGGAGCCTGCCCCGTCCCACAGGGCCCCCGCCAGGTCTCAGATCTTGCGCCACCCAGATGCGCAGGCGGGCGGTGGCGCCGGGAGGATGTTCTGTAGCACTTGGACACGTTCAACGGCGGCCTGACCCAGGCTGGTCa                                >  2:3646608/1‑150 (MQ=255)
                                                                                                                            aaCTGACTGGAGCCTGCGCCGTACCACCGGGGCCCCGCAAGGGAGAAGGGCTTGCGCCATCCGGATGGGCACGCGGGCGGTGGGGCCGTTAGAATGTTCTGTAGCAATTGGACACATTAAACGGCGTGAGGTCCCAGGCGGGTCAGGTAg                           >  2:4155925/1‑150 (MQ=255)
                                                                                                                                gactgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgccgc                       <  2:495973/150‑1 (MQ=255)
                                                                                                                                   tgGAGCCTGCGCGGTCCCACACGGACCCCCCAAGGTCTCACATCTTGCGCAAACCAGAAGCACACGCAGGCTGTGGCGCCGTGAGAAAGTTATGTAGCAAAAGGACACATTAAACCGCGTCAATAACCACGCTGGTCAGGTAGCCGCCaa                    >  2:1666417/1‑149 (MQ=255)
                                                                                                                                    ggAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATc                   <  1:4991333/150‑1 (MQ=255)
                                                                                                                                      accttGCGCGGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCg                 >  1:1105231/5‑150 (MQ=255)
                                                                                                                                       gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCgg                <  1:1608424/150‑1 (MQ=255)
                                                                                                                                        cttGCGCGGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGc               >  2:2557330/3‑150 (MQ=255)
                                                                                                                                         ttGCGCGGTCCCTGTGGGACCGCGCAGGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGc               >  2:3005834/2‑149 (MQ=255)
                                                                                                                                          tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTg             >  2:1635849/1‑150 (MQ=255)
                                                                                                                                            cgcgGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGAt           <  1:2738361/150‑1 (MQ=255)
                                                                                                                                             gcgGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATc          >  1:2524689/1‑150 (MQ=255)
                                                                                                                                                  cccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa     <  2:106993/150‑1 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCTAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  2:2436397/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  2:4280479/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  2:3140729/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  2:1241398/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  1:1817324/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  2:39415/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTGGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTCGCCGCCATCCGGCTGATCGGTCAg    >  2:167937/1‑150 (MQ=255)
                                                                                                                                                    cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc   >  1:5199090/1‑150 (MQ=255)
                                                                                                                                                    cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc   >  1:1201080/1‑150 (MQ=255)
                                                                                                                                                     acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGcc  >  2:3144378/1‑150 (MQ=255)
                                                                                                                                                     acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGcc  >  2:3767833/1‑150 (MQ=255)
                                                                                                                                                     acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGcc  >  1:1466436/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCC  >  NC_002947/233683‑233981

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: