Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,433:1 +G 5.7% 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 5.7% 478.8 / 22.2 106intergenic (+71/+20)yjbB/pepEputative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base . (53/47);  new base G (3/3);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 107 total aligned reads displayed.

TGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAG  >  NC_000913/4229286‑4229575
                                                                                                                                                    |                                                                                                                                              
tGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGCCAGGCGTCCCATCAGGCACTAGGAACAAGTTGCCGGAT‑c                                                                                                                                               >  2:1102830/1‑149 (MQ=255)
         tGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CAt                                                                                                                                             <  4:456630/142‑1 (MQ=255)
         tGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CAt                                                                                                                                             >  3:456630/1‑142 (MQ=255)
                         ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                                                                                >  1:1108631/1‑139 (MQ=255)
                         ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaaaa                                                                                                                      >  4:190006/1‑149 (MQ=255)
                           aCAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                    <  2:1412576/149‑1 (MQ=255)
                            cAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                             >  2:1454471/1‑139 (MQ=255)
                            cAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                             <  1:1454471/139‑1 (MQ=255)
                             aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                  <  3:190006/149‑1 (MQ=255)
                                   cTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                                                                                            >  1:541436/1‑149 (MQ=255)
                                      gAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCtt                                                                                                         >  4:878426/1‑149 (MQ=255)
                                                               gggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                 >  4:915209/1‑98 (MQ=25)
                                                               gggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                 <  3:915209/100‑3 (MQ=25)
                                                                    gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGATCAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                                                                                            <  3:1245616/100‑1 (MQ=255)
                                                                    gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                                                                                            >  4:1245616/1‑100 (MQ=255)
                                                                         ttAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                             <  4:714846/94‑1 (MQ=255)
                                                                         ttAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                             >  3:714846/1‑94 (MQ=255)
                                                                            atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                                                                                        <  3:291409/112‑1 (MQ=255)
                                                                             tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATGAg                                                                  >  1:732821/1‑149 (MQ=255)
                                                                             tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg                                                                  >  1:1439999/1‑149 (MQ=255)
                                                                                 gtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                              <  2:732821/149‑1 (MQ=255)
                                                                                 gtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                              <  2:1439999/149‑1 (MQ=255)
                                                                                 gtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                              <  1:1102830/149‑1 (MQ=255)
                                                                                  tgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggtg                                                             >  2:77183/1‑149 (MQ=255)
                                                                                     aTCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTc                                                          >  2:800178/1‑149 (MQ=255)
                                                                                         tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  1:382893/129‑1 (MQ=255)
                                                                                         tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  2:382893/1‑129 (MQ=255)
                                                                                           tAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACgg                                                                                                                   >  1:824668/1‑86 (MQ=255)
                                                                                           tAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACgg                                                                                                                   <  2:824668/86‑1 (MQ=255)
                                                                                             ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  3:459137/125‑1 (MQ=255)
                                                                                             ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  4:459137/1‑125 (MQ=255)
                                                                                             ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca                                                                      <  3:1145975/129‑1 (MQ=255)
                                                                                             ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca                                                                      >  4:1145975/1‑129 (MQ=255)
                                                                                             ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAcc                                                  >  3:386122/1‑149 (MQ=255)
                                                                                               ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                    >  2:981193/1‑81 (MQ=255)
                                                                                               ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                    <  1:981193/81‑1 (MQ=255)
                                                                                               ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCa                                                >  4:254119/1‑149 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  3:609348/122‑1 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  1:1131782/122‑1 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  4:609348/1‑122 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  2:1131782/1‑122 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                              <  1:669408/134‑1 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                              >  2:669408/1‑134 (MQ=255)
                                                                                                cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                               >  2:289475/1‑149 (MQ=255)
                                                                                                 tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                         >  1:463968/1‑106 (MQ=255)
                                                                                                 tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                         <  2:463968/106‑1 (MQ=255)
                                                                                                  gATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                             >  4:1194746/1‑149 (MQ=255)
                                                                                                    tAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                           <  1:77183/149‑1 (MQ=255)
                                                                                                     aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                          >  4:316671/1‑101 (MQ=255)
                                                                                                     aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                          <  3:316671/101‑1 (MQ=255)
                                                                                                     aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtag                                          <  4:386122/149‑1 (MQ=255)
                                                                                                       gACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagcg                                        >  4:1466812/1‑149 (MQ=255)
                                                                                                        aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  3:1049191/1‑114 (MQ=255)
                                                                                                        aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  4:1049191/114‑1 (MQ=255)
                                                                                                        aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                <  4:1321338/124‑1 (MQ=255)
                                                                                                        aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CAAGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                >  3:1321338/1‑124 (MQ=255)
                                                                                                           cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    <  2:541436/149‑1 (MQ=255)
                                                                                                           cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    <  1:289475/149‑1 (MQ=255)
                                                                                                           cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    <  3:254119/149‑1 (MQ=255)
                                                                                                                aGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTa                               <  4:1345889/149‑1 (MQ=255)
                                                                                                                 gCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTAc                              >  4:745074/1‑149 (MQ=255)
                                                                                                                   gTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTc                            >  4:301450/1‑149 (MQ=255)
                                                                                                                     cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          >  4:302076/1‑149 (MQ=255)
                                                                                                                     cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          >  4:487231/1‑149 (MQ=255)
                                                                                                                       cATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                        <  1:800178/149‑1 (MQ=255)
                                                                                                                        aTCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtgt                       >  3:1399664/1‑149 (MQ=255)
                                                                                                                           aGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGAt                    >  3:1243108/1‑149 (MQ=255)
                                                                                                                            ggCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    <  2:747551/132‑1 (MQ=255)
                                                                                                                            ggCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    >  1:747551/1‑132 (MQ=255)
                                                                                                                             gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                     >  2:611851/1‑146 (MQ=255)
                                                                                                                             gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                     <  1:611851/146‑1 (MQ=255)
                                                                                                                              cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                 <  3:487231/149‑1 (MQ=255)
                                                                                                                               aCTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  2:1180971/1‑91 (MQ=255)
                                                                                                                               aCTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  1:1180971/91‑1 (MQ=255)
                                                                                                                                 tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt              >  3:872550/1‑149 (MQ=255)
                                                                                                                                 tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt              <  4:872550/149‑1 (MQ=255)
                                                                                                                                 tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt              <  3:1194746/149‑1 (MQ=255)
                                                                                                                                  aTGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt               <  2:1001903/147‑1 (MQ=255)
                                                                                                                                  aTGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt               >  1:1001903/1‑147 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    <  1:1143550/123‑1 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                    >  2:1143550/1‑123 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          <  1:1088448/133‑1 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          >  4:622161/1‑133 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          >  2:1088448/1‑133 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          <  3:622161/133‑1 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCCTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                              <  3:1047183/113‑1 (MQ=255)
                                                                                                                                     aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCCTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                              >  4:1047183/1‑113 (MQ=255)
                                                                                                                                      aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCtt         >  4:1320296/1‑149 (MQ=255)
                                                                                                                                       cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          <  1:422202/83‑1 (MQ=255)
                                                                                                                                       cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                          >  2:422202/1‑83 (MQ=255)
                                                                                                                                       cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTc        <  4:1243108/149‑1 (MQ=255)
                                                                                                                                        aaGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTcc       >  1:151725/1‑149 (MQ=255)
                                                                                                                                          attGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagcgt                                       <  3:952547/113‑1 (MQ=255)
                                                                                                                                           ttGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc           >  1:951228/2‑142 (MQ=255)
                                                                                                                                           ttGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc           <  2:951228/141‑1 (MQ=255)
                                                                                                                                           aTGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt    <  3:1466812/149‑1 (MQ=255)
                                                                                                                                           aTGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt    >  2:1204540/1‑149 (MQ=255)
                                                                                                                                             gCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          >  1:661130/1‑125 (MQ=255)
                                                                                                                                             gCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAg  >  1:861453/1‑149 (MQ=255)
                                                                                                                                                ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                          >  2:555687/1‑122 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAG  >  NC_000913/4229286‑4229575

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