Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 155,546 A→C 10.4% 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,5460AC10.4% 220.2 / 15.0 77intergenic (‑251/+32)yjbB/ygdQNa+/Pi‑cotransporter/transport protein
Reads supporting (aligned to +/- strand):  ref base A (30/39);  new base C (3/5);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGCC  >  NC_002947/155401‑155692
                                                                                                                                                   |                                                                                                                                                  
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGTAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGaa                                                                                                                                                   <  5:95322/147‑1 (MQ=255)
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGaa                                                                                                                                                   <  5:317529/147‑1 (MQ=255)
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGaa                                                                                                                                                   <  3:71177/147‑1 (MQ=255)
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGaa                                                                                                                                                   <  3:253195/147‑1 (MQ=255)
   ccAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGcc                                                                                                                                              >  5:226873/1‑149 (MQ=255)
    cAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGccc                                                                                                                                             <  5:251781/149‑1 (MQ=255)
      gCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCtt                                                                                                                                           <  3:135920/149‑1 (MQ=255)
       cTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttg                                                                                                                                          <  8:176913/149‑1 (MQ=255)
       cTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttg                                                                                                                                          <  6:226873/149‑1 (MQ=255)
       cTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttg                                                                                                                                          >  3:215048/1‑149 (MQ=255)
        ttCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGTCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgt                                                                                                                                         >  4:105630/1‑149 (MQ=255)
        ttCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtt                                                                                                                                           >  1:15270/1‑147 (MQ=255)
        ttCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtt                                                                                                                                           <  2:15270/147‑1 (MQ=255)
        ttCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgt                                                                                                                                         >  8:349811/1‑149 (MQ=255)
          cgcgTGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCttttt                                                                                                                                        <  6:257673/148‑1 (MQ=255)
          cGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGTAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCGCttgtt                                                                                                                                        <  2:129034/148‑1 (MQ=255)
          cGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                        <  4:215048/148‑1 (MQ=255)
           gCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgttg                                                                                                                                       <  1:24205/148‑1 (MQ=255)
            cATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCAGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtttttaa                                                                                                                                      <  1:298618/149‑3 (MQ=255)
            cATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTtgtg                                                                                                                                     >  6:61584/1‑149 (MQ=255)
             aTGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                    <  3:168067/149‑1 (MQ=255)
             aTGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                                    <  7:50846/149‑1 (MQ=255)
                gCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAACAGGGCCCTGCGAAGGGCCCTTGTTGtgctg                                                                                                                                  <  7:269878/148‑1 (MQ=255)
                gCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACCGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgctg                                                                                                                                  <  3:320207/148‑1 (MQ=255)
                       tttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCt                                                                                                                          >  5:104307/1‑149 (MQ=255)
                           tttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTgg                                                                                                                       >  1:73460/1‑149 (MQ=255)
                           tttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtttttaat                                                                                                                                     >  3:175507/1‑131 (MQ=255)
                           tttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtttttaat                                                                                                                                     <  4:175507/134‑4 (MQ=255)
                            ttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGcc                                                                                                                     >  5:108127/1‑149 (MQ=255)
                             tGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCt                                                                                                                    >  7:253074/1‑149 (MQ=255)
                              gAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCt                                                                                                                    >  7:129777/1‑148 (MQ=255)
                               aaaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTa                                                                                                                  >  7:272122/1‑149 (MQ=255)
                                aaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                                 <  3:363737/149‑1 (MQ=255)
                                 aCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtttttaatgcc                                                                                                                                  >  6:49440/1‑125 (MQ=255)
                                 aCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtttttaatgcc                                                                                                                                  <  5:49440/131‑7 (MQ=255)
                                   aGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGa                                                                                                              <  8:253074/149‑1 (MQ=255)
                                       cTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGc                                                                                                                      <  3:314543/137‑1 (MQ=255)
                                       cTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGc                                                                                                                      >  4:314543/1‑137 (MQ=255)
                                         tGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCgg                                                                                                         >  2:287025/1‑148 (MQ=255)
                                         tGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCgg                                                                                                         >  5:241619/1‑148 (MQ=255)
                                              cGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTc                                                                                                   >  8:332312/1‑149 (MQ=255)
                                                 aTGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCtt                                                                                                <  6:241619/149‑1 (MQ=255)
                                                 aTGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCtt                                                                                                <  6:108127/149‑1 (MQ=255)
                                                    aCTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGTAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTgc                                                                                              <  4:232998/148‑1 (MQ=255)
                                                             gggggCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                                    <  4:119421/149‑1 (MQ=255)
                                                              ggggCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCGTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAcc                                                                                   <  2:358727/149‑1 (MQ=255)
                                                                ggCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                                    >  5:247488/1‑146 (MQ=255)
                                                                ggCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                                    <  6:247488/146‑1 (MQ=255)
                                                                   tgctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGtt                                                                               >  5:3363/1‑148 (MQ=255)
                                                                     gctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttg                                                                                                                                          <  2:274639/88‑1 (MQ=255)
                                                                      ctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCa                                                                             <  6:3363/148‑1 (MQ=255)
                                                                      ctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCAt                                                                            >  8:182121/1‑149 (MQ=255)
                                                                        tCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGAGCAGTTTCACCGGTTCATGa                                                                          >  3:205149/1‑149 (MQ=255)
                                                                              gCAAGTCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTg                                                                    <  3:105630/149‑1 (MQ=255)
                                                                              gCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgttg                                                                                                                                       <  2:372359/82‑1 (MQ=255)
                                                                              gCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTg                                                                    <  7:332312/149‑1 (MQ=255)
                                                                                         caGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGTCCCGGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                         <  8:293742/149‑1 (MQ=255)
                                                                                         caGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                         >  5:297583/1‑149 (MQ=255)
                                                                                          tctGACAACGCGCAGGAAAAAACGCACATAAAAAAGGCTTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTTACCGGTATTGTCCTTGTGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                         <  1:259612/146‑1 (MQ=255)
                                                                                          aGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCgg                                                        <  5:61584/149‑1 (MQ=255)
                                                                                              aCAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                         >  2:259612/1‑144 (MQ=255)
                                                                                                        aGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                           >  5:357193/1‑149 (MQ=255)
                                                                                                         ggAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                           <  1:183688/148‑1 (MQ=255)
                                                                                                         ggAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                           >  2:183688/1‑148 (MQ=255)
                                                                                                            aaaaaCGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACACGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCAt                                       >  1:155965/1‑149 (MQ=255)
                                                                                                               aaCGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTACATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCg                                    >  2:60664/1‑149 (MQ=255)
                                                                                                                    acaGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGt                               >  7:133666/1‑149 (MQ=255)
                                                                                                                      aGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGCTCCTTGCCCTCCTTGCGGGCCCTGGCGGTGCGCAGGCGGATGtt                              >  4:290324/1‑148 (MQ=255)
                                                                                                                        aaaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGAt                           <  2:171126/149‑1 (MQ=255)
                                                                                                                              ggCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTc                     <  6:297583/149‑1 (MQ=255)
                                                                                                                               gCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTcc                    <  2:73460/149‑1 (MQ=255)
                                                                                                                                   tAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTCCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCg                >  3:176502/1‑149 (MQ=255)
                                                                                                                                      aaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGtt                              >  6:376467/1‑132 (MQ=255)
                                                                                                                                      aaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGtt                              <  5:376467/132‑1 (MQ=255)
                                                                                                                                      aaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTCGCGGGCCATGGCGGTGCGCATGCGGATGtt                              <  5:267852/132‑1 (MQ=255)
                                                                                                                                      aaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTAGCGGGCCATGGCGGTGCGCATGCGGATGtt                              >  6:267852/1‑132 (MQ=255)
                                                                                                                                       cagtgCCCTGGGAAGGGGCCTTGTTGTGGTGCACTTCCTTGGAATTATTGACCCGTAATGTCCTTGCGCAGTATCACCGGTTCATGATCGTTGCCATCCTTGCGGGCCATTCCGGTGCGCATGCGGATTTTGATGGCTTCCACCGCCAg            <  4:27827/145‑1 (MQ=255)
                                                                                                                                        aGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGc           <  6:357193/149‑1 (MQ=255)
                                                                                                                                                 cGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAgg    <  7:182121/147‑1 (MQ=255)
                                                                                                                                                 cGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAgg    <  6:371985/147‑1 (MQ=255)
                                                                                                                                                 cGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGcc  <  6:104307/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGCC  >  NC_002947/155401‑155692

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