Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 155,546 A→C 5.6% intergenic (‑251/+32) yjbB ← / ← ygdQ Na+/Pi‑cotransporter/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947155,5460AC 5.6% 306.9 / 7.2 106intergenic (‑251/+32)yjbB/ygdQNa+/Pi‑cotransporter/transport protein
Reads supporting (aligned to +/- strand):  ref base A (47/53);  new base C (4/2);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 107 total aligned reads displayed.

CTGTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGG  >  NC_002947/155398‑155690
                                                                                                                                                      |                                                                                                                                                
  gTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAgg                                                                                                                                               <  1:1444781/150‑1 (MQ=255)
     gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGcc                                                                                                                                            >  2:1105575/1‑150 (MQ=255)
      atctcTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGCTGGCGCAGGCTTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGccc                                                                                                                                           <  2:652285/146‑1 (MQ=255)
       cAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCt                                                                                                                                          <  2:1850274/150‑1 (MQ=255)
           ttCGCATGGCAGTCATTTATTTGAAACTTGATTGGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCtagat                                                                                                                                      >  1:652285/1‑146 (MQ=255)
           ttCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                      <  1:304865/150‑1 (MQ=255)
               cATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGccc                                                                                                                                           >  1:1418538/1‑141 (MQ=255)
               cATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGccc                                                                                                                                           <  2:1418538/141‑1 (MQ=255)
               cATGGCAGTCATTTATTTGAAACAGATCTTGAAAAGGCCGAGTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                  <  1:122392/150‑1 (MQ=255)
               cATGGCAGTCATTTATTTGAAACAGATCTTAAAACGAATGACTGCCTGTGGGGGCTGCTTCGCGTGCAAGCCATCACAGTGACAACGCGCAGTAAAAAAAAAAAAAAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                  <  2:2232746/150‑1 (MQ=255)
               cATGGCAGTAATTCTTTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                  <  2:134505/150‑1 (MQ=255)
                  ggCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgctgc                                                                                                                               <  1:695517/150‑1 (MQ=255)
                   gCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCa                                                                                                                              >  1:2147137/1‑150 (MQ=255)
                       tCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGcc                                                                                                                                            >  1:972124/1‑132 (MQ=255)
                         atttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAGGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCt                                                                                                                        >  1:1548633/1‑150 (MQ=255)
                          tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                  >  1:1428699/1‑139 (MQ=255)
                          tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                  <  2:1428699/139‑1 (MQ=255)
                          tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                       <  1:1818956/150‑1 (MQ=255)
                          tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                       >  1:2285376/1‑150 (MQ=255)
                          tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                       >  2:1967214/1‑150 (MQ=255)
                          tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                       <  2:625163/150‑1 (MQ=255)
                           ttatctGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGAAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGGACTTCCTTg                                                                                                                      <  2:788383/150‑1 (MQ=255)
                           ttatctGAAACAGATCTTGAAAAGAATGACTCCCCGTGGGGGCTGCCTCGCGGGCAAGCCACCACAGTGACAACGCCCAGGAAAAAACGCAAAGAAAAAAGGCATTAAAAAGGGCCCGGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                      <  1:1647665/150‑1 (MQ=255)
                             atttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTtgtg                                                                                                                                   <  2:935979/135‑1 (MQ=255)
                             atttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTtgtg                                                                                                                                   >  1:935979/1‑135 (MQ=255)
                             atttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGc                                                                                                                    <  2:696613/150‑1 (MQ=255)
                                 gAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                      >  2:1857336/1‑128 (MQ=255)
                                 gAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                      <  1:1857336/128‑1 (MQ=255)
                                 gAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTa                                                                                                                <  2:567364/150‑1 (MQ=255)
                                  aaaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                               >  1:1876979/1‑150 (MQ=255)
                                  aaaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                               <  2:2147137/150‑1 (MQ=255)
                                  aaaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                               <  2:1541005/150‑1 (MQ=255)
                                   aaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAtt                                                                                                              >  2:1729076/1‑150 (MQ=255)
                                       aGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGAc                                                                                                           <  2:1876979/150‑1 (MQ=255)
                                         aTCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGAccc                                                                                                         <  1:226146/150‑1 (MQ=255)
                                               aaaCGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATg                                                                                                   >  2:2271231/1‑150 (MQ=255)
                                                aaCGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATg                                                                                                   >  2:718098/1‑150 (MQ=255)
                                                   gAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCt                                                                                               >  1:217879/1‑150 (MQ=255)
                                                         cTGCCCGTGGGGGCTGCTACGCGGGAAAGCCACCACAGTAAAAAAAAAAAAAAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCa                                                                                         <  2:1042400/150‑1 (MQ=255)
                                                          tGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTgg                                                                                                                     <  1:1462882/121‑1 (MQ=255)
                                                          tGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTgg                                                                                                                     >  2:1462882/1‑121 (MQ=255)
                                                             ccGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCAGCGGGTGGGCCTTATTGACCCGGAATGTCCTTGCGCAGttt                                                                                     >  1:122366/1‑150 (MQ=255)
                                                             ccGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGttt                                                                                     <  1:464579/150‑1 (MQ=255)
                                                             ccGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGtt                                                                                      <  2:2069598/149‑1 (MQ=255)
                                                             ccGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGtt                                                                                      >  1:2069598/1‑149 (MQ=255)
                                                                   gggCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGAc                                                                                                           >  1:127787/1‑122 (MQ=255)
                                                                   gggCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGAc                                                                                                           <  2:127787/122‑1 (MQ=255)
                                                                        gctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                      <  2:828480/106‑1 (MQ=255)
                                                                        gctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCAt                                                                          >  2:63787/1‑150 (MQ=255)
                                                                          ttCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAAAGTCCTTGCGGGGATTCATCGGtcccaca                                                                        >  2:1885005/1‑144 (MQ=255)
                                                                               gggCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCtt                                                                   <  1:1885005/150‑1 (MQ=255)
                                                                                ggCAAGCCACCACAGTGACAACGCGCAGGAAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCtt                                                                   <  1:718098/150‑1 (MQ=255)
                                                                                 gCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGc                                                                 <  1:63787/150‑1 (MQ=255)
                                                                                  cAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGcc                                                                <  2:2285376/150‑1 (MQ=255)
                                                                                   aaGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGcc                                                                                                                                            >  2:1015831/1‑73 (MQ=39)
                                                                                   aaGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGcc                                                                                                                                            <  1:1015831/73‑1 (MQ=255)
                                                                                   aaGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCa                                                               >  2:947861/1‑150 (MQ=255)
                                                                                     gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTg                                                                  <  1:2108332/145‑1 (MQ=255)
                                                                                     gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTg                                                                  >  2:2108332/1‑145 (MQ=255)
                                                                                     gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATc                                                             >  2:1730087/1‑150 (MQ=255)
                                                                                     gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCATTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTTCCATc                                                             >  1:130993/1‑150 (MQ=255)
                                                                                          cacaGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGc                                                        <  1:1729076/150‑1 (MQ=255)
                                                                                           acaGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                       <  2:1489369/150‑1 (MQ=255)
                                                                                               tGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGcc                                                   >  2:774735/1‑150 (MQ=255)
                                                                                               tGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGcc                                                   >  2:1625134/1‑150 (MQ=255)
                                                                                                   aaCGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATgg                                               >  1:269851/1‑150 (MQ=255)
                                                                                                   aaCGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATgg                                               >  2:2048731/1‑150 (MQ=255)
                                                                                                     cgcgCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCg                                             >  2:1880971/1‑150 (MQ=255)
                                                                                                      gcgcAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCgg                                            <  2:1733649/150‑1 (MQ=255)
                                                                                                        gcAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                       >  1:1934903/1‑137 (MQ=255)
                                                                                                        gcAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTg                                          <  1:890541/150‑1 (MQ=255)
                                                                                                                    cGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTTCGGGCCATGGCGGTGCGCATGCgggt                               <  2:217879/150‑3 (MQ=255)
                                                                                                                     gCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATg                              >  2:1445922/1‑150 (MQ=255)
                                                                                                                         aGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGa                          >  2:2160406/1‑150 (MQ=255)
                                                                                                                          gAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGAt                         <  1:947861/150‑1 (MQ=255)
                                                                                                                          gAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGAt                         >  1:611066/1‑150 (MQ=255)
                                                                                                                           aaaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATg                        >  2:1142279/1‑150 (MQ=255)
                                                                                                                            aaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATgg                       <  1:2048731/150‑1 (MQ=255)
                                                                                                                             aaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGc                      <  1:1967214/150‑1 (MQ=255)
                                                                                                                                aGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGa                                                                        <  1:65/97‑1 (MQ=255)
                                                                                                                                aGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGa                                                                        >  2:65/1‑97 (MQ=255)
                                                                                                                                aGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTc                   <  2:2178901/150‑1 (MQ=255)
                                                                                                                                    aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                      <  2:1890679/47‑1 (MQ=255)
                                                                                                                                    aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                      >  1:1890679/1‑47 (MQ=255)
                                                                                                                                    aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCAcc               <  1:1445922/150‑1 (MQ=255)
                                                                                                                                       aacaaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGcc            >  1:1555308/4‑150 (MQ=255)
                                                                                                                                       aaaaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                         <  1:350373/121‑1 (MQ=255)
                                                                                                                                       aaaaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                         >  2:350373/1‑121 (MQ=255)
                                                                                                                                           aGGGCCCGGCGAAGGGCCCTTGTTGTTCTGCTCTTCCTTGGCCTTATTGACCCGGAATGTCCTTGTGCAGTTTCACTGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGCTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCg        <  2:1125176/150‑1 (MQ=255)
                                                                                                                                            gggCCCTTCGCAGGGCCTTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCga       >  2:1191025/1‑150 (MQ=255)
                                                                                                                                            gggCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGCTCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCg                                  >  1:707123/1‑123 (MQ=255)
                                                                                                                                            gggCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCg                                  <  2:707123/123‑1 (MQ=255)
                                                                                                                                            gggCCCTGCGAAGGGCCCTTGTTGTGCTGCCTTTTAATGGGAGCAGCGACCCGCAATTTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCga       <  1:2271231/150‑1 (MQ=255)
                                                                                                                                               cccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCGTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGaa    >  2:2018910/1‑150 (MQ=255)
                                                                                                                                               cccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTCATGGCTTCCACCGCCAGCGAGaa    <  1:1880971/150‑1 (MQ=255)
                                                                                                                                                ccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                                  >  1:2096168/1‑71 (MQ=255)
                                                                                                                                                ccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                                  <  2:2096168/71‑1 (MQ=255)
                                                                                                                                                 cTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAgg  >  2:632332/1‑150 (MQ=255)
                                                                                                                                                  tGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCg              <  2:1413413/137‑1 (MQ=255)
                                                                                                                                                  tGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCg              >  1:1413413/1‑137 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
CTGTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGG  >  NC_002947/155398‑155690

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