Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,095 C→T intergenic (+229/+527) rtcR → / ← glpG DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,0950CT96.2% 160.4 / ‑3.9 55intergenic (+229/+527)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (23/28);  minor base G (2/0);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTG  >  NC_000913/3559861‑3560337
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
aGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATc                                                                                                                                                                                                                                    <  1:3994/251‑1 (MQ=255)
    ggctggccgctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAgg                                                                                                                                                                                                                                <  1:14388‑M1/35‑1 (MQ=255)
     gctggccgctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGc                                                                                                                                                                                                                               >  1:27246‑M1/216‑251 (MQ=255)
            gctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTg                                                                                                                                                                                                                        <  1:321175‑M1/43‑1 (MQ=255)
              tgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTggg                                                                                                                                                                                                                      <  1:5790‑M1/45‑1 (MQ=255)
                    aatcccgcctgaagggaaagctgcacgtaatcatcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactagcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATaa                                                                                                                                                                                                                >  1:326274‑M1/201‑251 (MQ=255)
                     atcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATaa                                                                                                                                                                                                                <  1:245481‑M1/51‑1 (MQ=255)
                      ccccgcctgaagggaaagctgcacgtaagcagcaagcgatatacgcagcgaattgagcggcataacctgaatccgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATaaaa                                                                                                                                                                                                              >  2:5790‑M1/199‑251 (MQ=255)
                          gcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaacCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCg                                                                                                                                                                                                           <  1:287960‑M1/56‑1 (MQ=255)
                                agggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggcagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATTGAGCGTCg                                                                                                                                                                                                    <  2:293851‑M1/63‑1 (MQ=255)
                                    aaagatgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTc                                                                                                                                                                                                >  2:14388‑M1/185‑251 (MQ=255)
                                       gctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCATAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGcc                                                                                                                                                                                             >  1:140466‑M1/182‑251 (MQ=255)
                                               taatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgcggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGCTAAACTCGAGCGTCGCTTCGCCAATGATcc                                                                                                                                                                                     <  1:79882‑M1/78‑1 (MQ=255)
                                                                     gcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagcttcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGTCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTc                                                                                                                                                               >  1:26636‑M1/152‑251 (MQ=255)
                                                                               cggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCg                                                                                                                                                     >  1:124146‑M1/142‑251 (MQ=255)
                                                                                   ataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaagaagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGaa                                                                                                                                                 >  1:209217‑M1/138‑251 (MQ=255)
                                                                                       cctgaatctgatgcagcacctggcacggctgggacggaattcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagatggagtcattacaCGATGCCGCTTATCTCCAGAATGCAGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGAGCCCGACATCGCGGCTGCGTAATTACCCAccggcgagactgtg                                                                                                                                             <  2:113197‑M1/118‑6 (MQ=255)
                                                                                           aatctgaggcagcacctggcacggctgggacggaagtcgatgtcgttctcaaaatcggtggcgctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGCTAAAATCGAGCGTCACTTCGCCAATGATCCCGCCCTCGCGGATGCGTAATTCGCCACCGGGCAGAATGATGcgg                                                                                                                                         <  2:272281‑M1/122‑2 (MQ=255)
                                                                                             tctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacatacaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCACGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGaa                                                                                                                                        <  1:323006‑M1/123‑1 (MQ=255)
                                                                                                gaggcagcacctggcacgggtgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagacagcgggcattagctgaacataaaacactagcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGaaaa                                                                                                                                      >  2:268390‑M1/125‑249 (MQ=255)
                                                                                                         cctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTAc                                                                                                                           >  1:178268‑M1/116‑251 (MQ=255)
                                                                                                           tggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcctgacaaagtcatcgggccttatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAAATGGGAGATAAAATCGAGCGTCGCTACGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGGCACCGGAGAGAATGATGCGAAAATCTTCTTTTAcc                                                                                                                          <  1:324049‑M1/137‑1 (MQ=255)
                                                                                                                   ctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACgtg                                                                                                                 <  1:244051‑M1/146‑1 (MQ=255)
                                                                                                                           gaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCg                                                                                                                                                     <  1:243752‑M1/110‑1 (MQ=255)
                                                                                                                           gaaggcgctgtcgttctcaaaatcgggggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCg                                                                                                                                                     >  2:243752‑M1/98‑207 (MQ=255)
                                                                                                                                cgctgtcgttctcaaaagcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTga                                                                                                    >  1:229575‑M1/93‑251 (MQ=255)
                                                                                                                                    gtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGAtt                                                                                                >  1:221570‑M1/89‑251 (MQ=255)
                                                                                                                                      cgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGAttgt                                                                                              <  2:227350‑M1/165‑1 (MQ=255)
                                                                                                                                                        agctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAAt                                                                            >  1:78987‑M1/69‑251 (MQ=255)
                                                                                                                                                        agctgcatgacaaagtcatcgggcattagctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCGCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAAt                                                                            >  1:9407‑M1/69‑251 (MQ=255)
                                                                                                                                                             catgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCtg                                                                       >  1:172918‑M1/64‑251 (MQ=255)
                                                                                                                                                              atgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCtgt                                                                      <  1:291109‑M1/189‑1 (MQ=255)
                                                                                                                                                                  caaagtcatcgggcattatctgaacagaaatcactatcaataagttggaggcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCGTCGCCAATGAGCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGGGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGAtt                                                                  >  2:245481‑M1/59‑251 (MQ=255)
                                                                                                                                                                     agtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGCAAATCTTCTTTTACCCTCAACGTGTTAGccacggtgggaatgttggtgacaaggtccaaattgatgtgtttgtg                                                               <  2:23847‑M1/196‑23 (MQ=255)
                                                                                                                                                                        catcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGc                                                              <  1:98332‑M1/197‑1 (MQ=255)
                                                                                                                                                                         accgggaatcatctgaacataaaactctatcaagaaggtggaggcattaccCTATGCCGCTTATCGCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAACCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCAGCAACGTGTTAGCAACGTTGATATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAgcaat                                                           >  2:288710‑M1/52‑249 (MQ=255)
                                                                                                                                                                           cgggcattatctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAgtgc                                                         <  1:37491‑M1/202‑1 (MQ=255)
                                                                                                                                                                                   atctgaacataacacactatcaatacacaggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTCATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCCACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGCGCAGTGCGTGCGCt                                                  <  1:110249‑M1/209‑1 (MQ=255)
                                                                                                                                                                                     ctgaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAcc                                               <  1:49278‑M1/212‑1 (MQ=255)
                                                                                                                                                                                       gaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAcc                                               >  1:293851‑M1/38‑249 (MQ=255)
                                                                                                                                                                                       gaacataaaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAcc                                               >  1:205084‑M1/38‑249 (MQ=255)
                                                                                                                                                                                       gaacataaaacactatcaagaagtgggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCAGCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAcc                                               >  2:332010‑M1/38‑249 (MQ=255)
                                                                                                                                                                                              aaacactatcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGc                                      <  1:114454‑M1/221‑1 (MQ=255)
                                                                                                                                                                                                   ctatcaataagtgggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGCGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGGGAGAATCTTCTTTTACCATCAACGTGGTAGCCACGTTGAGATTGTTGGTGACAAGGCGCCCATTGCTGTGATTGACCGGTGAGTGCGCTACCGCCTCcggcggggt                                 <  2:261499‑M1/226‑5 (MQ=255)
                                                                                                                                                                                                     atcaataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCgtggt                                 >  1:34359‑M1/24‑249 (MQ=255)
                                                                                                                                                                                                       caataagttggagtcattaccCAATGCCGCTTATCTCCAGAATGCAGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGCAATTCGCCACCGGCGAGAATGATGCGCAAATCTTCTTTTACCATCAACGTGTTAGAAACGTTGAGATTGTTGGTGACAATGCGCCAATTGCTGTGCTTGAGCAGTGCGTGCGctaccgcttccgccggggtgacc                             <  2:222976‑M1/230‑12 (MQ=255)
                                                                                                                                                                                                        aataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTCGCAACGTTGAGATTGTTGGTGACAATGCGCAAAtt                                                                           <  1:343682‑M1/184‑1 (MQ=255)
                                                                                                                                                                                                         ataagttggagtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGat                           <  1:274741‑M1/232‑1 (MQ=255)
                                                                                                                                                                                                            aggggggagcattaccCGATTCCGCTTTGCCCCCGAAAGCCGAAAGCAAGGCGGAACTGGGAAATAAAACCGGGCGTCGCTTCGCCAATGAGCCCGCCATCGCGTCTGCTTAATTCGCCATCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGGGTTAGCAACGTTGAGAGTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCTGCGTGGTGCCGATATc                        >  2:186309‑M1/17‑251 (MQ=255)
                                                                                                                                                                                                                   gtcattaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAAc                 >  1:227350‑M1/10‑251 (MQ=255)
                                                                                                                                                                                                                      attaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGc              <  1:28688‑M1/245‑1 (MQ=255)
                                                                                                                                                                                                                      attaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAAGCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGc              >  2:274741‑M1/7‑251 (MQ=255)
                                                                                                                                                                                                                       ttaccCGCTGCCGCTTATCTCCAGAATGCCGTAATCATGGCGGAACTGGGAGCTACAATCGCGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGTTGCGAAAATCTTCTTTTACCATCAACGTGTTCGCAACGTTGATATTTTTGGTGACATTGTGCACATTGCTGTGATTGAGCAGTGCGTGTGCTACCGCTTCCGGCGTGGTGCCGTTATCgataaaaag               <  2:10613‑M1/244‑4 (MQ=255)
                                                                                                                                                                                                                       ttaccCGATGCCGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCg             <  1:285439‑M1/246‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTTATCTCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTg  >  1:204956/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
AGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTG  >  NC_000913/3559861‑3560337

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