Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,225,170 G→A F154F (TTC→TTT yohC ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,225,1700GA84.0% 117.0 / 14.5 50F154F (TTC→TTTyohCputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (22/20);  total (26/24)
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.57e-01

GATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGC  >  NC_000913/2224925‑2225402
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
gATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCt                                                                                                                                                                                                                                     >  2:353866/1‑251 (MQ=255)
 aTTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTc                                                                                                                                                                                                                                    <  1:141135/251‑1 (MQ=255)
    tACTGATGCTCGAGACGCCAGCATAAGGAGCAAAAAGTGCGAGCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACGTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGATCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAATGCCTGCCTCTTTGTTGATATGCAGAAAGCTCgg                                                                                                                                                                                                                                  <  2:93104/250‑1 (MQ=255)
        gAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAAt                                                                                                                                                                                                                               <  2:137471/249‑1 (MQ=255)
                cGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTaa                                                                                                                                                                                                                     <  1:321382/251‑1 (MQ=255)
                       atcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGt                                                                                                                                                                                                              <  1:47285/251‑1 (MQ=255)
                         catAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAg                                                                                                                                                                                                            <  2:403772/251‑1 (MQ=255)
                            aaCGAGCAAAAAGTGCGATCGTCAAAGGAGACTCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCg                                                                                                                                                                                                         >  2:383343/1‑251 (MQ=255)
                                    aaaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGAAGGTAGCCGGTGTaaa                                                                                                                                                                                                  >  2:42257/1‑250 (MQ=255)
                                          gcgctCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGGACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCGTCCTCTTGGTTGATATGCAGGAAGCTCGGAATACCAAGTTAAAGCAGGGAGCCGGTGTAAAACAgcg                                                                                                                                                                                            <  2:268672/246‑1 (MQ=255)
                                           cGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGcc                                                                                                                                                                                          <  1:106535/251‑1 (MQ=255)
                                             aTCGTCAAAGGCGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCAc                                                                                                                                                                                        >  1:53051/1‑251 (MQ=255)
                                                gTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCATGGAGTCGGGGTAAAACATCTCCACCGt                                                                                                                                                                                     <  1:334749/251‑1 (MQ=255)
                                                gTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGt                                                                                                                                                                                     >  1:45252/1‑251 (MQ=255)
                                                        agaATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGAcc                                                                                                                                                                              >  1:341219/1‑250 (MQ=255)
                                                                     aaTGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAacc                                                                                                                                                                 >  1:26429/1‑250 (MQ=255)
                                                                           gcATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCTGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAaccagacc                                                                                                                                                            <  1:385453/249‑1 (MQ=255)
                                                                             aTTGCTGACATTTACGCCAGCAATGCAACGGCAAAGAACTTTCTAGAACAAACGATAACCGTAGACCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTg                                                                                                                                                        <  2:81345/251‑1 (MQ=255)
                                                                                gCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTATAGAACACACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGAAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGAt                                                                                                                                                     >  1:269596/1‑251 (MQ=255)
                                                                                             ccAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATTACGGTTAGGGCGATCAGCACGTCCATTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCCTCCTCTTTGTTGATATTCAGAAAGCGCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGGGCCGGGCAGCGCACACAACCAGAACAGTGGATAAAGCGCCACCAg                                                                                                                                        <  1:432744/251‑1 (MQ=255)
                                                                                              cAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGccacc                                                                                                                                          >  1:156989/1‑248 (MQ=255)
                                                                                                  aaTGCAGCGTCAGAGAACTTTCTAGAACAACCGATAACCGTATCCCAATACTATCACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACCCCAATCGCCAGCGTCGACCTAGAGACGCTAAGGCCTCCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAACGCGCCACCAGACCGc                                                                                                                                   >  2:155086/1‑251 (MQ=255)
                                                                                                  aaTGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGCGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGc                                                                                                                                   >  1:54643/1‑251 (MQ=255)
                                                                                                   aTGCAACGTCAAAGAACTTTCTAGAACAAACGATAACAGTAGCCACATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAACAGAACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGACCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCACTGGATAAAGCGCCACCAGCCCGCt                                                                                                                                  >  1:205282/1‑251 (MQ=255)
                                                                                                          gTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATCGTATAACGGTTAGGGCGAGAAGCACTTCCAGTACCAGCACGCCACTCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATAACAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAAc                                                                                                                           >  2:405563/1‑251 (MQ=255)
                                                                                                          gTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAACACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAAc                                                                                                                           <  1:237503/251‑1 (MQ=255)
                                                                                                                      tCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGc                                                                                                                <  1:26356/250‑1 (MQ=255)
                                                                                                                        tAGAACAACCGATCACCGTCGCCCCATAGTATAACGGTTGGGGCGAGCAGCACTTCCAGTACCAGCACTCCAATCGCCAGCGTCGAACTAGCGGAGCTAAGGCCTTCCTCTTTGCTGATCTTCAGAAAGCTAGGACTACCAAGGTAAAGCAGGAAGCCGGTGCAACACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGAGGACAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTGGCAAc                                                                                                             >  2:289894/1‑251 (MQ=255)
                                                                                                                          gAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGt                                                                                                           <  1:283994/251‑1 (MQ=255)
                                                                                                                          gAACAAACGATAACAGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGTCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGCTAAAGAGCCACCAGACCGCTTAAAAACAGCGGAGTGGCAACGt                                                                                                           >  2:290202/1‑251 (MQ=255)
                                                                                                                           aaCAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGt                                                                                                           >  2:326528/1‑250 (MQ=255)
                                                                                                                               aaCGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGTTGTAAAACAGCTCCACCGTGCCGACCATCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGGTGCAACGTAGccc                                                                                                      <  1:206282/251‑1 (MQ=255)
                                                                                                                                                   aGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTgcgc                                                                                  >  1:31437/1‑251 (MQ=255)
                                                                                                                                                                gggcgtgaaaatcTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGACTTCCTATTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTCGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACCACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAACACAGCGGAGTTGCACCGTAGCCCGCGACGACCAGGCAGTGCGCTAGTGACGGACGc                                                                     >  2:387391/13‑251 (MQ=255)
                                                                                                                                                                    gAGCAGCACTTCAAGTACCAGCAGGCCAATCGCCCGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGATGATCTTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGCGCCGACCAGCGCACACAACCCGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCg                                                                 >  2:274369/1‑251 (MQ=255)
                                                                                                                                                                        agcaCTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATa                                                             <  1:42257/251‑1 (MQ=255)
                                                                                                                                                                               ccAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTAcgc                                                      <  1:168132/251‑1 (MQ=255)
                                                                                                                                                                                cAGTACCAGCACGCCAATCGCAAGCGTCGCACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTAcgcg                                                     >  2:71130/1‑251 (MQ=255)
                                                                                                                                                                                     ccAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATc                                                >  1:423900/1‑251 (MQ=255)
                                                                                                                                                                                           cGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATccacc                                            >  2:144726/1‑249 (MQ=255)
                                                                                                                                                                                             ccAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGAGGCTGCGGAGAATTACGCGCTATCCGCCAGAt                                        <  2:269596/251‑1 (MQ=255)
                                                                                                                                                                                             ccAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGGGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGAt                                        <  1:405563/251‑1 (MQ=255)
                                                                                                                                                                                              cAATCGCCAGAGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAATTTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCCGTGGATAAAGCGCCACCCGACCGCTTCACAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTCCGCAAGTGACGGACGCTGCGGATAATTACGCGCCACCCACCAgcgg                                       >  2:169505/1‑248 (MQ=255)
                                                                                                                                                                                                aTCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGGCCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGTAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGt                                                                              <  2:152689/210‑1 (MQ=255)
                                                                                                                                                                                                aTCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGCAGACCATGCAGTGCGCAAgc                                                                              >  1:152689/1‑209 (MQ=255)
                                                                                                                                                                                                         gTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGATGACCCGCCCCAt                            <  1:353866/251‑1 (MQ=255)
                                                                                                                                                                                                                tAGAGAAGCTAAGGCATTCCTCTTTGTTGATATTAAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCTCCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACAAGACCGCTTAAAAACAGCGTAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGAGCAAGTGACAGACGCTGCGGATAATTACGCGCCAGACACCAGATGACCCGCACCATCACcgcg                     >  2:311575/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaGCTAAAGCTTTCCTCTTTCTTGATATTCATAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTCAAACAGCGCCACCGTGACGACCAGCGCACACACCCAGACCAGTGGATAAAGCGCCACCAGACCGCTTACCACCAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGATGACCCGCCCCATCCCCGCGAcccc                >  2:364768/1‑251 (MQ=255)
                                                                                                                                                                                                                      aGCTAAGGCCTTCCTCTTTGTTGATATTCAGAAAGCTCGGAAGACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCCCCAGACCGCTTAAAACCAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGATGACCCGCCCCATCACCGCGACCcca               <  2:45252/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cgctttggttattttCAGAAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGGGGATAAAGCGCCCCCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGATGACCCGCCCCATCGCCGCGACCCCAGCCAGCATAACGc  <  2:53051/238‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
GATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGATATTCAGGAAGCTCGGAATACCAAGGTAAAGCAGGTAGCCGGTGTAAAACAGCGCCACCGTGCCGACCAGCGCACACAACCAGACCAGTGGATAAAGCGCCACCAGACCGCTTAAAAACAGCGGAGTTGCAACGTAGCCCGCGAAGACCATGCAGTGCGCAAGTGACGGACGCTGCGGATAATTACGCGCCATCCACCAGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGC  >  NC_000913/2224925‑2225402

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