Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,399,929 A→T L17Q (CTG→CAG)  mreC ← cell shape determining protein MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,399,9290AT97.7% 161.3 / NA 44L17Q (CTG→CAG) mreCcell shape determining protein MreC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (23/20);  minor base G (0/1);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTA  >  NC_000913/3399683‑3400168
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
tGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGCGCGCCACCAGCACCGCCTGAAt                                                                                                                                                                                                                                             >  2:301981/1‑251 (MQ=255)
tGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGCTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAAt                                                                                                                                                                                                                                             >  1:343840/1‑251 (MQ=255)
    tGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAgg                                                                                                                                                                                                                                         >  2:167866/1‑251 (MQ=255)
             cAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGt                                                                                                                                                                                                                                >  2:223851/1‑251 (MQ=255)
                        agTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGcc                                                                                                                                                                                                                     >  1:176009/1‑251 (MQ=255)
                           tCACGGTTTTTCAGCAACGGTTCCGGACGTAACGCCCGGTTTTCAAGTTCTAATTGGGCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCAc                                                                                                                                                                                                                   <  2:82867/250‑1 (MQ=255)
                                gTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCt                                                                                                                                                                                                               <  1:143840/249‑1 (MQ=255)
                                      cAGCAACAGTTCCGGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAAtt                                                                                                                                                                                                        <  2:248301/250‑1 (MQ=255)
                                      cAGCAACAGTTCCGGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAAtt                                                                                                                                                                                                        <  1:223851/250‑1 (MQ=255)
                                       aGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTgg                                                                                                                                                                                                      >  1:236610/1‑251 (MQ=255)
                                              gTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCAt                                                                                                                                                                                                >  2:102923/1‑250 (MQ=255)
                                                 ccTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATTCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCCCGGCTAAAAATTGGCTTCATAAGt                                                                                                                                                                                            >  2:286703/1‑251 (MQ=255)
                                                      cGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                                                                                                                                       >  1:65854/1‑251 (MQ=255)
                                                                           tCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAGAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcagg                                                                                                                                                                  >  2:183113/1‑251 (MQ=255)
                                                                               aTTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATCTAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCTCCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcaggcagg                                                                                                                                                              <  1:126626/251‑1 (MQ=255)
                                                                                              ccAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTc                                                                                                                                               >  2:300619/1‑251 (MQ=255)
                                                                                              ccAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTc                                                                                                                                               >  2:238465/1‑251 (MQ=255)
                                                                                              ccAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTc                                                                                                                                               <  2:167549/251‑1 (MQ=255)
                                                                                                aGCGTCGGCGATACGCCATCCAGCAATTCACGAGAAGCATTGGAAACAAAGAAGAAAGGACGGACGGCGGGATCCATAGAAGTACGGATTTGACTGAACACCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGGCGGAATAAGACGAATCTGTAGCGACGGTCCACGGCGAAAAATTGGCTTCATGAGTTATGCGTAGTGTCGTATCAGGCCATGCATGGTAAAAATACGCTTccc                                                                                                                                             <  2:32034/251‑1 (MQ=255)
                                                                                                                 aTCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCgctt                                                                                                                             >  2:198435/1‑247 (MQ=255)
                                                                                                                 aTCCAGAAATTCACGAGGAGAATTGGAAACAAATGAGAAAGGACTGACGGCGGTGTCCATATAAGTACGGATTTGACTGACCGTCCCCAGGCGGCTGTCGGCAAAAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACGTCCCCTGCCTGCATCCGATTa                                                                                                                            <  2:335249/251‑1 (MQ=255)
                                                                                                                          ttCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTCCTCTTCGCt                                                                                                                   >  2:117605/1‑251 (MQ=255)
                                                                                                                          ttCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCt                                                                                                                   >  2:354296/1‑251 (MQ=255)
                                                                                                                                      aTTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTcgcc                                                                                                       <  1:241330/251‑1 (MQ=255)
                                                                                                                                               aaattAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCCCCGCCTGAATAAGGCGAAGCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCAt                                                                                              <  1:232914/247‑1 (MQ=255)
                                                                                                                                                  gTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCCGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTc                                                                                           >  2:146980/1‑251 (MQ=255)
                                                                                                                                                  gTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTc                                                                                           >  1:6125/1‑251 (MQ=255)
                                                                                                                                                    aGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACGCTTCCCCTGCCTGCAGCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCg                                                                                          <  1:164642/250‑1 (MQ=255)
                                                                                                                                                               gcggcggTATCCGTATAAGTACGGATTGGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgc                                                                              <  2:274960/250‑1 (MQ=255)
                                                                                                                                                                         ccATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgc                                                                    >  1:323422/1‑251 (MQ=255)
                                                                                                                                                                               aaGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcg                                                               <  1:12591/250‑1 (MQ=255)
                                                                                                                                                                                                    tCCCCAGGCGGCTGTCGGCAATAATAATGACGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGt                                         >  1:79951/1‑251 (MQ=255)
                                                                                                                                                                                                        cAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATGTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTc                                     <  1:29397/251‑1 (MQ=255)
                                                                                                                                                                                                          ggcggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg                                   <  2:356685/251‑1 (MQ=255)
                                                                                                                                                                                                           gcggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGc                                  <  2:63360/251‑1 (MQ=255)
                                                                                                                                                                                                             ggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTTTCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGca                                 >  1:292530/1‑250 (MQ=255)
                                                                                                                                                                                                             ggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaa                                <  1:198435/251‑1 (MQ=255)
                                                                                                                                                                                                                 gTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCaacaac                            >  1:98717/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCTCCTCCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACt                        <  1:51850/251‑1 (MQ=255)
                                                                                                                                                                                                                      cAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTg                       <  1:167866/251‑1 (MQ=255)
                                                                                                                                                                                                                             aatGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGc                 <  1:300619/250‑1 (MQ=255)
                                                                                                                                                                                                                                       gcCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTcc      >  2:45500/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccaGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACCCTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTa  <  2:189127/251‑1 (MQ=255)
                                                                                                                                                                                                                                            caGCACCGCCTGAATACGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTa  >  1:293193/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
TGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTA  >  NC_000913/3399683‑3400168

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