Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,837 C→T intergenic (‑96/+562) yjdN ← / ← yjdM PF06983 family protein YjdN/zinc ribbon domain‑containing protein YjdM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,8370CT98.6% 202.1 / NA 71intergenic (‑96/+562)yjdN/yjdMPF06983 family protein YjdN/zinc ribbon domain‑containing protein YjdM
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (28/41);  minor base G (0/1);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCG  >  NC_000913/4325590‑4326078
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ccGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCTTTTCGCCGAAGTTGAGTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGATCGAAATTTATACAGATGAGAGGCTGCCTGATAAGACGCGCGAGCGTCGCAGCAGGTAGt                                                                                                                                                                                                                                                 <  2:559240/251‑1 (MQ=255)
ccGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGCGCGTCGCATCAGGTAGt                                                                                                                                                                                                                                                 >  1:560875/1‑251 (MQ=255)
   gAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCgg                                                                                                                                                                                                                                              >  1:40914/1‑251 (MQ=255)
     aGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGc                                                                                                                                                                                                                                             >  1:426748/1‑250 (MQ=255)
                    gcTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATgcg                                                                                                                                                                                                                             >  1:480949/1‑251 (MQ=255)
                       gTCCGGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGGCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATtcggcg                                                                                                                                                                                                                          <  2:458846/251‑1 (MQ=255)
                         ccTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTGTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAAGCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGt                                                                                                                                                                                                                         <  2:549398/250‑1 (MQ=255)
                            gcgcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTaa                                                                                                                                                                                                                       <  2:548908/249‑1 (MQ=255)
                             cgcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACg                                                                                                                                                                                                                    <  1:69278/251‑1 (MQ=255)
                                tGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAGATTTATACAGGTGAGAGGCGGCCTGATTAGACGCGCGAGCGTCGCATCAGGTAGTAGGACCTGTTGCCGGATGAGGCGTAAACGCCt                                                                                                                                                                                                                 >  2:196170/1‑251 (MQ=255)
                                 gATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                                                                                 <  1:290827/250‑1 (MQ=255)
                                        gCATGTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                                                                                         <  1:362856/251‑1 (MQ=255)
                                         catgtCGCCGAAGATGATTTTATAGAGCAGTTCCGCTCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACTGTTACCGGCAAAAGAGAGGTAGGGCCTTAACGGCATGATGGTGACCGCAGTTAAGATAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACTCGCGAGCGTCGCATCAGGGAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                                                                                        <  1:398774/247‑1 (MQ=38)
                                         cATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                                                                                        >  1:432277/1‑251 (MQ=255)
                                              cGCCGAAGTTGATGGTATAGAGAAGTTCCGCGCCCAACGTACGGTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTCACGGCATTATGTTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACTGATGAGAGTCGGCCTGATAAGTCGCGTGAGCGTCGCATCAGGTAGTCGGCACTGTTGCTGGATGCGGCGTAAACGCCTAATCCGGCCTACg                                                                                                                                                                                                    <  2:121215/250‑1 (MQ=255)
                                               ggcgaatttGATTTTATAGAGCAGTTCCGCGCCCAACGTACGGTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                                                  <  2:213257/243‑1 (MQ=255)
                                                  gAAGTTGACTTTATAGAGCAGTTCCGCTTCCAACGTACGTTGAGAATAGGCAATCGCGTCGTAACAGTTGCCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGTCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCTAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg                                                                                                                                                                                                 <  1:565312/249‑1 (MQ=38)
                                                             tataGGGCAGTTCCGCGCCCCACGTACGTTGATAATAGGCAAACGCGTCGGAATAGTGACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGGTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACGGAGGGTGCGGGAAtt                                                                                                                                                                                     <  2:97593/250‑1 (MQ=38)
                                                              atatggtaggtctgcgcCCATCGTACGTTGATAATAGGCACTCGCGTCGGAACAGTTTCCGGCAAAAGAGAGGTAGGGACGTGACGTAATGGTGGTGACCTCAGTTAAGAGAATCGAGTTAAGTGTAGTTCGAAATTGATACAGATGAGATGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGTGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                                                                                                                    <  2:145547/237‑1 (MQ=25)
                                                                    cAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGcc                                                                                                                                                                             <  1:500970/251‑1 (MQ=255)
                                                                    cAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGACGGAACAGTTACCGGCAAAAGTGAGGTAGGGACTTAACGGCATGATGGAGACCTCAGTTAAGAAAAGACAGTTAAGTGTAGTTCGCAATTTATACAGTGGAGAGGCGGCCTGATAAAACGCGCGCGCGTCGCCACAGGCAGTCGGCACGGTTTCCGGAAGCGGCGGAAATGCCTTATACGGCCGTCGGAggggtggaggattttagggcc                                                                                                                                                                             >  2:142911/1‑233 (MQ=255)
                                                                         ccGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                                                                         <  1:389465/250‑1 (MQ=255)
                                                                          cgcgCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAACGAGAGGTAGGGACTTAACGGCATGCTGGTGACATCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGGTGCGGCGTAAACGCCTTATCCGGCCTCCGGAGGGAGAGGGAATTTGTAGGCCTGATa                                                                                                                                                                        >  2:473512/1‑250 (MQ=255)
                                                                            cgcCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGCAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATa                                                                                                                                                                        >  1:524644/1‑248 (MQ=255)
                                                                             gcCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGCCTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCCAGCGTCGCATCAGGTAGTCGGCCCTGTTGCCGGATGGGGCGTCAACCCCTTATCCGGCCTACGGCGGGTGCGGGAATTTGTAGGCCTgctacgac                                                                                                                                                                    >  1:475851/1‑244 (MQ=255)
                                                                              cccAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                                                                                   >  1:83517/1‑251 (MQ=255)
                                                                                cacagtacggTGATAAGAGGCTACCGTGTCGGAACAGTGCCCGGCAAAAGAGAGGTAGGTACTTAACGGCATGATGGCGAACTCAGATAAGAGAAGCCATATAAGTGTGGATCGAAATTGATACAGACGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTGAGGGTGCGGGAATGTGTAGGCTCGATACGAcgcg                                                                                                                                                                 <  2:554180/247‑1 (MQ=18)
                                                                                        gtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                                                                                                         <  1:561478/251‑1 (MQ=255)
                                                                                         ttGATAATAGGCAATCGCGTCGGAACAGTTACCGGCCAAAGCGAGGAAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTAATAAGACGCGCGAGCGTCGCCTCAGGTAAGCGACACTGTTGCCGGATGAGGCTAAAACGCCTTATACGGCCTTCGGAGGGTGCGGGGTTTATTAGGCCTGATCAGAAGAGCAAGCGTCg                                                                                                                                                        >  2:502365/1‑251 (MQ=255)
                                                                                          tgagaagaGGCAATCGGGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                                                                                                                        <  2:250178/243‑1 (MQ=255)
                                                                                          tGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                                                                                                                        <  1:104095/250‑1 (MQ=255)
                                                                                              aataGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                                                                                                    <  1:407575/250‑1 (MQ=255)
                                                                                                  ggccctatgtgcGGAACAGTTACCGACAGAATATAGTTAGTGATTTAACGGCATGAGAGTTACCCCATTTAGGAGAAACCAGTTAAGTGTTGTACGAAAGTTAGACAGAGGGGATGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTATTCGGTACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCa                                                                                                                                               <  2:479122/240‑1 (MQ=9)
                                                                                                      agcgcgTCGGAACATTTATCGGCAAAAGAGAGGTAGGGACTGAACGGCATGATGGTGACCTCAGGTAAGAGAAGCCAGTTAAGTGTATTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                                                                                                                                           <  2:451811/249‑1 (MQ=25)
                                                                                                         gcgTCGGAACAGTTAACGGCAAAAGAGAGGGAGGGATTTAACGGAATGATGGTGACCTCAGGTAAGAGAAGCCAGTTAAGTGGAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGc                                                                                                                                         <  2:382544/250‑1 (MQ=21)
                                                                                                               gAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTg                                                                                                                                  <  1:186906/251‑1 (MQ=255)
                                                                                                                            cAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggcgg                                                                                                                     >  1:390004/1‑250 (MQ=255)
                                                                                                                              aaaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGAATAGACGCGCAAGCCTCGCCTCAGGCGGTCGGCACCCTTGCCGGATgcggcgga                                                                                                                    >  2:133234/1‑248 (MQ=25)
                                                                                                                                agagagGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAAc                                                                                                                 <  1:311974/251‑1 (MQ=255)
                                                                                                                                      gTAGGGACGTAACGGAATGAGAGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGGAGTTCGAAGTTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGGCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGccgta                                                                                                           <  2:244626/251‑4 (MQ=18)
                                                                                                                                          ggACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                       >  1:416495/1‑251 (MQ=255)
                                                                                                                                          ggACTGAACGGCAGGATGGTGAACTCAGTGAAGAGAAGCCTGTTTAGGGTTGTACTAAATTTACGAGGATGAGTGGCTGCCTGTTAAGACGCGCTAGCGTTGCATCAGGTAGTCGGCACTGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTATGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                       <  2:290473/251‑1 (MQ=14)
                                                                                                                                                    gCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTAGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCCGTCGGCACCATTGCCGGATGCGGCGTAAAAGCCTTCTCCGGCCtcagg                                                                                               >  2:18239/1‑245 (MQ=255)
                                                                                                                                                            gTGACCTCAGTTAAGAGAAGCCAGTTAAGGGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGATGCGCGAGCGTCGCAGCAGGTAGTCGGTACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCTTACGGAGGGTGCGGGAATTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCggg                                                                                     <  1:471712/251‑1 (MQ=18)
                                                                                                                                                                       tAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGc                                                                         >  1:514345/1‑251 (MQ=255)
                                                                                                                                                                         agagAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                                                                       >  1:519957/1‑251 (MQ=255)
                                                                                                                                                                            gaAGCCAGTTAAGGGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                    <  1:414291/251‑1 (MQ=255)
                                                                                                                                                                             aaGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATa                                                                   <  1:235078/251‑1 (MQ=255)
                                                                                                                                                                               gCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATa                                                                   >  1:424312/1‑249 (MQ=255)
                                                                                                                                                                                     tAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGgcgcgc                                                           >  1:24223/1‑251 (MQ=255)
                                                                                                                                                                                       aGTGTAGTACGAAATTTATACAGAGGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCaa                                                         <  1:434738/251‑1 (MQ=255)
                                                                                                                                                                                         tgtAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                       >  1:482620/1‑251 (MQ=255)
                                                                                                                                                                                            aGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                    <  1:107117/251‑1 (MQ=255)
                                                                                                                                                                                             gTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                   <  1:138380/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaTTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                              <  1:560879/251‑1 (MQ=255)
                                                                                                                                                                                                         taCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGAt                                                                                                                            <  1:329682/167‑1 (MQ=255)
                                                                                                                                                                                                         taCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGGAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCAGCATTGCCGGAt                                                                                                                            >  2:329682/1‑167 (MQ=255)
                                                                                                                                                                                                          aCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                       <  1:133234/250‑1 (MQ=255)
                                                                                                                                                                                                          aCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                                      <  1:473512/251‑1 (MQ=255)
                                                                                                                                                                                                           cAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                <  1:14078/193‑1 (MQ=255)
                                                                                                                                                                                                                   agGCGGCCTGATAAGACGCTCGAGCGTCGAATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTg                             <  1:265037/251‑1 (MQ=255)
                                                                                                                                                                                                                      cggcCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                          >  1:128179/1‑251 (MQ=255)
                                                                                                                                                                                                                      cggcCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGGTGCGGAGTAAACGCCTGATCCGGCCTCGGGGGGGGGGGGGAAATTTGAGGACTGGAAAGACGCGCAAGCGGAGCATCAgagaggcggaaccgtagccg                          >  2:435358/1‑232 (MQ=11)
                                                                                                                                                                                                                          cTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGACGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATg                      >  1:34094/1‑251 (MQ=255)
                                                                                                                                                                                                                            gATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCACCGTTGCCGGATgcg                    <  1:148758/251‑1 (MQ=255)
                                                                                                                                                                                                                                aGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGt                >  1:340803/1‑251 (MQ=255)
                                                                                                                                                                                                                                aGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGACGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGt                >  1:385984/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cgcgAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAAGCCTGCTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACg           >  1:394941/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gggccGCAGCAGGTAGTCGGTACTGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGATGGTGCGGGAAATTGTAGGCCTGATAAGACGCGCAAGCGGTGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAGCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTa      <  2:196575/247‑1 (MQ=11)
                                                                                                                                                                                                                                              cgcctCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGTCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATcc   <  1:293739/246‑1 (MQ=255)
                                                                                                                                                                                                                                              cGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTTCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTCGGCACCGTGGCCGGATGCGGCGTAAACGCCTTATCCg  >  1:456712/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCG  >  NC_000913/4325590‑4326078

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