Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,433:1 +G 6.8% intergenic (+71/+20) yjbB → / ← pepE putative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,229,4331.G 6.8% 509.6 / 30.4 118intergenic (+71/+20)yjbB/pepEputative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base . (46/64);  new base G (2/6);  total (48/70)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-03

Only 100 of 125 total aligned reads displayed.

CTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACC  >  NC_000913/4229285‑4229578
                                                                                                                                                     |                                                                                                                                                 
cTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATGAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAt                                                                                                                                                    >  1:633410/1‑149 (MQ=255)
   aaCTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTTATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CAt                                                                                                                                                <  4:578508/149‑1 (MQ=255)
    aCTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAtgcat                                                                                                                                                <  4:143056/149‑5 (MQ=255)
     cTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGGCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATga                                                                                                                                              >  3:79380/1‑148 (MQ=255)
          tGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGCAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCg                                                                                                                                         <  3:727098/149‑1 (MQ=255)
            ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGc                                                                                                                                       <  1:727405/149‑1 (MQ=255)
            ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGc                                                                                                                                       >  1:76465/1‑149 (MQ=255)
            ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGc                                                                                                                                       >  2:770757/1‑149 (MQ=255)
            ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTAGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGc                                                                                                                                       >  2:870668/1‑149 (MQ=255)
            ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCCTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGc                                                                                                                                       <  2:766808/149‑1 (MQ=255)
              ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCat                                                                                                                                     <  3:103551/149‑2 (MQ=255)
               tGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                    <  4:468693/149‑3 (MQ=255)
                          ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaaaa                                                                                                                         >  2:149089/1‑149 (MQ=255)
                              aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                     <  2:76465/149‑1 (MQ=255)
                              aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                     <  2:675039/149‑1 (MQ=255)
                              aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                     <  2:577296/149‑1 (MQ=255)
                              aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                     <  4:389277/149‑1 (MQ=255)
                                       gAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCtt                                                                                                            >  3:840558/1‑149 (MQ=255)
                                            gCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                                                                                       <  2:636613/149‑1 (MQ=255)
                                                    aagatgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATGAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGTAACCGc                                                                                               <  2:633410/149‑1 (MQ=255)
                                                       atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                       <  1:10185/122‑1 (MQ=255)
                                                       atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                       >  2:10185/1‑122 (MQ=255)
                                                       atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                            <  1:148785/149‑1 (MQ=255)
                                                       atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                            <  1:184704/149‑1 (MQ=255)
                                                       atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                            <  3:521026/149‑1 (MQ=255)
                                                       atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                            <  1:604007/149‑1 (MQ=255)
                                                        tgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                                                                                           >  1:859408/1‑149 (MQ=255)
                                                         gaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttctt                                                                                          >  1:236466/1‑149 (MQ=255)
                                                           agGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCa                                                                                        <  2:222413/149‑1 (MQ=255)
                                                               cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAg                                                                                    <  1:149089/149‑1 (MQ=255)
                                                                 ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc                                                                                  <  4:322985/149‑1 (MQ=255)
                                                                 ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc                                                                                  <  3:102662/149‑1 (MQ=255)
                                                                   actaGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACACGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCtt                                                                                <  1:378085/146‑1 (MQ=255)
                                                                     gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAt                                                                                                                                                    >  3:452110/1‑80 (MQ=255)
                                                                     gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAt                                                                                                                                                    <  4:452110/80‑1 (MQ=255)
                                                                     gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaa                                                                              >  2:249272/1‑149 (MQ=255)
                                                                      attATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                             <  3:188350/147‑1 (MQ=255)
                                                                      aGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                             <  1:686210/149‑1 (MQ=255)
                                                                      aGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                             <  4:517189/149‑1 (MQ=255)
                                                                        tATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAaca                                                                           >  1:180912/1‑149 (MQ=255)
                                                                         aTTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGTATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                          >  4:176857/1‑149 (MQ=255)
                                                                         aTTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                          >  1:61604/1‑149 (MQ=255)
                                                                         aTTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                          >  3:195371/1‑149 (MQ=255)
                                                                          ttAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca                                                                         <  2:236466/149‑1 (MQ=255)
                                                                             atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCa                                                                                     <  3:134640/134‑1 (MQ=255)
                                                                             atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCa                                                                                     >  4:134640/1‑134 (MQ=255)
                                                                             atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATaa                                                                      >  3:634734/1‑149 (MQ=255)
                                                                              tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                    <  3:521141/86‑3 (MQ=21)
                                                                              tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                    >  4:521141/1‑84 (MQ=21)
                                                                               atTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgt                                                                    <  3:428438/149‑1 (MQ=255)
                                                                                ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                   >  2:437612/1‑149 (MQ=255)
                                                                                ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                   >  2:895437/1‑149 (MQ=255)
                                                                                ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                   >  1:198414/1‑149 (MQ=255)
                                                                                ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                   <  2:450209/149‑1 (MQ=255)
                                                                                   tgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggtg                                                                <  1:536680/149‑1 (MQ=255)
                                                                                     aaTCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCgg                                                                                                     >  1:242131/1‑110 (MQ=255)
                                                                                       tCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCg                                                            <  3:464646/149‑1 (MQ=255)
                                                                                        ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCgg                                                           <  2:859408/149‑1 (MQ=255)
                                                                                         cTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                          <  2:603661/149‑1 (MQ=255)
                                                                                          tgtgGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGc                                                         >  1:771569/1‑149 (MQ=255)
                                                                                          tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGct                                                                                              <  1:834876/112‑1 (MQ=255)
                                                                                          tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGct                                                                                              >  2:834876/1‑112 (MQ=255)
                                                                                          tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt                                                                                 <  1:388459/125‑1 (MQ=255)
                                                                                               gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAccc                                                    <  1:249272/149‑1 (MQ=255)
                                                                                                 cTGATAAGACGCGTCAAGCGTCGCATCAGGCAGTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                                  >  4:462866/1‑149 (MQ=255)
                                                                                                     tAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                              >  4:364437/1‑149 (MQ=255)
                                                                                                     tAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                              >  2:239405/1‑149 (MQ=255)
                                                                                                      aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtag                                             >  2:739979/1‑149 (MQ=255)
                                                                                                      aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtag                                             >  1:590781/1‑149 (MQ=255)
                                                                                                       aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                            <  2:136434/149‑1 (MQ=255)
                                                                                                        gACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagcg                                           >  2:207822/1‑149 (MQ=255)
                                                                                                         aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCgg                                                           >  4:662505/1‑132 (MQ=255)
                                                                                                         aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCgg                                                           <  3:662505/132‑1 (MQ=255)
                                                                                                         aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                 >  2:35444/1‑142 (MQ=255)
                                                                                                          cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTg                                         >  3:156885/1‑149 (MQ=255)
                                                                                                              tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCt                                     >  2:291506/1‑149 (MQ=255)
                                                                                                                 aGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTa                                  <  4:669065/149‑1 (MQ=255)
                                                                                                                      cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                             <  2:185573/149‑1 (MQ=255)
                                                                                                                      cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                             <  4:446753/149‑1 (MQ=255)
                                                                                                                      cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                             <  3:835764/149‑1 (MQ=255)
                                                                                                                       gCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACt                            >  4:713344/1‑149 (MQ=255)
                                                                                                                           cAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                            >  1:183300/1‑129 (MQ=255)
                                                                                                                           cAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                            <  2:183300/129‑1 (MQ=255)
                                                                                                                           cAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                        <  3:364437/149‑1 (MQ=255)
                                                                                                                           cAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                        <  1:239405/149‑1 (MQ=255)
                                                                                                                              gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGAccc                     >  1:494524/1‑149 (MQ=255)
                                                                                                                              gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                     >  4:365963/1‑149 (MQ=255)
                                                                                                                              gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGATGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCATCGTCGCGTGACCTTTACTCACTGTGATcc                     >  2:330358/1‑149 (MQ=255)
                                                                                                                               cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                    <  4:79380/149‑1 (MQ=255)
                                                                                                                               cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                    <  3:749638/149‑1 (MQ=255)
                                                                                                                               cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                    <  3:491915/149‑1 (MQ=255)
                                                                                                                                 cTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                  <  3:429508/149‑1 (MQ=255)
                                                                                                                                      aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                      <  4:400639/140‑1 (MQ=255)
                                                                                                                                      aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                      >  3:58020/1‑140 (MQ=255)
                                                                                                                                      aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                      >  3:400639/1‑140 (MQ=255)
                                                                                                                                      aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                      <  4:58020/140‑1 (MQ=255)
                                                                                                                                         aaGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTcc          >  2:636515/1‑149 (MQ=255)
                                                                                                                                             tGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTa      <  3:365963/149‑1 (MQ=255)
                                                                                                                                               ccGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGa    >  1:861267/1‑149 (MQ=255)
                                                                                                                                                 ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAcc  >  1:143945/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACC  >  NC_000913/4229285‑4229578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑