Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,838,902 C→T 11.1% L219L (CTG→TTG)  yicL → EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,9020CT11.1% 165.9 / 3.0 73L219L (CTG→TTG) yicLEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base C (18/46);  new base T (7/2);  total (25/48)
Fisher's exact test for biased strand distribution p-value = 6.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

ACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCATCGGTGATTTTGATTTCAATGGA  >  NC_000913/3838680‑3839122
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  
caccccGGTTATCCGTCGTCATTTTCGATCTCCCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCg                                                                                                                                                                                                                               <  2:415436/223‑1 (MQ=255)
       gTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAGCAACCTTTGTCGTTAACGGCAGTTTGATTCTGGCGttttt                                                                                                                                                                                                                          >  1:90655/1‑226 (MQ=255)
          attcgtcGTCATTATCGGTCCCTCCTTGCGCGTTGTTCTGGGGCCTTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGTTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑tt                                                                                                                                                                                                                    <  2:394668/222‑1 (MQ=255)
             cgcattCATTTACGAACTCTCCGGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACGTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTtg                                                                                                                                                                                                                  <  1:424276/222‑1 (MQ=255)
                cGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGTGTATGCTGATTGGCGGGCTGATTCTGTTGCCTTTTTATGCCAGACACGGAACAAACTTTGTCGTTAACGGCAGGTTGATTTTGGCTGTTTTTTA‑ATtggtgg                                                                                                                                                                                                              >  1:369836/1‑228 (MQ=255)
                cGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGGTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTCACGGCAGGTTGATTCTTGGGTTTTTTTA‑TTttgtgg                                                                                                                                                                                                              >  1:575527/1‑228 (MQ=255)
                  tCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAAATTTGTCGTTAACGGCCGTTTGGTTCTGGCGTTTTTTTA‑TTTtgtggt                                                                                                                                                                                                             >  1:208018/1‑227 (MQ=255)
                          gATCTCTCCTGCCGCCTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTAc                                                                                                                                                                                                    <  1:378881/228‑1 (MQ=255)
                             ctctcCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTTCCTTTTTATGCCCGACAAGGACCAAACATTGTCGTTAACGGCGGGTTTATTTTGGGGCTTTTTTT‑TTTGTGGGGCATTgtgaagcg                                                                                                                                                                                                 >  2:204961/1‑221 (MQ=255)
                              tcctcTGCCGTGTTTTTCTGGGGTATTTCCTTGGCATTTGCTGCGGCATTCTAGTCCACCTTCCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTTATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTcc                                                                                                                                                                                                <  2:369836/224‑1 (MQ=255)
                                tcCTGCCGCTTTTGTCTTGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATTCCACCTATCCCTCAACGCTAATTGCCCGGTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCt                                                                                                                                                                                              <  1:308129/228‑1 (MQ=255)
                                cccTGCCGTTTTTTTCTGGGGCATTGCCCCGGCCTTTGCTGCTGCATTCTATACCGCCTATTCCTCAACGTTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCt                                                                                                                                                                                              <  2:90655/227‑1 (MQ=255)
                                     ccGCTTTGGTCTGGGGCATTGCCTGGGCATTTGCTGCTGCATTCTATACCACCTCTCCGCCAACGATAATTGCTCGCTATGGCTCGTTACCAGTCGTCGGCGGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACAt                                                                                                                                                                                         <  2:572384/228‑1 (MQ=255)
                                     ccGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTTTCGTTAACGGCCGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTTGTACGGTCCTGACa                                                                                                                                                                                          >  1:482550/1‑227 (MQ=255)
                                          ttgttTTGGGGCATTGCCTGGGCTGTTGCGGTTTCCTTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGGCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTGCGTCCCTGACATTTAGt                                                                                                                                                                                    <  2:365437/228‑1 (MQ=255)
                                          ttgttCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGGCTTATTCTGTTGCCTTTTTATGCCAGACAAGGGACAAAGCTTTTCGGTTACGGGCGGTTTTTTTGGGCGTGTTTTTACTTTGTGGGTCATT‑GTACGGcccccgactttatt                                                                                                                                                                                    >  2:307096/1‑217 (MQ=255)
                                                       ttGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCCTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCCGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACACTTTGTC‑ttgtactcaaag                                                                                                                                                                      >  1:404382/1‑217 (MQ=255)
                                                          ccTGGGCATTTGCTGCTGGCTTCTATACCCCCTGTCCCTCAACGGTAATTGCCCGGTATGGCACGTTTCCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAg                                                                                                                                                                   <  1:224254/228‑1 (MQ=255)
                                                                    tgctgctgcATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTTTTGCGTTTTTTTC‑TTTGTTGGTCATTGGTACATCTCTGACATTTACAC‑TTTACCTGAAAGAAGCTCaatt‑aat                                                                                                                                                        >  1:296076/1‑228 (MQ=255)
                                                                     gctgctgcATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTTTTGGGGTTTTTTA‑TTTGTTGGTCCTTGGTAAGTCCCCGTACTTTAGTC‑TTTTACTTGAAGGAGCCCaatt‑aatt                                                                                                                                                       >  1:567725/1‑228 (MQ=255)
                                                                                   taCCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGCCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTTGGGGTCCTTGGTACGTCCCTGAACTTTAGTC‑TTTACTCTAAAGGAGCACATTT‑AATTGTC‑GG‑TCCaaaacc                                                                                                                                       >  1:277439/1‑226 (MQ=255)
                                                                                   taCCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGAAAAGGAACAAACTTTGTCGTTAACGGCAGTTTGGTTCTGGCGTTTTTTTA‑TTTGGTGGGCCTTGGTCAGTCCCTGGCATTTAATC‑CTTACCTGGAATGAGCACAATT‑AATTTGC‑GG‑TCCtaaaac                                                                                                                                       >  1:300748/1‑226 (MQ=255)
                                                                                     ccaccTACCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCATTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCa                                                                                                                                     <  1:500913/228‑1 (MQ=255)
                                                                                        ccTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAACATTTGTCGTTAACGGCAGTTTGCTTTTGGCGTTTTTTTA‑TTTTGTGGTCAATGGTACGTCCCTGACCTTTTGTC‑TGTACCTTAAAGGAGCACAATT‑TATTCGC‑GT‑GCCAAAAGCCacgc                                                                                                                                  >  1:564950/1‑225 (MQ=255)
                                                                                        ccTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTAcc                                                                                                                                                                            >  1:357953/1‑189 (MQ=255)
                                                                                         cTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCt                                                                                                                                                                           <  2:357953/189‑1 (MQ=255)
                                                                                          tATCCCTGATCGCTAATTGCCCGGTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTGGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCAtt                                                                                                                                <  2:115905/228‑1 (MQ=255)
                                                                                           aTCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTTTCGTTAACGGCTGTTTTTTTTTTGCGTTTTTTTT‑TTTTGTGTTCATTTGTACGTTCCCTACCTTTTGTC‑TGTTCCTTAAAGGGGTGCCATT‑TATTGTCGGG‑TCCaaacccacacttt                                                                                                                                >  2:648343/1‑218 (MQ=255)
                                                                                           aTCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTCTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGGTTTTTTA‑TTTGGTGGGCACTGGTACGTTCCTTGCCTTTATTCTTATACCTAAAGGGGGCA‑AATT‑CATGTGC‑GGATCCAAAA‑CCAGCActt                                                                                                                               >  1:374201/1‑228 (MQ=255)
                                                                                               ttcaattttAATGGGCCGCTATGGCACGTTACCAGTCGTCGTCTGGATTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGGCAAGGAACAAGCTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAg                                                                                                                           <  2:175981/220‑1 (MQ=255)
                                                                                               cTCAATGTTTATTGCCCGTTTTGGCCCGTTACCTGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAg                                                                                                                           <  2:296076/228‑1 (MQ=255)
                                                                                               cTCAACGCTAATTGCCCGGTATGGCTCGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTTGGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAg                                                                                                                           <  1:272287/228‑1 (MQ=255)
                                                                                                tCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCCGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTGTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGc                                                                                                                          <  1:105087/228‑1 (MQ=255)
                                                                                                  aaCGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCAAGAAAAGGAACAAACTTTGGCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTTCATGGTACTCCCCTGCCCTTTAGTC‑TGTACCTGAAAGGAGGCCAATT‑AATTGGC‑GG‑TCCACAAACCAGCATTTTGAGCtg                                                                                                                        >  1:406778/1‑228 (MQ=255)
                                                                                                   aCGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGcac                                                                                                                                                                <  1:77354/190‑1 (MQ=255)
                                                                                                   aCGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCGGTTGCCTTTTTATGCCAGGCAAGGAACAAACTTTGTCGTTATCGGCAGTTTGATTCTGGCGTTTTTTTA‑TTCGGTGGTCATTGGTACGTCCCt                                                                                                                                                                                              <  1:422262/161‑1 (MQ=255)
                                                                                                     gCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGc                                                                                                                          <  1:202514/223‑1 (MQ=255)
                                                                                                       ttaTTGCCCGCTTTGGCACGTTACCAGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  2:482550/226‑1 (MQ=255)
                                                                                                       tAATTGCCCGTTTTGGCATGTTCCCCGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  2:552625/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGGTTGGGCAAGTGAACAGGCGTCGGCTGGAGTATGCTGATTGGCGGTCGGATTATTGTGCCTTTTTATGCCAGACGAGGAACAAACTTTGGCGATAACGGCAGTTTGAGTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑CGTACCTGAAAGGAGCACAATT‑AAGTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:648343/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGCTTGGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:132120/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGCTTAGTCCTGTTGCCAGTCGTGGGCTGGAGTATGCCGATTGGCGGGTTGATTCGGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATCCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTAGGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCGGTGCAg                                                                                                                   <  2:404382/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGATATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:22206/228‑1 (MQ=255)
                                                                                                         actGCCCGCTATGGCCACTTTCCAGTCTTCGGCTGGATTATGCTGATTGGCCGTCTGATTCTGTTGCCTTTTTATTCCAGACAAGGAACAATCTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑GTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGaa                                                                                                                 <  1:409381/226‑1 (MQ=255)
                                                                                                         aTTTCCCGTTATGGCACGTTACCAGTCGTCGGGTGGAGTATGCTTATTGGCGGTCTTATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACGTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGaa                                                                                                                 <  2:567725/228‑1 (MQ=255)
                                                                                                                gttttggatcGTTCCCAGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATTGCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTaa                                                                                                          <  2:494941/219‑1 (MQ=255)
                                                                                                                          gTTTCCAGTTGTCGGCTGGGGTATGCTTATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCt                                                                                                <  2:487179/228‑1 (MQ=255)
                                                                                                                            tACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTACTTGGGAATCCCCTGCATTTTATTC‑TTTACCTGAAAGGGTGCCAATTAAATTGGG‑GG‑TCCAAAACCCGGC‑TTTTTAGCTTTGCGAAAACAATAAGAAgagcgcaaa                                                                                              >  2:150447/1‑225 (MQ=255)
                                                                                                                             tcccgttgtcgGCTGGGTTTGGCTGGTTGGGCGTCTGATGTCGTTGCCGTTTTCTGCCAGAGCAGGACCGAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGTATTTTGAGCTGTGCAGAACCATTAAGTAGCGGGCTAct                                                                                             <  2:124269/221‑1 (MQ=255)
                                                                                                                             gccAGTCGTCGGCTGGAGGATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAATGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTAc                                                                                              <  1:274961/226‑1 (MQ=255)
                                                                                                                                  tcgtcgGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACGTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:345341/201‑1 (MQ=255)
                                                                                                                                    gtcgGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCCTGCATTTTATC‑CGTACCTTTAAGGGGGCAAATT‑CATTTGC‑GG‑GCCCACAACCCGCATTTTTGCCTTTGCAGAACCCTTTAGAAacgcgccaaaccctttt                                                                                      >  2:122296/1‑217 (MQ=255)
                                                                                                                                      cgGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGTTTCTGGCGTTTTTTTT‑TTTTGGTGTCATTTGGACGTTCCTGACATTTTTGC‑TGTATCTGAAATGAGCCCAATT‑AATTTGC‑GG‑TCCCAAAGCCCGCCTTTTTGGCCGGGCAGAAACCTTAAAGAGAGCGCTACCCCCtttttt                                                                                    >  2:48841/1‑225 (MQ=255)
                                                                                                                                       tgtttgAGTATCGTGATGGGCGGTCTGTTTCTGGTTGTTTTTTATGCCGTACAAGGAACAAACTTTGGCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTTAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTttgctg                                                                                   <  2:51460/223‑1 (MQ=255)
                                                                                                                                             aGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGaa                                                                                                                 <  1:115655/192‑1 (MQ=255)
                                                                                                                                                  gCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTcc                                                                                                                                                                                                <  1:168462/112‑1 (MQ=255)
                                                                                                                                                           gcgTCTTATTCTTTTGCCTTTTTATGCCAGACAAGGAAAAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTca                                                               <  1:204961/225‑1 (MQ=255)
                                                                                                                                                             gTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGGTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTTGC‑GG‑GCCAAAAACCAGCATTTTTAGCTGTGACGAACCATTAAGTAGCGCGCTACTCTCCTTTCTGTTGTTGGGGA‑TCaaggtccca                                                             >  1:439758/1‑219 (MQ=255)
                                                                                                                                                               ttgtttCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTa                                                          <  2:204990/224‑1 (MQ=255)
                                                                                                                                                                ttaTTCTGTTGCCTTTTTATGCCAGTCAAAGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCCCGTTCACAttcc                                                         <  1:149311/226‑3 (MQ=255)
                                                                                                                                                                    tCTGGTGCTTTTTTAGGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑CCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGa                                                     <  2:208018/228‑1 (MQ=255)
                                                                                                                                                                        ttGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTAAGTCCCTGCCATTTAGTC‑TGTACCTGAAAGGAGCCCAATT‑AATTGGC‑GG‑TCCCAAAGCCATCATTTTGAGCTTGGAAGAACCATTAAGTAGCGCGCTACTCCCTTTGCGTTTGCTGGGGAGTCA‑GTTCACAGTAACCGAccag                                                 >  2:199436/1‑225 (MQ=255)
                                                                                                                                                                             tttttATTCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTTTTGCTTTTTTTTA‑TTTGGTGGTCATTGGTGCGTCCCGGCCATTTAGTC‑TGGACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGTA‑TCACGTTCACATTACCGGACTGGCTggg                                            <  2:406778/228‑1 (MQ=255)
                                                                                                                                                                                    gccattaatgtaacaacTTTTGTTGTTTACGGCAGTTTGATTTTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACgctg                                     <  2:300748/217‑1 (MQ=255)
                                                                                                                                                                                     ccAGACAAGGAACAAACTTTGGCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACgctgc                                    <  1:197763/228‑1 (MQ=255)
                                                                                                                                                                                         aCAAGGAACAAACTTTGTCTTTACCGGTAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGTTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCCGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGAt                                <  1:129858/228‑1 (MQ=255)
                                                                                                                                                                                         aCAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGAt                                <  1:48841/228‑1 (MQ=255)
                                                                                                                                                                                                 caacttttGGCGTTAATGGGTGTTTTGTTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGTTACTTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCAt                        <  2:277439/223‑1 (MQ=255)
                                                                                                                                                                                                        tGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCa                                                                                                                    >  1:279854/1‑130 (MQ=255)
                                                                                                                                                                                                         gTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  2:279854/130‑1 (MQ=255)
                                                                                                                                                                                                              tAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTTCCGGACTGTCTGGGGACGCTGCTGGTTCTGTCATCGGggattttgt            >  1:506286/1‑226 (MQ=255)
                                                                                                                                                                                                               aaCGGCAGTTTTATTCTGGCGTTTTTTGA‑TTTGGTGGTCATTGGTTCGTCCCTGACATTTAGTC‑GGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCGTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCATCGGtgattttgat           <  1:256558/227‑1 (MQ=255)
                                                                                                                                                                                                                       tttGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑GCCAAAAGCCAGCCTTTTTAGCTGTGCAGAACCATTAAGTAAGCCGCTACCCTCTTTGCTGTTGCTGGGGG‑TAACGTTACCATTACCGGACTGGCTGGGTAAGCCGCGGGTTTTTTCAACGGGGGTTTTGTTTTAAATGGa  >  2:249348/1‑228 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  
ACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTA‑TTTGGTGGTCATTGGTACGTCCCTGACATTTAGTC‑TGTACCTGAAAGGAGCACAATT‑AATTGGC‑GG‑TCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGA‑TCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCATCGGTGATTTTGATTTCAATGGA  >  NC_000913/3838680‑3839122

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