Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 8.1% intergenic (+19/+21) metE → / ← ysgA 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3330AT 8.1% 399.7 / 19.1 145intergenic (+19/+21)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (69/64);  new base T (6/6);  total (75/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 149 total aligned reads displayed.

CGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACAC  >  NC_000913/4015184‑4015482
                                                                                                                                                     |                                                                                                                                                     
   cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                    >  1:1148769/1‑150 (MQ=255)
     aaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacca                                                                                                                                                  >  2:1978079/1‑150 (MQ=255)
        cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCAccc                                                                                                                                               <  1:1295035/150‑1 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                                                                             <  1:2125083/150‑1 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                                                                             <  1:885859/150‑1 (MQ=255)
          cATTCCGGCACAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                                                                             >  1:1199027/1‑150 (MQ=255)
           aTTCCGGCAGAGCGCCTGTGTGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCAAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                                                                            <  2:1657647/150‑1 (MQ=255)
           aTTCCGGCAGAGCGCCTGTGGGTCCACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCTCTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                                                                            <  2:1069666/150‑1 (MQ=255)
           aTTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGt                                                                                                                                            <  1:1401885/150‑1 (MQ=255)
           aTTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGt                                                                                                                                            >  2:1145936/1‑150 (MQ=255)
            ttCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTc                                                                                                                                           >  1:486205/1‑150 (MQ=255)
            ttCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTc                                                                                                                                           >  1:986623/1‑150 (MQ=255)
              ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttg                                                                                                                                         <  2:1705996/150‑4 (MQ=255)
                ggCAGAGCGCCTGTGGGTCAACCCGGACTGTGGTCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCtttt                                                                                                                                       >  1:1509212/1‑150 (MQ=255)
                   cgagCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                                                                    <  2:460837/149‑1 (MQ=255)
                    gagCGCCTGGGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                                                                   <  2:345458/150‑1 (MQ=255)
                        gcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTAc                                                                                                                               <  2:1240836/150‑1 (MQ=255)
                          gtgtgGGTCAACCCGGACTGTGGCCGGAAACCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCATAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAg                                                                                                                             <  2:1859927/149‑1 (MQ=255)
                          cTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAg                                                                                                                             >  1:717526/1‑150 (MQ=255)
                              gggTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCTGAAACCCGCGCGGCACTGGCGCACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGCGTAAAATCCAAACCGGGTGGTAATACACCCCGGTCTTTTCTCATTACAGCGAc                                                                                                                         >  2:416388/1‑150 (MQ=255)
                                gTCACCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGCctt                                                                                                                       >  1:822977/1‑150 (MQ=255)
                                gTCAACCCGGACTGTGGCCTTAAAACGCGCGGCTGGCCAGATACACGCGCGGCGCTGGCGAACATGGTGCTGGCGGCGCAGAACTTGCGTCGGGTGTACAATCCAAACCGGGTGGTAATAGCAACCGGTCGTATGTCATTACAGCGActt                                                                                                                       >  2:410906/1‑150 (MQ=255)
                                gTCAACCCGGACTGTGGCCTAAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActt                                                                                                                       >  1:447815/1‑150 (MQ=255)
                                    aCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCGGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACttcttc                                                                                                                   <  2:2171981/150‑1 (MQ=255)
                                    aCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACttcttc                                                                                                                   >  2:921976/1‑150 (MQ=255)
                                      ccGGTCGGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTccc                                                                                                                 <  1:863799/150‑1 (MQ=255)
                                         gACTGTGGCCTGAAAACGCGCGGCTGGCCAGCAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATTCCACCCGGTCTTTTCTCATTACAGCGACTTCTTccaccc                                                                                                              >  1:361733/1‑146 (MQ=255)
                                             gtgGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAc                                                                                                          >  2:2215617/1‑150 (MQ=255)
                                             gtgGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAc                                                                                                          <  1:1156750/150‑1 (MQ=255)
                                               ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                        >  1:590434/1‑150 (MQ=255)
                                               ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACAACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                        >  2:643619/1‑150 (MQ=255)
                                                   tGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGGCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTTCCGGGCCTTCTGCCCACCAGACTGCtga                                                                                                    >  1:417295/1‑148 (MQ=255)
                                                     aaaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                  >  2:404447/1‑150 (MQ=255)
                                                      aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAc                                                                                                 >  1:1672372/1‑150 (MQ=255)
                                                       aaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                                                                                >  2:323437/1‑150 (MQ=255)
                                                        aCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCa                                                                                               <  1:1393888/150‑1 (MQ=255)
                                                         cgcgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                                                                              >  2:930909/1‑150 (MQ=255)
                                                            gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTcc                                                                                           <  2:986623/150‑1 (MQ=255)
                                                                  ggCCAGAAACCCGCGCGGCACTGGCTAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATa                                                                                     <  1:1606440/150‑1 (MQ=255)
                                                                  ggCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATa                                                                                     >  2:1570966/1‑150 (MQ=255)
                                                                    ccAGAAAACCGCGCGGCACTGGAGAACATGGTGCAGTCGGCGAAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATAACAACCGGTCTTATATCATTACAGCTACTTATTCACCACATACTGCTTAAACAATTACATCATACg                                                                                   >  2:582942/1‑150 (MQ=255)
                                                                         aaCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTTGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgcc                                                                              <  2:831272/150‑1 (MQ=255)
                                                                           cccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                            >  2:1932580/1‑150 (MQ=255)
                                                                               cgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCAAGAATACGCTgccagccaa                                                                        >  2:696808/1‑149 (MQ=255)
                                                                                 cgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACCGCGACTTCTTCCCACCCTACTGCTGAAACCATTCCAGCATTCGCTGCCAGCCAact                                                                      >  1:229281/1‑150 (MQ=255)
                                                                                 cgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTACAGCATACGCTGCCATCCAtct                                                                      >  2:573161/1‑150 (MQ=255)
                                                                                   gCACTGTCGAACATGGTTCAGGCGTGGCAGAATTTGCGTCGGGGGTAAAATCCAAACCGGGGTGTAATACCACCCGTTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                                                                    <  2:921053/150‑1 (MQ=255)
                                                                                     aCTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTGAACCATTCCAGCATACGCTGCCAGCCATCTTCTg                                                                  <  2:426991/150‑1 (MQ=255)
                                                                                     aCTGGCGAACATGGGGCATGCGGCGCAGAACTTGTGTCGGGGGTAAACTACAAACCGGGTGGAAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg                                                                  <  1:213844/150‑1 (MQ=255)
                                                                                      cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc                                                                 <  1:1663503/150‑1 (MQ=255)
                                                                                        ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg                                                               <  2:627674/150‑1 (MQ=255)
                                                                                        ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg                                                               >  2:712312/1‑150 (MQ=255)
                                                                                        ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg                                                               >  2:856323/1‑150 (MQ=255)
                                                                                            aaCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGc                                                           <  2:1337128/150‑1 (MQ=255)
                                                                                             aCATGGTGCAGGCGGCGCAGAACTGGCGTCGGGGGTAAATTCCAAACCGGGTGGTAATACCACCGGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTATACCATTCCAGGAAACGCTGCCAGCCATCTTCTGCAGATGCg                                                          <  1:1008856/150‑1 (MQ=255)
                                                                                              aaTGGTGAAGGAGTCGCAGAACTTGCGTCGTGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTGCTTCCCACCATCCTGCTTAATCCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg                                                         <  1:93732/149‑1 (MQ=255)
                                                                                               aTGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc                                                        <  1:1184688/150‑1 (MQ=255)
                                                                                                tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa                                                       >  1:786531/1‑150 (MQ=255)
                                                                                                tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa                                                       >  2:1569101/1‑150 (MQ=255)
                                                                                                   tGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTACACCATTCCAGCATACGCTGCCAGCCAGCTTCTGCAGATGCGGCATGa                                                    >  1:1971377/1‑150 (MQ=255)
                                                                                                   tGCAGGCGGCGCAGAACTTGAGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTATCTCATTACAGCGACTTCTTACCACCATACTGCTTAAACCATTACATCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAtca                                                    >  1:1943557/1‑148 (MQ=255)
                                                                                                     cAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCATCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTTCAGATGCGGCATGATa                                                  >  2:1474753/1‑150 (MQ=255)
                                                                                                       ggcggcGCAGAACTTGCGTCGTGGGTAAAATACAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGc                                                <  1:40270/150‑1 (MQ=255)
                                                                                                          ggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCg                                             >  2:1229443/1‑150 (MQ=255)
                                                                                                          ggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAGGATAGCTCg                                             >  2:1898589/1‑150 (MQ=255)
                                                                                                            cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGCCTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCgtg                                           >  2:200449/1‑148 (MQ=255)
                                                                                                             gcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGc                                          <  1:1471618/150‑1 (MQ=255)
                                                                                                              cAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCg                                         <  1:2205160/150‑1 (MQ=255)
                                                                                                               aGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGa                                        >  1:2014256/1‑150 (MQ=255)
                                                                                                                gAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGAt                                       >  1:29515/1‑150 (MQ=255)
                                                                                                                gAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGAt                                       <  2:1506910/150‑1 (MQ=255)
                                                                                                                  aCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATaa                                     <  2:2212781/150‑1 (MQ=255)
                                                                                                                   cTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAAt                                    >  2:1244360/1‑150 (MQ=255)
                                                                                                                   cTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCAGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACTCTGCCAGCCATCTTCTGCAGATTCGGCATGATAGCTCGGGCGATAAt                                    >  1:422830/1‑150 (MQ=255)
                                                                                                                    ttGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATc                                   >  1:1916295/1‑150 (MQ=255)
                                                                                                                    ttGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACACGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATATCTCGGGCGATAATc                                   >  1:1308870/1‑150 (MQ=255)
                                                                                                                     tGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCa                                  >  2:974327/1‑150 (MQ=255)
                                                                                                                          cGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGtt                             <  2:1199027/150‑1 (MQ=255)
                                                                                                                            ggggTAAATTCAAAAAAGGGTGGTATTACCACCCGGTCTTTTCTCATTACATCGATTTCTTCGCACCATACTGCTTAAACCATTCCAGCATACGCTGGCTGCCATCTTCTGCAGATGCGGAATGTTAGATCGGGAGATAATCTGCGTTGa                           <  2:659003/150‑1 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTGGTAATACCCCGCGGTGTTTTCTCATTACAGCGACTTCTTCCCCCCATACTGCTTAAACACTTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGATCGGGCGAGAATCAGCGTTGAatg                        >  1:308194/1‑150 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatg                        >  2:37042/1‑150 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatg                        <  2:933277/150‑1 (MQ=255)
                                                                                                                                tAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgc                       <  1:1100014/150‑1 (MQ=255)
                                                                                                                                 aaaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgca                      >  2:1513444/1‑150 (MQ=255)
                                                                                                                                   aaTCCAAACCGGGTGTTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgcatg                    <  1:1669649/150‑1 (MQ=255)
                                                                                                                                   aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCTTGATAGCTCGGGCGATAATCAGCGTTGAatgcatg                    <  1:1084851/150‑1 (MQ=255)
                                                                                                                                       ccaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCCTACTGCTTAAACCATTACAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCg                <  2:153649/148‑1 (MQ=255)
                                                                                                                                       cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCg                >  2:1136495/1‑150 (MQ=255)
                                                                                                                                        agaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATCATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCgg               >  1:2023868/3‑150 (MQ=255)
                                                                                                                                         gaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCCCCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGc              >  1:1921805/2‑150 (MQ=255)
                                                                                                                                         aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTTCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGc              <  2:149182/150‑1 (MQ=255)
                                                                                                                                           ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGt            >  2:676172/1‑150 (MQ=255)
                                                                                                                                           ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGt            >  2:116087/1‑150 (MQ=255)
                                                                                                                                           ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGt            >  2:539315/1‑150 (MQ=255)
                                                                                                                                           ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACAATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGTCGATAATCATAGTTGAATGCATGCCcgtcgt            >  2:746552/1‑150 (MQ=255)
                                                                                                                                            cGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTc           <  2:2053790/150‑1 (MQ=255)
                                                                                                                                                tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCggg       <  2:486205/150‑1 (MQ=255)
                                                                                                                                                tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCAGCCAGCCATCTTCTGCAGATGCGGCATGATAGCACGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCggg       >  1:617601/1‑150 (MQ=255)
                                                                                                                                                 ggtAATACCACCGGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGACAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGt      <  2:2096046/150‑1 (MQ=255)
                                                                                                                                                     aTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTacac  >  2:274298/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACAC  >  NC_000913/4015184‑4015482

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