Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,578 C→T S180S (TCG→TCA)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,5780CT93.6% 125.2 / 0.3 47S180S (TCG→TCA) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (20/24);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

AGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG‑GCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCA  >  NC_000913/4623351‑4623804
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAg                                                                                                                                                                                                             <  2:35495/251‑1 (MQ=255)
agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGGCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAg                                                                                                                                                                                                             <  2:102134/251‑1 (MQ=255)
  cagaggcaggaATGCCGGAGCCGGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGTTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTGGGCCGGGGAGATGATGTCCGCTTCGTCACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATTCCCGTCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCATTTCGCCGCCAGtt                                                                                                                                                                                                           <  2:433274/242‑1 (MQ=255)
    gATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTATGGCGTTTTGTCCGGGGAGATGATTTACGCTTCTACACGCTCGACATAATGGGCGCAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGCGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGACTTTCGCCGCCAGtttt                                                                                                                                                                                                         >  1:278045/1‑250 (MQ=255)
                  gtacccTGACCCAATTTCCACTCCCTGCGTCTCTTGCGGGCAAGGGTTTCGGCGAATTTTGGCAAGTCTGTCGTGTTGTCCGGTGAGATGATTTCCGCTTCTACACGCTCGCCATAATTGGC‑AAAAAGTGCCTTGGTTATTGCCTCGCAAACCTGCTCATTGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGTGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                                                                                                                                                                           <  2:160321/247‑1 (MQ=255)
                   gAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCtt                                                                                                                                                                                                                       >  1:101310/1‑222 (MQ=255)
                   gAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCtt                                                                                                                                                                                                                       <  2:101310/222‑1 (MQ=255)
                    aGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTg                                                                                                                                                                                          >  1:398472/1‑250 (MQ=255)
                     gccgtaCCGAAGTTCCAGTCCCGGCGACTCTGACGGGCAAAGGTGTCGGCGAAGTTTGGCAAGTCTGGCGTTTGGTCCGGGGAGAGGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCTGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTgag                                                                                                                                                                                        <  2:247402/246‑1 (MQ=255)
                       cTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTg                                                                                                                                                                                          >  1:382593/1‑247 (MQ=255)
                             gAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTg                                                                                                                                                                                >  1:119055/1‑251 (MQ=255)
                                          cggCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGTCAATCCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGGTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATTCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCtaag                                                                                                                                                                   <  2:417266/249‑4 (MQ=255)
                                            gCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGt                                                                                                                                                                 <  1:359050/251‑1 (MQ=255)
                                                       gggCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAgc                                                                                                                                                      <  1:23498/251‑2 (MQ=255)
                                                           aaaGGTTTCGGCGAAGTTGGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAgcagca                                                                                                                                                  <  2:382593/251‑1 (MQ=255)
                                                                        aaGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCgtggtg                                                                                                                                     >  1:94594/1‑251 (MQ=255)
                                                                         aGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGtggtgg                                                                                                                                    >  1:257418/1‑251 (MQ=255)
                                                                             tGGCAAGTCTGGCGTTTTTTCCGGGGAGATGATTTCCGCTTCTCCACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAATCCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAg                                                                                                                                 <  1:207973/250‑1 (MQ=255)
                                                                              ggCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGcc                                                                                                                               >  2:454594/1‑251 (MQ=255)
                                                                                cAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGCGCAGCAAGGTGCCGTGGTGGAAGccac                                                                                                                             >  2:231025/1‑249 (MQ=255)
                                                                                                 ccGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAACCCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATACCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGCGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTcgcg                                                                                                             >  2:35211/2‑250 (MQ=255)
                                                                                                    gggAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGAt                                                                                                           >  1:303828/1‑249 (MQ=255)
                                                                                                     ggAGATGATTTCCGCGTCTACACGCTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGGTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAgg                                                                                                        <  2:36253/251‑1 (MQ=255)
                                                                                                                          acgcTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGCGATAGTTTGCCAGGCGGCTAAGGGCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCggg                                                                                   >  2:422683/2‑250 (MQ=255)
                                                                                                                          acgcTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGTAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGt                                                                                   >  1:258385/2‑251 (MQ=255)
                                                                                                                             cTCGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGc                                                                                >  1:133151/1‑251 (MQ=255)
                                                                                                                              ccGCCATAATGGGCGAAAAAG‑GCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGcc                                                                               <  1:333554/250‑1 (MQ=255)
                                                                                                                                                aaG‑GCCTTGGTTATGACCTCGCAAACCTGCTCATGGGTGATCACCGGCAGCAGCTCGGGGAGGTTGGTCAAGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTc‑ccacc                                                            >  2:188172/1‑251 (MQ=255)
                                                                                                                                                        tgggtATGGCCGCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATc                                                       <  1:231025/246‑1 (MQ=255)
                                                                                                                                                               ggCCTCGCAAACCTGCTCAGGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCCCGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTAGCCGGATTGAGATAGTTTGCCAGGCTGCTAAGGGCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACgcc                                               <  2:303828/251‑3 (MQ=255)
                                                                                                                                                                 ccTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTAc                                                  <  1:474315/246‑1 (MQ=255)
                                                                                                                                                                 ccTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTAACCACCCGACCGTTAcaactc                                             >  2:474315/1‑246 (MQ=255)
                                                                                                                                                                    cGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCg                                                                                 >  1:367660/1‑212 (MQ=255)
                                                                                                                                                                    cGCAAACCTGCTCATGGGTGAGCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCg                                                                                 <  2:367660/212‑1 (MQ=255)
                                                                                                                                                                               tCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCCCCAGATCGTTACGTCCGGACGCTTCGgcact                               >  2:147779/1‑248 (MQ=255)
                                                                                                                                                                                cATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGTTAATGCCTTTCGCCGCCAGTTTCGTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCt                               <  2:457897/250‑1 (MQ=255)
                                                                                                                                                                                  tGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGAc                            >  1:42375/1‑251 (MQ=255)
                                                                                                                                                                                       gATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGccaa                       <  1:22384/251‑3 (MQ=255)
                                                                                                                                                                                              ggacgcAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAAGGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTGGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTGCGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGtt                <  2:203140/246‑1 (MQ=255)
                                                                                                                                                                                                cagcAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTc               >  1:138339/4‑250 (MQ=255)
                                                                                                                                                                                                        cGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCAtt       <  2:119055/250‑1 (MQ=255)
                                                                                                                                                                                                          gggtggtgggTCACGCGGTAACCCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTGCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTGCGCCCGGACGCGTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAg    <  2:95717/243‑1 (MQ=255)
                                                                                                                                                                                                            gatgTTGGTCACGCGGGAACCCACTGATGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTGTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGca  <  2:94594/248‑1 (MQ=255)
                                                                                                                                                                                                               tgtggggcacggggAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGc                                                                                                                                            <  1:183022/99‑1 (MQ=255)
                                                                                                                                                                                                               tgtggggcacggggAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGc                                                                                                                                            >  2:183022/12‑110 (MQ=255)
                                                                                                                                                                                                                 tGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCggc                                                                         <  1:60283/175‑3 (MQ=255)
                                                                                                                                                                                                                 tGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCggc                                                                         >  2:60283/1‑173 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
AGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG‑GCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCA  >  NC_000913/4623351‑4623804

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: ‑