Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 155,543 G→T 5.3% intergenic (‑248/+35) 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,5430GT 5.3% 247.9 / 3.4 76intergenic (‑248/+35)yjbB/ygdQNa+/Pi‑cotransporter/transport protein
Reads supporting (aligned to +/- strand):  ref base G (30/42);  new base T (2/2);  total (32/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 = 2.24e-01

TCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAG  >  NC_002947/155401‑155686
                                                                                                                                                |                                                                                                                                               
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACTCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAgg                                                                                                                                           <  1:232167/149‑1 (MQ=255)
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAg                                                                                                                                            <  3:219620/148‑1 (MQ=255)
tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCGGCTTCGCGGGCTAGCCACCACAGTG‑‑ACAACGCGCATGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAg                                                                                                                                            <  5:139007/148‑1 (MQ=255)
  gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGc                                                                                                                                         >  3:223928/1‑149 (MQ=255)
  gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGc                                                                                                                                         >  7:263040/1‑149 (MQ=255)
  gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCACGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGc                                                                                                                                         >  6:111580/1‑149 (MQ=255)
    cAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGccc                                                                                                                                       <  6:251345/149‑1 (MQ=255)
    cAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGACTGGccc                                                                                                                                       >  4:253645/1‑149 (MQ=255)
     aGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGTCTCAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGcc                                                                                                                                        <  1:105934/149‑1 (MQ=255)
     aGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGAGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCt                                                                                                                                      <  7:329141/149‑1 (MQ=255)
      gCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACCCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCtt                                                                                                                                     <  7:272289/149‑1 (MQ=255)
          cGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGACAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgt                                                                                                                                   <  8:202314/147‑1 (MQ=255)
          cGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                  <  7:207878/148‑1 (MQ=255)
          cGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                  <  5:323046/148‑1 (MQ=255)
          cGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                  <  5:188260/148‑1 (MQ=255)
            cATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTtgtg                                                                                                                               <  7:168132/149‑1 (MQ=255)
             aTGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACTCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                              <  8:333109/149‑1 (MQ=255)
                    tCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACtt                                                                                                                       <  5:323787/149‑1 (MQ=255)
                      atttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAggg                                                                                                                                          <  3:289521/128‑1 (MQ=255)
                      atttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAggg                                                                                                                                          >  4:289521/1‑128 (MQ=255)
                        ttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                   >  4:154762/1‑149 (MQ=255)
                        ttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                   >  5:196689/1‑149 (MQ=255)
                          atttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCACGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                  <  5:111580/148‑1 (MQ=255)
                           tttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGc                                                                                                                >  5:227600/1‑149 (MQ=255)
                               aaaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAACGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTa                                                                                                            <  6:262537/149‑1 (MQ=255)
                                aaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTa                                                                                                            <  8:274835/148‑1 (MQ=255)
                                aaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTa                                                                                                            <  6:196689/148‑1 (MQ=255)
                                aaCAGATCTTGAAAAGAATGACTGCCCGTGGGGTCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCAGAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                           <  3:253645/149‑1 (MQ=255)
                                 aCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGcc                                                                                                               >  4:220688/1‑144 (MQ=255)
                                 aCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGcc                                                                                                               <  3:220688/144‑1 (MQ=255)
                                 aCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAtt                                                                                                          <  4:223928/149‑1 (MQ=255)
                                        ttGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCATGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCg                                                                                                    >  4:100587/1‑148 (MQ=255)
                                        ttGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCAACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCg                                                                                                    >  5:176239/1‑148 (MQ=255)
                                        tagAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTTAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTCGCCTTATTGAccct                                                                                                    >  8:265443/3‑147 (MQ=255)
                                          gAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGAccc                                                                                                     <  8:27283/145‑1 (MQ=255)
                                           aaaCGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAAt                                                                                                >  3:348911/1‑149 (MQ=255)
                                            aaCGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAAt                                                                                                <  8:176438/148‑1 (MQ=255)
                                            aaCGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAAt                                                                                                <  3:147090/148‑1 (MQ=255)
                                                 ttctcTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCtt                                                                                          <  7:201012/145‑1 (MQ=255)
                                                 aTGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCATGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCtt                                                                                          <  3:100587/149‑1 (MQ=255)
                                                    aCTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTgcg                                                                                       >  1:348120/1‑149 (MQ=255)
                                                    aCTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTgc                                                                                        <  6:210217/148‑1 (MQ=255)
                                                      tGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGAGCAGGAAAAAACGCACAGAAAAAAGGAATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCAATTCATTGGCCTTATTGACCCGGATTGTCATTgagaa                                                                                     >  8:201012/1‑145 (MQ=255)
                                                           gTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTc                                                                                >  2:275772/1‑149 (MQ=255)
                                                             gggggCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                              >  3:215468/1‑149 (MQ=255)
                                                              gtggCTGCATCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGCAAAACACGCACAGAAAAAAGGCCTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                              <  8:263040/146‑1 (MQ=255)
                                                                 gctgctTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGt                                                                          <  6:227600/149‑1 (MQ=255)
                                                                   tgctTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGtt                                                                         >  8:189973/1‑148 (MQ=255)
                                                                     ctTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTc                                                                        <  2:348120/147‑1 (MQ=255)
                                                                     ctTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCAt                                                                      <  5:167746/149‑1 (MQ=255)
                                                                      ttCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCAt                                                                      >  7:67449/1‑148 (MQ=255)
                                                                        cgcgGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGa                                                                    <  7:285094/148‑1 (MQ=255)
                                                                         gcgGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATc                                                                  <  4:71535/149‑1 (MQ=255)
                                                                             gCAAGCCACCACAGGG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCt                                                                <  4:234461/147‑1 (MQ=255)
                                                                                      cacaGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTg                                                     >  1:164730/1‑149 (MQ=255)
                                                                                        caGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCATTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                   <  6:176239/149‑1 (MQ=255)
                                                                                                 aaCGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATg                                            >  6:120333/1‑149 (MQ=255)
                                                                                                         ggAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgcg                                    <  4:348911/149‑1 (MQ=255)
                                                                                                              aaaCGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGc                               <  5:120333/149‑1 (MQ=255)
                                                                                                                  gCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGAt                           >  6:60922/1‑149 (MQ=255)
                                                                                                                        aaaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGAt                     >  8:204408/1‑149 (MQ=255)
                                                                                                                         aaaaaGGCATTAAAAAGGGCCCTGCGAAAGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATg                    <  1:125897/149‑1 (MQ=255)
                                                                                                                            aaGGCATTAAAAACGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCTGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCt                 <  4:101845/149‑1 (MQ=255)
                                                                                                                                 aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCAc            >  5:22163/1‑149 (MQ=255)
                                                                                                                                  ttAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCAcc           >  1:354628/1‑149 (MQ=255)
                                                                                                                                    aaaaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCTATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCg          <  1:275772/148‑1 (MQ=255)
                                                                                                                                    aaaaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACAGGTTCATGATCCTTGCCATGCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGc         >  3:195243/1‑149 (MQ=255)
                                                                                                                                      caaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGc         >  2:48137/2‑147 (MQ=255)
                                                                                                                                      caaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCAAGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGc         >  6:282058/2‑147 (MQ=255)
                                                                                                                                       aaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAg      >  4:116158/1‑149 (MQ=255)
                                                                                                                                       aaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGACATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAg      >  3:31769/1‑149 (MQ=255)
                                                                                                                                        aGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGTTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGc     <  7:265443/149‑1 (MQ=255)
                                                                                                                                          ggCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCga   <  7:189973/149‑1 (MQ=255)
                                                                                                                                             ccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATg                    >  6:223403/1‑129 (MQ=255)
                                                                                                                                             ccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATg                    <  5:223403/129‑1 (MQ=255)
                                                                                                                                             ccTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCgag  >  3:241720/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTG‑‑ACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAG  >  NC_002947/155401‑155686

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: