Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,011 C→T 7.6% A259V (GCG→GTG)  tnaB → tryptophan:H(+) symporter TnaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,0110CT 7.6% 164.0 / 3.6 63A259V (GCG→GTG) tnaBtryptophan:H(+) symporter TnaB
Reads supporting (aligned to +/- strand):  ref base C (33/24);  new base T (2/3);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

TATTCCCAGCCTGATTATTTGCTATGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTT  >  NC_000913/3890868‑3891152
                                                                                                                                                |                                                                                                                                             
tatTCCCAGCCTGATTATTTGCTATGGAAAACGCAAAGATGAGTTAATCAAAAGCGTGGTATAGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGc                                                                                                                                               <  2:380324/143‑1 (MQ=255)
     ccAGCCTGAGTATTGGCTATGTAAAACGCAAAGATAAGTTACTCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGAtt                                                                                                                                           <  2:304957/143‑1 (MQ=255)
      cAGCCTGATTATTTGCTATGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTa                                                                                                                                          >  1:490899/1‑143 (MQ=255)
                 tttGCTATGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCCCTATTGCACCATGGGGAATATTCCGCGAGAAAGATTTAAGGCGCGTATCTCCTCAcgc                                                                                                                               >  1:678199/1‑143 (MQ=255)
                  ttGCTATGGAAATCGCGAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATGCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAggcg                                                                                                                              <  2:560374/143‑1 (MQ=255)
                   tGCTATGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAggcgg                                                                                                                             >  1:547272/1‑143 (MQ=255)
                   tGCTATGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAggcgg                                                                                                                             >  2:547671/1‑143 (MQ=255)
                        tGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGGGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTACGGCGATTATCTCCTCAGGCGGCatcc                                                                                                                        >  1:67356/1‑140 (MQ=255)
                          gAAAGCGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTAGTGCACCATGGGGAATATTCCGCGAGAAAGCGTTAAGGCGATTATCTCCTCAGGCGGCAACGtt                                                                                                                      <  2:796574/143‑1 (MQ=255)
                             aaCGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCGCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGAt                                                                                                                   <  1:223984/143‑1 (MQ=255)
                                gCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCGTGGGGAATATTCCGCGAGAACGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGACTCg                                                                                                                >  2:424257/1‑143 (MQ=255)
                                gCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGGATATTCCGCGAGCAAGCTTTAAGGCGATTATCTCCTCAGGCGGCCACGTTGATTCg                                                                                                                >  2:377074/1‑143 (MQ=255)
                                 cAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTCTTGCACCATGGGGAATATTCCGCGAGAAAGCTTTCAGGCGCTTATCTCCTCAGGCGCCAACGTTGATCCGc                                                                                                               >  2:409684/1‑143 (MQ=255)
                                          gTTCATCAAAGGCGGGGTATTTGGTTCGCTGCGGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCCCCATGGGCAATATTCCGCCAGATAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAAt                                                                                                      <  2:625916/143‑1 (MQ=255)
                                          gTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTCTTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGACGCAACGTTGATTAGCAAGTAAAAt                                                                                                      >  1:802225/1‑143 (MQ=255)
                                          gTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAACGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCCGGTGAAAt                                                                                                      >  2:492094/1‑143 (MQ=255)
                                          gTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAAt                                                                                                      >  2:368653/1‑143 (MQ=255)
                                           ttcaccAACAGCGTGGTATTCGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCCTGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATc                                                                                                     <  1:167029/138‑1 (MQ=255)
                                             aagcAAAAGCGTGGTATGTGGTTCGGTGCTGGCGCTGTTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCGTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGt                                                                                                   <  2:205703/140‑1 (MQ=255)
                                                cAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCCGGCGGCAACGTTGGTTCGCCGGTGAAATCGTTcc                                                                                                >  1:3692/1‑143 (MQ=255)
                                                cAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCCACGCCGCTTCGCTGGTTAACTCGTTcc                                                                                                >  1:606818/1‑143 (MQ=255)
                                                     gCGTGGTTTTTGGTTCTCTCCTGACGCTGGGGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGc                                                                                           <  2:302754/143‑1 (MQ=255)
                                                        tGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCAcc                                                                                        >  1:393556/1‑143 (MQ=255)
                                                         ggTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCa                                                                                       >  1:340077/1‑143 (MQ=255)
                                                          gTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCCACGTTGATTCGCTGGTGCAATCGTTCCTCGGCCCCaa                                                                                      >  1:769375/1‑143 (MQ=255)
                                                            aTTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCTGCAACGTTGATTCGCTGGTGAAATAGTTCATCGGAACCAAAc                                                                                    >  1:249084/1‑143 (MQ=255)
                                                            aTTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAAc                                                                                    <  2:27304/143‑1 (MQ=255)
                                                               tGGTTTGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGc                                                                                                               <  1:278734/113‑1 (MQ=255)
                                                               tGGTTTGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGc                                                                                                               >  2:278734/1‑113 (MQ=255)
                                                                       tgctgGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGTGGAATATTCCGCGAGAAAGCTTGAAGTCGATTATCTCCTCAGGCGGCAACGTGGCTTCGCTGGTGTAATTGTTGCTCGCAACACAACACCACgcgttt                                                                         >  2:442019/1‑138 (MQ=255)
                                                                       tgctgGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTGCGCGAGAAAGCTTTAAGGCGATTATCTCCTCCGGCGGCAACGTTGATTTGCTGGTGCAATAGTTCGTCGGGAGCAAACAGCACGGCAtt                                                                         >  1:251431/1‑143 (MQ=255)
                                                                       tgctgGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTACGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCTGCACCAAACAGCACGGCAtt                                                                         >  1:852540/1‑143 (MQ=255)
                                                                        gctgGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAGACAGCACGGCAttg                                                                        >  1:680630/1‑142 (MQ=255)
                                                                         ctgGCGCTGGTGATTTATCTCTTCTGTCTCTATTGCACCATGGTGAATATTCCGCGAGAGAGCTTTAAGGTGATTATCTCCTCAGGCGGCAATGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTAt                                                                       <  2:377497/143‑1 (MQ=255)
                                                                            gcgcTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATACTCCGCGAGAAAGCTTTAAGGCGATTATCCCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGa                                                                    <  1:724933/143‑1 (MQ=255)
                                                                            gcgcTGGTCATTTATCTCTTCTGGCTCTATTGCAGCCGGGGGAACATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGa                                                                    <  1:30398/143‑1 (MQ=255)
                                                                                   tGATTTACCTCCTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGc                                                             <  1:175649/143‑1 (MQ=255)
                                                                                   ggATTTAGCTCTTCTGGCTTCTGTTCACTGTGGGCAATATGCCGCGAGATAGCTTTGAGGGGATTATCTCCTCAGGTGTAAACGTTGAGTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGc                                                             <  2:352865/142‑1 (MQ=255)
                                                                                       ttATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGTGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGACACCAAACAGCACGGCATTATCGAGTTTTGCctgc                                                         >  2:865401/1‑143 (MQ=255)
                                                                                       ttATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGTGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGACACCAAACAGCACGGCATTATCGAGTTTTGCctgc                                                         >  2:861705/1‑143 (MQ=255)
                                                                                          tctcTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTgg                                                      >  1:879988/1‑143 (MQ=255)
                                                                                             cttctGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGgtgt                                                   <  1:384232/143‑1 (MQ=255)
                                                                                              ttctGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGtt                                                  >  1:478340/1‑143 (MQ=255)
                                                                                              ttctGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGCGtt                                                  >  1:665126/1‑143 (MQ=255)
                                                                                                  ggCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCCGGTGAAATCGTTCCTCGGCACCCAACCGCACGGCATTATCGAGTTTTGCCTGCTGGTGCtctcg                                              >  2:711310/1‑142 (MQ=255)
                                                                                                      tctATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAAc                                           >  1:9053/1‑143 (MQ=255)
                                                                                                       ctATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACt                                          >  1:815412/1‑143 (MQ=255)
                                                                                                         aTTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTa                                        >  1:565621/1‑143 (MQ=255)
                                                                                                         aTTGCACCATGGGGAATAGTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTa                                        <  2:864577/143‑1 (MQ=255)
                                                                                                                cATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTg                                 >  1:104523/1‑143 (MQ=255)
                                                                                                                 acggggAATATGCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGc                                <  1:526442/141‑1 (MQ=255)
                                                                                                                    gggAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTGACTTAGCTGTTGCCAg                             >  1:85839/1‑143 (MQ=255)
                                                                                                                      gcatatTCCGCGAGAGAGCTTTAAGGCGAGTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACCGCACGGCATTATCTAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAgtt                           <  2:719034/141‑1 (MQ=255)
                                                                                                                              cgcgAGACAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCttt                   <  1:586284/143‑1 (MQ=255)
                                                                                                                              cgcgAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCttt                   <  1:547671/143‑1 (MQ=255)
                                                                                                                              cgcgAGAAAGCTTTAAGGCGATTAGCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCttt                   <  1:482525/143‑1 (MQ=255)
                                                                                                                               gcgAGAAAGCGTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTg                  <  2:115988/143‑1 (MQ=255)
                                                                                                                                    aaaGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGTCTGCTGGTGTGCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTc             <  2:393556/143‑1 (MQ=255)
                                                                                                                                         tttAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGGGTCACGCt        >  1:487046/1‑143 (MQ=255)
                                                                                                                                          ttAAGGTGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGACACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTg       <  1:865401/143‑1 (MQ=255)
                                                                                                                                          ttAAGGTGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGACACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTg       <  1:861705/143‑1 (MQ=255)
                                                                                                                                           tAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTgg      <  1:672392/143‑1 (MQ=255)
                                                                                                                                            aaGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTCACTTAGCTGTTGCCAGTTCGTTCTATGGTATCACGCTggg     >  2:222152/1‑143 (MQ=255)
                                                                                                                                               gCGATTATCTCCGCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGgtt  <  1:750952/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TATTCCCAGCCTGATTATTTGCTATGGAAAACGCAAAGATAAGTTAATCAAAAGCGTGGTATTTGGTTCGCTGCTGGCGCTGGTGATTTATCTCTTCTGGCTCTATTGCACCATGGGGAATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTT  >  NC_000913/3890868‑3891152

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