Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,154,407 A→T 9.8% I271N (ATC→AAC)  yqhC ← transcriptional activator of yqhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,154,4070AT 9.8% 141.7 / 2.1 52I271N (ATC→AAC) yqhCtranscriptional activator of yqhD
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base T (2/4);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGC  >  NC_000913/3154184‑3154632
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                 
tctgttGGGGTTAATAAAAAAAAATACTTTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGGATTTAATTCCCCTGCATCTCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGatcat                                                                                                                                                                                                                               <  1:650118/224‑1 (MQ=255)
 tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCCTGCCGCCATGGATGatcacc                                                                                                                                                                                                                              >  1:375421/1‑226 (MQ=255)
            aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCCCACCGAAGTAACGTTTAAACTCACGCCTAAATTGCGATGCGCTTTCATAGCCGACTCGCCTCGCTGCTCCCCTGGCCTACATCCCGTCATGGGTGTTCACCATCCCCGCCt                                                                                                                                                                                                                    >  2:208354/1‑227 (MQ=255)
             ttAAAAAAAAATACTGTAACGCCTGACCATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCtt                                                                                                                                                                                                                   <  1:208354/227‑1 (MQ=255)
                               aCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAAttcttc                                                                                                                                                                                                <  1:513458/228‑1 (MQ=255)
                                     gACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCCTAGCCGACGCGCATCGCTCCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCCCCTTATGCTAACGGAAATTCTTCAAATAc                                                                                                                                                                                          >  2:198253/1‑228 (MQ=255)
                                          ttttCCCCGTTCCCGGTTTCTGTGCCGGGAACCTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCCAACTTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCaa                                                                                                                                                                                     <  2:200757/228‑1 (MQ=255)
                                             tCCCCGTTCCCGGTTGCTTTCCCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCGCGGCTAAATTGCGATGCGCTGTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTGATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACgg                                                                                                                                                                                  <  1:36697/228‑1 (MQ=255)
                                                     cccGGTTTCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTAc                                                                                                                                                                          <  1:251871/228‑1 (MQ=255)
                                                     cccGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTAc                                                                                                                                                                          <  1:155050/228‑1 (MQ=255)
                                                           gcatGTACCGGGAACGTATGTAATTCCCCGGCATCGCCCGCACACGTGCCGCATCGTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCGTAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGa                                                                                                                                                                    <  1:436365/225‑1 (MQ=255)
                                                                tACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGCCAGAc                                                                                                                                                               >  1:85125/1‑228 (MQ=255)
                                                                tACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGGTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGCCAGAc                                                                                                                                                               >  1:696525/1‑228 (MQ=255)
                                                                          gTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCGTAGCCGACGCGCCTCGCTGCTGCGCTGGCATTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAAACGGTAATTCCTCAAATACTGCAACGGCGAGGAACTGGGCACAGACTTaaaaatgc                                                                                                                                                     >  2:436365/1‑224 (MQ=255)
                                                                          gTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTCCGCTGTCCTTCATGCCGTCATGGATGATCATCCTCCGCGCCTTTTGCAGTCGGTAATTCTTCAAATACTGCACCGGCGCGGGACCGGGGACAAACTTAAAAtcaa                                                                                                                                                     >  2:574542/1‑225 (MQ=255)
                                                                                     cccTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGGCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGcc                                                                                                                                          >  1:416972/1‑228 (MQ=255)
                                                                                     cccTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGTGGGAACGcc                                                                                                                                          >  1:68616/1‑228 (MQ=255)
                                                                                         gccgcGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCCCACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATa                                                                                                                                      <  1:462095/224‑1 (MQ=255)
                                                                                               cccGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGACGGCCGCTACGc                                                                                                                                   >  1:585406/1‑225 (MQ=255)
                                                                                               cccGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGccc                                                                                                                                 >  1:510145/1‑225 (MQ=255)
                                                                                               cccGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGACTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGTTGATCATCATCCGCGCCTTATGCCGCCGGTAATTCTTCACATACTGCAACGGCCAGGTACTTGTCGCAGAATTAAAATTATGGGGGAACGCCGATaacctc                                                                                                                                 >  2:167970/1‑222 (MQ=255)
                                                                                                    aTTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCCGACGGTAATTCTTCAAAt                                                                                                                                                                                            >  2:461447/1‑163 (MQ=255)
                                                                                                    aTTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCCACGCCGAGGTACTGGTGCCAGACTTCAAATTATGGTGGAACGCCGATACGCTCATCtt                                                                                                                            >  1:360774/1‑227 (MQ=255)
                                                                                                     ttCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATa                                                                                                                                                                                           <  1:461447/163‑1 (MQ=255)
                                                                                                          gCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTCAAAATTATGGGGGAACGCCGTTACGCCCATGTTGGCTTc                                                                                                                      >  1:217169/1‑227 (MQ=255)
                                                                                                                     ccccGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTg                                                                                                          >  1:448230/1‑228 (MQ=255)
                                                                                                                            gTCACACCGAAGTAACGTTTAAACCCGCGGCTAAAGTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCTTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACg                                                                                                   <  1:697083/228‑1 (MQ=255)
                                                                                                                               acacCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCCCGCCTTCTGCAGCCGGTACTTCTTCAAATACTGCAACGGCGAGGTACTGGTGACCGCCTTAAAATTATGGTGGAACGCCGACACGCCCATGTTGGCTTCTGCCCCAAGTGGCTCGagcccc                                                                                                >  2:410388/1‑223 (MQ=255)
                                                                                                                                      cgTAACGTTTAAACTCCGGGGTCAATGGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGTTGGCCTTCATGGCGGCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGtttt                                                                                         <  2:416972/227‑1 (MQ=255)
                                                                                                                                        tAACTTTTAAACCCACGGCTAAATTGCGATGCGCTTTCATAGCCGAGGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCg                                                                                       <  2:375421/228‑1 (MQ=255)
                                                                                                                                         aaCGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCgg                                                                                      <  1:167970/228‑1 (MQ=255)
                                                                                                                                         aaCGTTTAAACCCGCGGTTAAATTGCGGTGGGTTTTCCAGGCCGACCCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCgg                                                                                      <  2:85125/228‑1 (MQ=255)
                                                                                                                                          aCGTTTAAACCCACGGCGTAATTGCGTTGGGGTTTCATTGCCGACGCGCATCGCTGCTGCGCTGGCCTTCAGTCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTGCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGgt                                                                                     <  2:274655/228‑1 (MQ=255)
                                                                                                                                           cGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATCCCGTCATGGATGAACATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATCCTGCCACGGCGAGGTACTGGTGACAGACTTCAAATTATGGTGGAACGCCGACACGCTCATCTTGGCTTCTTCCGCCAGTTGCTGGACGCCCCGGATTTCggcg                                                                                    >  2:96995/1‑226 (MQ=255)
                                                                                                                                               tAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGACCGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCCCAGGTTTTCGGTGaatt                                                                                >  1:70939/1‑228 (MQ=255)
                                                                                                                                               tAAACTCACGGCTAAAGTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGtatt                                                                                <  1:198253/228‑1 (MQ=255)
                                                                                                                                                      aCGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGTTCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTa                                                                                                                                                            >  1:446679/1‑145 (MQ=255)
                                                                                                                                                      aCGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATAAGCTCATGTTGGCTTCTGCCGCCAGTTGATGGAACCTCAGGTTTTCGGGGTATTTTTtctc                                                                         >  2:377224/1‑228 (MQ=255)
                                                                                                                                                       cGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGTTCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTaa                                                                                                                                                           <  2:446679/145‑1 (MQ=255)
                                                                                                                                                             aaTTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAACTTATGGGGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCACGTTTTCGGGGTATTTATTTTCAATCCGt                                                                  >  2:12326/1‑228 (MQ=255)
                                                                                                                                                                    aTGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCCAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTTTTCTCAATCCGTTTCAGCa                                                           >  1:679466/1‑228 (MQ=255)
                                                                                                                                                                              aaggCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCAAGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATGCTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTaa                                                  <  1:410388/224‑1 (MQ=255)
                                                                                                                                                                                   ctaCGCCCATCGCTGCTGGGCTGGCCGTCATGCCCTCATGGATGATCATCATCCGCGCCGTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATTCGCTCATGTTGGCCTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGAc                                            <  2:68616/226‑1 (MQ=255)
                                                                                                                                                                                      cgcgCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGGTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTg                                          >  1:364519/1‑227 (MQ=255)
                                                                                                                                                                                               tctctgcGTTGGCCTTACATCCTTCATGTACTTTCCTCACCCGGGCCGTATGCAGACGGTAATTCTTAAAATCCTGCAACGGCGAGGTACTGGTGACAGATTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCt                                <  2:696525/225‑1 (MQ=255)
                                                                                                                                                                                               ctttttcgcTGGCTTTAATCCCGTCATGGATGTTACTCAGCCGCGCTTTTTGCAGATGGTAATTCTTCAAATGCTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCt                                <  2:585406/222‑1 (MQ=255)
                                                                                                                                                                                                     cgcTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGAc                          >  1:385972/1‑228 (MQ=255)
                                                                                                                                                                                                            cTTTATGCCGTTAATGATGTTCATCAGCCGGGCCTTATGCAGACGGTAATTCTTCAAATTCTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGGGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAg                   <  2:631689/228‑1 (MQ=255)
                                                                                                                                                                                                             ttCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGtttt                                                                                         <  2:402473/157‑1 (MQ=255)
                                                                                                                                                                                                             ttCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGtttt                                                                                         >  1:402473/1‑157 (MQ=255)
                                                                                                                                                                                                                cATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCgg                                                                                      <  1:685527/158‑1 (MQ=255)
                                                                                                                                                                                                                cATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGGTTTCgg                                                                                      >  2:685527/1‑158 (MQ=255)
                                                                                                                                                                                                                     cGTCATGGATGATCATCATCCGCGCCTTATGCAGCCGGTAATTCTTCAAATACTGCAACGGCGAGGTCCTAGTGACAGACTTAAACTTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTTCCGCCAGTTGCTCGACGCTCCAGCTTTCGGTGTATTTATTCTCCACCCGGTTCAGCACTCGCCTAATCCGACTTAAGTGATTCTTGCGACGTACCACTGCCAGt            >  2:294646/1‑227 (MQ=255)
                                                                                                                                                                                                                          tGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAAcgcg      >  1:527673/1‑228 (MQ=255)
                                                                                                                                                                                                                               gatCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTAAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCGTGCGACTAACCAGCGCCAGTAACGccccgc  >  1:43451/1‑227 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                 
CTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAGTAACGTTTAAACTCACGGCTAAATTGCGATGCGCTTTCATAGCCGACGCGCATCGCTGCTGCGCTGGCCTTCATGCCGTCATGGATGATCATCATCCGCGCCTTATGCAGACGGTAATTCTTCAAATACTGCAACGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGC  >  NC_000913/3154184‑3154632

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