Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6610TC43.3% 86.7 / 12.5 60I44V (ATC→GTC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (28/6);  new base C (1/25);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 3.19e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTAT  >  NC_000913/3705511‑3705806
                                                                                                                                                      |                                                                                                                                                 
gCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGAt                                                                                                                                                   <  2:826348/151‑1 (MQ=255)
 cACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATc                                                                                                                                                  <  2:181084/151‑1 (MQ=255)
          ttCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcg                                                                                                                                         >  2:343561/1‑151 (MQ=255)
          ttCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcg                                                                                                                                         >  2:235047/1‑151 (MQ=255)
           tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcg                                                                                                                                         <  1:448424/150‑1 (MQ=255)
            cccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgcc                                                                                                                                       >  2:597382/1‑151 (MQ=255)
                        tGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcat                                                                                                                           >  1:829766/1‑151 (MQ=255)
                          cGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcatc                                                                                                                          >  1:116373/1‑150 (MQ=255)
                               tttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                     >  1:357882/1‑150 (MQ=255)
                                 tCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATc                                                                                                                  >  2:952154/1‑151 (MQ=255)
                                  aaTTGAAATGCCTAGACCGCCATGCAAACCGCCCAAAATGAAAGGGAGATCCTGCTCCCACATCGGTTTATCTAAGCCGAGTCCAGCCAGCAGCGGCGCGGGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCg                                                                                                                 <  1:101740/150‑1 (MQ=255)
                                    ttaaaaTCCCTAGACCCCCATGAATACCCCCCAAAATGTAATGGAGATCCTGCGCCCACACCGTTTTATCAAAGCCGATTTCACCCAGCAGCGGCGCGTGCCGCTGGGGGGAGACCCCCCGTTCCCCCGCCATGATCATCACCGGATCGc                                                                                                                <  2:757509/147‑1 (MQ=255)
                                      aaaaTCCCTGGACCGCAATGCAAACCCCCCAAAATTTAAGGAAGATACTGCCGCCACATCGGTTTATCTAACCCGATTCCAGCCAGCAGCGGCGCGTGCCGCGCTGGGAAGATCCCACGTGCGCCCGCCATGATCATCACCGGATCGCCCg                                                                                                             <  2:326057/150‑1 (MQ=255)
                                       aaaTTCCTAGACCCCCTTGCTTAACGCCTAAATTGTAAGGGAAATCCTCCTGCCACATCGGTTTATATAAGCCGATTCCAGCCAGCCGCTGCGCGTGACGTTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCgg                                                                                                            <  2:279456/151‑1 (MQ=255)
                                        aaTGCCTAGATCGCCATGCATAACCCCCAAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                           <  2:140844/151‑1 (MQ=255)
                                        aaTCCCTAGCCCCCCATGCATACCCCCCAAATTGTAAGGGAGATCCTGCTCCCACATCTTTTTATCTAACCCGTTTTCACCCAGCGGCTGCGCTTGCCGCTGTGGGGAGCCCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                           <  1:1011772/151‑1 (MQ=255)
                                         aTGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                           >  2:925356/1‑150 (MQ=255)
                                           gCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                        >  2:612604/1‑151 (MQ=255)
                                              tAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtg                                                                                                     >  2:413730/1‑151 (MQ=255)
                                                      cATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa                                                                                              >  2:972467/1‑150 (MQ=255)
                                                           aTAACCCCCAAAATGTAATGGAGATCCTGCTGCCACACCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCa                                                                                         <  2:330186/150‑1 (MQ=255)
                                                            tAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaa                                                                                        >  2:482462/1‑150 (MQ=255)
                                                             accGCCCAAAATATAAGGGGGATTCTGCCGCCACACCGTTTTATCTACCCGGATTAAGGCCAGCAGCTGCGCGGGCCGCTCGGGGGAGCCCCCCGTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCa                                                                                         <  1:310754/146‑1 (MQ=255)
                                                             aaCGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAt                                                                                      >  2:456870/1‑151 (MQ=255)
                                                             aaCGCCCAAAATGTAAGGGAGTTACTGCTCCCACATCGTTTTATCAAGGCCGAGTTCACCCGGCAGCTGCGCGGGCCGCTGTGGGGATACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCa                                                                                         <  2:878630/148‑1 (MQ=255)
                                                               cccccaAAATGTAAGGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                     <  2:220413/148‑1 (MQ=255)
                                                                        tGTAATGGAGATACTGCCCCCCCATCGGTTTATCAAGGCCGTGTTCCGCCGGGCGCGGCCCGGGCCGCTGGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                            <  1:885803/150‑1 (MQ=255)
                                                                           aagggAGATCCTCCTCCCCCATCGTTTTACCTAAGCCGAGTTCAGCCAGCAGCTGGGCGTGCCGCTGGGGGGAGACCCCCCGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                          <  2:29415/146‑1 (MQ=255)
                                                                              ggAGATCCTGCTCCCACATGGTTTTATCTAAGCCGATTCCAGCCAGCGGCGGCGCGGGCCGCTGGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAcc                                                                     <  1:413709/151‑1 (MQ=255)
                                                                                agaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                    >  1:176451/1‑150 (MQ=255)
                                                                                  aTACTGCCGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCCCGTGACGCTCTGGCGAGATCCCACGTTCGCCCGCCATGATCATCACCGCCTCGCCCGGGCTCATGTCGACAAACGCAATTGCGAGAAGGGTAATACCGAt                                                                  >  1:597868/1‑150 (MQ=255)
                                                                                     ctgctcCCACATGGGTTTATAAAAGCCGATTACGGCCAGCAGCGGCGGGGGCCGCTGGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                <  1:72837/149‑1 (MQ=255)
                                                                                       gctccCACATCGTTTTATCAAAGCCGAGTCCAGCCAGCAGCGGCGCGTGACGTTGGGGGGAGACCCCCCGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                <  2:253023/147‑1 (MQ=255)
                                                                                        ctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                           >  2:350261/1‑151 (MQ=255)
                                                                                           ccACATCGTTTTATCAAAGCCGATTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                        <  2:105383/151‑1 (MQ=255)
                                                                                             acaTCGGTTTATCTAAGCCGGTTCCACCCAGCAGCCGGGCGGGCGGCTGGGGGGGGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACg                                                             <  1:840473/144‑1 (MQ=255)
                                                                                             acaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                      >  2:673549/1‑151 (MQ=255)
                                                                                                  gttttATCAAGGCCGAGTCCAGCCAGCAGCTGCGCGGGCCGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                  <  1:977766/148‑1 (MQ=255)
                                                                                                       aTCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTccc                                            >  2:919265/1‑151 (MQ=255)
                                                                                                           aaGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                        >  1:282811/1‑151 (MQ=255)
                                                                                                            aGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACg                                       >  1:709177/1‑151 (MQ=255)
                                                                                                                 aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                   >  1:941699/1‑150 (MQ=255)
                                                                                                                    tCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgaga                                >  2:541242/1‑150 (MQ=255)
                                                                                                                    ccacccagcGGCGGCGCGTGCCGCGGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgaga                                <  1:685542/146‑1 (MQ=255)
                                                                                                                      agccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGCGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt                              >  2:305008/1‑150 (MQ=255)
                                                                                                                              gcggcgcGTGGCGCTGTGGGGAGATCCCAGGGCCGACCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCa                      <  1:776067/147‑1 (MQ=255)
                                                                                                                                 gcgcGGGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCaa                     <  1:22800/148‑1 (MQ=255)
                                                                                                                                  cgcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACAt                  >  1:29699/1‑150 (MQ=255)
                                                                                                                                   gtgggccgcgggggggAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACa                   <  1:625747/136‑1 (MQ=255)
                                                                                                                                       gccgccggggggAGTTCCCCCGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc             <  2:572246/142‑1 (MQ=255)
                                                                                                                                        aCGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCgg           >  2:302231/1‑151 (MQ=255)
                                                                                                                                          gcgcgggggAGCCCCCACGTCCGCCCGCCATGATCATCACCGAATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc             <  1:886993/142‑1 (MQ=255)
                                                                                                                                           ctctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAtt        >  1:243051/1‑151 (MQ=255)
                                                                                                                                             ggggggAGACCCATCGTGCAACCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtct      <  1:631792/149‑1 (MQ=255)
                                                                                                                                              gggggAGACCCCTCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGa          <  2:68843/145‑1 (MQ=255)
                                                                                                                                              gggggAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtctc     <  1:371945/150‑1 (MQ=255)
                                                                                                                                               ggggAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct    >  2:242889/1‑151 (MQ=255)
                                                                                                                                               ggggAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGa          <  1:931117/145‑1 (MQ=255)
                                                                                                                                                gggAGCCCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCta   <  1:53580/151‑1 (MQ=255)
                                                                                                                                                   agaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtat  >  2:500700/1‑149 (MQ=255)
                                                                                                                                                   agaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtat  >  2:468222/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTAT  >  NC_000913/3705511‑3705806

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