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

ACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCATCGGTGATTTTGATTTCAATGGA  >  NC_000913/3838680‑3839122
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  
caccccGGTTATCCGTCGTCATTTTCGATCTCCCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCg                                                                                                                                                                                                                               <  2:415436/223‑1 (MQ=255)
       gTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAGCAACCTTTGTCGTTAACGGCAGTTTGATTCTGGCGttttt                                                                                                                                                                                                                          >  1:90655/1‑226 (MQ=255)
          attcgtcGTCATTATCGGTCCCTCCTTGCGCGTTGTTCTGGGGCCTTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGTTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTAtt                                                                                                                                                                                                                    <  2:394668/222‑1 (MQ=255)
             cgcattCATTTACGAACTCTCCGGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACGTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTtg                                                                                                                                                                                                                  <  1:424276/222‑1 (MQ=255)
                cGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGTGTATGCTGATTGGCGGGCTGATTCTGTTGCCTTTTTATGCCAGACACGGAACAAACTTTGTCGTTAACGGCAGGTTGATTTTGGCTGTTTTTTAATtggtgg                                                                                                                                                                                                              >  1:369836/1‑228 (MQ=255)
                cGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGGTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTCACGGCAGGTTGATTCTTGGGTTTTTTTATTttgtgg                                                                                                                                                                                                              >  1:575527/1‑228 (MQ=255)
                  tCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAAATTTGTCGTTAACGGCCGTTTGGTTCTGGCGTTTTTTTATTTtgtggt                                                                                                                                                                                                             >  1:208018/1‑227 (MQ=255)
                          gATCTCTCCTGCCGCCTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTAc                                                                                                                                                                                                    <  1:378881/228‑1 (MQ=255)
                             ctctcCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTTCCTTTTTATGCCCGACAAGGACCAAACATTGTCGTTAACGGCGGGTTTATTTTGGGGCTTTTTTTTTTGTGGGGCATTgtgaagcg                                                                                                                                                                                                 >  2:204961/1‑221 (MQ=255)
                              tcctcTGCCGTGTTTTTCTGGGGTATTTCCTTGGCATTTGCTGCGGCATTCTAGTCCACCTTCCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTTATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTcc                                                                                                                                                                                                <  2:369836/224‑1 (MQ=255)
                                tcCTGCCGCTTTTGTCTTGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATTCCACCTATCCCTCAACGCTAATTGCCCGGTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCt                                                                                                                                                                                              <  1:308129/228‑1 (MQ=255)
                                cccTGCCGTTTTTTTCTGGGGCATTGCCCCGGCCTTTGCTGCTGCATTCTATACCGCCTATTCCTCAACGTTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCt                                                                                                                                                                                              <  2:90655/227‑1 (MQ=255)
                                     ccGCTTTGGTCTGGGGCATTGCCTGGGCATTTGCTGCTGCATTCTATACCACCTCTCCGCCAACGATAATTGCTCGCTATGGCTCGTTACCAGTCGTCGGCGGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACAt                                                                                                                                                                                         <  2:572384/228‑1 (MQ=255)
                                     ccGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTTTCGTTAACGGCCGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTTGTACGGTCCTGACa                                                                                                                                                                                          >  1:482550/1‑227 (MQ=255)
                                          ttgttTTGGGGCATTGCCTGGGCTGTTGCGGTTTCCTTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGGCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTGCGTCCCTGACATTTAGt                                                                                                                                                                                    <  2:365437/228‑1 (MQ=255)
                                          ttgttCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGGCTTATTCTGTTGCCTTTTTATGCCAGACAAGGGACAAAGCTTTTCGGTTACGGGCGGTTTTTTTGGGCGTGTTTTTACTTTGTGGGTCATTGTACGGcccccgactttatt                                                                                                                                                                                    >  2:307096/1‑217 (MQ=255)
                                                       ttGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCCTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCCGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACACTTTGTCttgtactcaaag                                                                                                                                                                      >  1:404382/1‑217 (MQ=255)
                                                          ccTGGGCATTTGCTGCTGGCTTCTATACCCCCTGTCCCTCAACGGTAATTGCCCGGTATGGCACGTTTCCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAg                                                                                                                                                                   <  1:224254/228‑1 (MQ=255)
                                                                    tgctgctgcATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTTTTGCGTTTTTTTCTTTGTTGGTCATTGGTACATCTCTGACATTTACACTTTACCTGAAAGAAGCTCaattaat                                                                                                                                                        >  1:296076/1‑228 (MQ=255)
                                                                     gctgctgcATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTTTTGGGGTTTTTTATTTGTTGGTCCTTGGTAAGTCCCCGTACTTTAGTCTTTTACTTGAAGGAGCCCaattaatt                                                                                                                                                       >  1:567725/1‑228 (MQ=255)
                                                                                   taCCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGCCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTTGGGGTCCTTGGTACGTCCCTGAACTTTAGTCTTTACTCTAAAGGAGCACATTTAATTGTCGGTCCaaaacc                                                                                                                                       >  1:277439/1‑226 (MQ=255)
                                                                                   taCCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGAAAAGGAACAAACTTTGTCGTTAACGGCAGTTTGGTTCTGGCGTTTTTTTATTTGGTGGGCCTTGGTCAGTCCCTGGCATTTAATCCTTACCTGGAATGAGCACAATTAATTTGCGGTCCtaaaac                                                                                                                                       >  1:300748/1‑226 (MQ=255)
                                                                                     ccaccTACCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCATTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCa                                                                                                                                     <  1:500913/228‑1 (MQ=255)
                                                                                        ccTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAACATTTGTCGTTAACGGCAGTTTGCTTTTGGCGTTTTTTTATTTTGTGGTCAATGGTACGTCCCTGACCTTTTGTCTGTACCTTAAAGGAGCACAATTTATTCGCGTGCCAAAAGCCacgc                                                                                                                                  >  1:564950/1‑225 (MQ=255)
                                                                                        ccTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTAcc                                                                                                                                                                            >  1:357953/1‑189 (MQ=255)
                                                                                         cTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCt                                                                                                                                                                           <  2:357953/189‑1 (MQ=255)
                                                                                          tATCCCTGATCGCTAATTGCCCGGTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTGGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCAtt                                                                                                                                <  2:115905/228‑1 (MQ=255)
                                                                                           aTCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTTTCGTTAACGGCTGTTTTTTTTTTGCGTTTTTTTTTTTTGTGTTCATTTGTACGTTCCCTACCTTTTGTCTGTTCCTTAAAGGGGTGCCATTTATTGTCGGGTCCaaacccacacttt                                                                                                                                >  2:648343/1‑218 (MQ=255)
                                                                                           aTCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTCTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGGTTTTTTATTTGGTGGGCACTGGTACGTTCCTTGCCTTTATTCTTATACCTAAAGGGGGCAAATTCATGTGCGGATCCAAAACCAGCActt                                                                                                                               >  1:374201/1‑228 (MQ=255)
                                                                                               ttcaattttAATGGGCCGCTATGGCACGTTACCAGTCGTCGTCTGGATTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGGCAAGGAACAAGCTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAg                                                                                                                           <  2:175981/220‑1 (MQ=255)
                                                                                               cTCAATGTTTATTGCCCGTTTTGGCCCGTTACCTGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAg                                                                                                                           <  2:296076/228‑1 (MQ=255)
                                                                                               cTCAACGCTAATTGCCCGGTATGGCTCGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTTGGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAg                                                                                                                           <  1:272287/228‑1 (MQ=255)
                                                                                                tCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCCGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTGTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGc                                                                                                                          <  1:105087/228‑1 (MQ=255)
                                                                                                  aaCGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCAAGAAAAGGAACAAACTTTGGCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTTCATGGTACTCCCCTGCCCTTTAGTCTGTACCTGAAAGGAGGCCAATTAATTGGCGGTCCACAAACCAGCATTTTGAGCtg                                                                                                                        >  1:406778/1‑228 (MQ=255)
                                                                                                   aCGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGcac                                                                                                                                                                <  1:77354/190‑1 (MQ=255)
                                                                                                   aCGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCGGTTGCCTTTTTATGCCAGGCAAGGAACAAACTTTGTCGTTATCGGCAGTTTGATTCTGGCGTTTTTTTATTCGGTGGTCATTGGTACGTCCCt                                                                                                                                                                                              <  1:422262/161‑1 (MQ=255)
                                                                                                     gCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGc                                                                                                                          <  1:202514/223‑1 (MQ=255)
                                                                                                       ttaTTGCCCGCTTTGGCACGTTACCAGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  2:482550/226‑1 (MQ=255)
                                                                                                       tAATTGCCCGTTTTGGCATGTTCCCCGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  2:552625/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGGTTGGGCAAGTGAACAGGCGTCGGCTGGAGTATGCTGATTGGCGGTCGGATTATTGTGCCTTTTTATGCCAGACGAGGAACAAACTTTGGCGATAACGGCAGTTTGAGTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCCGTACCTGAAAGGAGCACAATTAAGTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:648343/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGCTTGGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:132120/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGCTTAGTCCTGTTGCCAGTCGTGGGCTGGAGTATGCCGATTGGCGGGTTGATTCGGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATCCTGGCGTTTTTTTATTTGGTGGTCATTGGTAGGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCGGTGCAg                                                                                                                   <  2:404382/228‑1 (MQ=255)
                                                                                                       tAATTGCCCGATATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:22206/228‑1 (MQ=255)
                                                                                                         actGCCCGCTATGGCCACTTTCCAGTCTTCGGCTGGATTATGCTGATTGGCCGTCTGATTCTGTTGCCTTTTTATTCCAGACAAGGAACAATCTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTAGTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGaa                                                                                                                 <  1:409381/226‑1 (MQ=255)
                                                                                                         aTTTCCCGTTATGGCACGTTACCAGTCGTCGGGTGGAGTATGCTTATTGGCGGTCTTATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACGTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGaa                                                                                                                 <  2:567725/228‑1 (MQ=255)
                                                                                                                gttttggatcGTTCCCAGTCGTCGGGTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATTGCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTaa                                                                                                          <  2:494941/219‑1 (MQ=255)
                                                                                                                          gTTTCCAGTTGTCGGCTGGGGTATGCTTATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCt                                                                                                <  2:487179/228‑1 (MQ=255)
                                                                                                                            tACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTACTTGGGAATCCCCTGCATTTTATTCTTTACCTGAAAGGGTGCCAATTAAATTGGGGGTCCAAAACCCGGCTTTTTAGCTTTGCGAAAACAATAAGAAgagcgcaaa                                                                                              >  2:150447/1‑225 (MQ=255)
                                                                                                                             tcccgttgtcgGCTGGGTTTGGCTGGTTGGGCGTCTGATGTCGTTGCCGTTTTCTGCCAGAGCAGGACCGAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGTATTTTGAGCTGTGCAGAACCATTAAGTAGCGGGCTAct                                                                                             <  2:124269/221‑1 (MQ=255)
                                                                                                                             gccAGTCGTCGGCTGGAGGATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAATGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTAc                                                                                              <  1:274961/226‑1 (MQ=255)
                                                                                                                                  tcgtcgGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACGTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  1:345341/201‑1 (MQ=255)
                                                                                                                                    gtcgGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCCTGCATTTTATCCGTACCTTTAAGGGGGCAAATTCATTTGCGGGCCCACAACCCGCATTTTTGCCTTTGCAGAACCCTTTAGAAacgcgccaaaccctttt                                                                                      >  2:122296/1‑217 (MQ=255)
                                                                                                                                      cgGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGTTTCTGGCGTTTTTTTTTTTTGGTGTCATTTGGACGTTCCTGACATTTTTGCTGTATCTGAAATGAGCCCAATTAATTTGCGGTCCCAAAGCCCGCCTTTTTGGCCGGGCAGAAACCTTAAAGAGAGCGCTACCCCCtttttt                                                                                    >  2:48841/1‑225 (MQ=255)
                                                                                                                                       tgtttgAGTATCGTGATGGGCGGTCTGTTTCTGGTTGTTTTTTATGCCGTACAAGGAACAAACTTTGGCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTTAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTttgctg                                                                                   <  2:51460/223‑1 (MQ=255)
                                                                                                                                             aGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGaa                                                                                                                 <  1:115655/192‑1 (MQ=255)
                                                                                                                                                  gCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTcc                                                                                                                                                                                                <  1:168462/112‑1 (MQ=255)
                                                                                                                                                           gcgTCTTATTCTTTTGCCTTTTTATGCCAGACAAGGAAAAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTca                                                               <  1:204961/225‑1 (MQ=255)
                                                                                                                                                             gTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGGTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTTGCGGGCCAAAAACCAGCATTTTTAGCTGTGACGAACCATTAAGTAGCGCGCTACTCTCCTTTCTGTTGTTGGGGATCaaggtccca                                                             >  1:439758/1‑219 (MQ=255)
                                                                                                                                                               ttgtttCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTa                                                          <  2:204990/224‑1 (MQ=255)
                                                                                                                                                                ttaTTCTGTTGCCTTTTTATGCCAGTCAAAGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCCCGTTCACAttcc                                                         <  1:149311/226‑3 (MQ=255)
                                                                                                                                                                    tCTGGTGCTTTTTTAGGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGCCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGa                                                     <  2:208018/228‑1 (MQ=255)
                                                                                                                                                                        ttGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTAAGTCCCTGCCATTTAGTCTGTACCTGAAAGGAGCCCAATTAATTGGCGGTCCCAAAGCCATCATTTTGAGCTTGGAAGAACCATTAAGTAGCGCGCTACTCCCTTTGCGTTTGCTGGGGAGTCAGTTCACAGTAACCGAccag                                                 >  2:199436/1‑225 (MQ=255)
                                                                                                                                                                             tttttATTCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTTTTGCTTTTTTTTATTTGGTGGTCATTGGTGCGTCCCGGCCATTTAGTCTGGACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGTATCACGTTCACATTACCGGACTGGCTggg                                            <  2:406778/228‑1 (MQ=255)
                                                                                                                                                                                    gccattaatgtaacaacTTTTGTTGTTTACGGCAGTTTGATTTTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACgctg                                     <  2:300748/217‑1 (MQ=255)
                                                                                                                                                                                     ccAGACAAGGAACAAACTTTGGCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACgctgc                                    <  1:197763/228‑1 (MQ=255)
                                                                                                                                                                                         aCAAGGAACAAACTTTGTCTTTACCGGTAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGTTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCCGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGAt                                <  1:129858/228‑1 (MQ=255)
                                                                                                                                                                                         aCAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGAt                                <  1:48841/228‑1 (MQ=255)
                                                                                                                                                                                                 caacttttGGCGTTAATGGGTGTTTTGTTCTGGCGTTTTTTTATTTGGTGGTCATTGTTACTTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCAt                        <  2:277439/223‑1 (MQ=255)
                                                                                                                                                                                                        tGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCa                                                                                                                    >  1:279854/1‑130 (MQ=255)
                                                                                                                                                                                                         gTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAg                                                                                                                   <  2:279854/130‑1 (MQ=255)
                                                                                                                                                                                                              tAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTTCCGGACTGTCTGGGGACGCTGCTGGTTCTGTCATCGGggattttgt            >  1:506286/1‑226 (MQ=255)
                                                                                                                                                                                                               aaCGGCAGTTTTATTCTGGCGTTTTTTGATTTGGTGGTCATTGGTTCGTCCCTGACATTTAGTCGGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCGTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCATCGGtgattttgat           <  1:256558/227‑1 (MQ=255)
                                                                                                                                                                                                                       tttGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGGCCAAAAGCCAGCCTTTTTAGCTGTGCAGAACCATTAAGTAAGCCGCTACCCTCTTTGCTGTTGCTGGGGGTAACGTTACCATTACCGGACTGGCTGGGTAAGCCGCGGGTTTTTTCAACGGGGGTTTTGTTTTAAATGGa  >  2:249348/1‑228 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                  
ACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAGTATGCTGATTGGCGGTCTGATTCTGTTGCCTTTTTATGCCAGACAAGGAACAAACTTTGTCGTTAACGGCAGTTTGATTCTGGCGTTTTTTTATTTGGTGGTCATTGGTACGTCCCTGACATTTAGTCTGTACCTGAAAGGAGCACAATTAATTGGCGGTCCAAAAGCCAGCATTTTGAGCTGTGCAGAACCATTAAGTAGCGCGCTACTCTCTTTGCTGTTGCTGGGGATCACGTTCACATTACCGGACTGGCTGGGAACGCTGCTGATTCTGTCATCGGTGATTTTGATTTCAATGGA  >  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: