Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,0900AG40.4% 83.3 / 10.8 57I83V (ATT→GTT) hisFimidazole glycerol phosphate synthase, catalytic subunit with HisH
Reads supporting (aligned to +/- strand):  ref base A (6/28);  new base G (23/0);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.08e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  NC_000913/2095945‑2096240
                                                                                                                                                 |                                                                                                                                                         
gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAg                                                                                                                                                      <  1:212948/151‑1 (MQ=255)
gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGtt                                                                                                                                                         >  2:184197/1‑148 (MQ=255)
  tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaa                                                                                                                                                       >  2:31072/1‑148 (MQ=255)
  tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGtt                                                                                                                                                    >  1:534443/1‑150 (MQ=255)
      aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTaa                                                                                                                                                       >  2:51437/1‑144 (MQ=255)
      aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTg                                                                                                                                               <  2:656673/151‑1 (MQ=255)
      aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCATGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGtc‑t                                                                                                                                                  >  1:618330/1‑148 (MQ=255)
       aaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGGT‑TCTgg                                                                                                                                              >  1:652708/1‑151 (MQ=255)
         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGaa                                                                                                                                            <  2:302747/151‑1 (MQ=255)
         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGaa                                                                                                                                            <  1:941645/151‑1 (MQ=255)
                gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCcg                                                                                                                                     <  2:557767/151‑1 (MQ=255)
                gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTGAATT‑TCTggga                                                                                                                                            >  1:643029/1‑142 (MQ=255)
                   gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCcgc‑g                                                                                                                                  <  1:275242/150‑1 (MQ=255)
                   gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCcgc‑g                                                                                                                                  <  1:570396/150‑1 (MQ=255)
                          gcgcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTGAGTC‑ACTGGGAAAAGGCGC‑CA‑AAAATCt                                                                                                                         >  1:522547/1‑151 (MQ=255)
                            gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTATGTC‑TCTGGGAGATACCGC‑GC‑AAAtttttt                                                                                                                       >  1:678391/1‑151 (MQ=255)
                            gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTC‑TCTTGGAGAAGGCCC‑GA‑AAAtttttt                                                                                                                       >  2:860869/1‑151 (MQ=255)
                             cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAAGC‑TCTGGAAGATGCCGC‑GA‑AAAtttttt                                                                                                                       >  2:794392/1‑150 (MQ=255)
                             cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCttt                                                                                                                       >  2:626832/1‑150 (MQ=255)
                                   aaCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTgg                                                                                                                <  1:278820/151‑1 (MQ=255)
                                    aCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAAGC‑TCTTGGAGATGTCGT‑GA‑AAACTCttttcattcg                                                                                                                >  1:505203/1‑143 (MQ=255)
                                     cTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAATT‑TCTGGGAGAATGCGC‑GA‑AAATtcttttcttttg                                                                                                                >  1:678235/1‑147 (MQ=255)
                                        gtgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCttt                                                                                                                  <  2:754449/144‑1 (MQ=255)
                                        gtgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCgg                                                                                                           <  1:451640/151‑1 (MQ=255)
                                        gtgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGATTAGTCTTTTGGAAG‑TGGCGCGGA‑AAATTTTTTCCTTTgggggg                                                                                                            >  1:834125/1‑147 (MQ=255)
                                              tACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGTC‑TCTGGGAGAAGGCGC‑GA‑AAATTTTTTTCTTTGGCGCGGGTaaa                                                                                                      >  1:690737/1‑150 (MQ=255)
                                               aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATaa                                                                                                       <  2:790336/148‑1 (MQ=255)
                                                 tatatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTC‑TCTGGGAGAAGCCGC‑GA‑AAAATCTTTTCTTTGGCGCGGGTAAAAttt                                                                                                  >  2:504456/2‑151 (MQ=255)
                                                 gagatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAAttt                                                                                                  <  1:609697/148‑1 (MQ=255)
                                                  atatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAAt                                                                                                    <  1:881236/148‑1 (MQ=255)
                                                  atatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTC‑TCTGGAAGAAGCCGCGGA‑AAATTTTTTCTTTTgggcggggaaaaatt                                                                                                   >  1:499468/1‑137 (MQ=255)
                                                    atCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTcc                                                                                                <  2:625827/150‑1 (MQ=255)
                                                    atCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTC‑TCTGGGAGATGCCGC‑CC‑AAATTTTTTTCTTTTGCGCGGATAAAATTTcc                                                                                                >  2:345505/1‑150 (MQ=255)
                                                     tCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAAATTTTTTCTTTTGCGCGGGTAAAAATTCCAt                                                                                              >  2:202339/1‑151 (MQ=255)
                                                            ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATC‑TCTGGAAGAAGCCCC‑GA‑AAATTCTTTTCTTTGGCGCGGGTAAAAATTTCATCaa                                                                                           >  2:85238/1‑147 (MQ=255)
                                                             tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAActc                                                                                        <  1:42539/149‑1 (MQ=255)
                                                                 gCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgc                                                                                  <  2:185949/151‑1 (MQ=255)
                                                                                     aaaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACatta                                                               <  1:324441/150‑1 (MQ=255)
                                                                                     aaaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGGC‑TCTGGGAGATGCCGC‑GC‑AAATTCTTTCCTTTTGCGCGGATAAAATTTCCATCAACCCTCCTGCGCTGGCAGACCCAACatt                                                                >  2:281719/1‑149 (MQ=255)
                                                                                             ggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTC‑TCTGGAAGATGCCCC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAACCAACATTAATTACTCGc                                                      >  2:217968/1‑151 (MQ=255)
                                                                                                 tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTGATAT‑TCTGGGGGATGGCCC‑GA‑AAATTCTTTTCTTTTGCGCGGATAAAATTTTCTTTAAGTCTCCTGCGCTGGGAAACCCAACAATTATTTCTTTCCTg                                                   >  2:920635/1‑150 (MQ=255)
                                                                                                    cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCTCCTGGc                                                 <  1:920538/149‑1 (MQ=255)
                                                                                                    cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGcc                                                <  2:805614/150‑1 (MQ=255)
                                                                                                      cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGATAAGGC‑TCTTGGAGAAGGCGC‑GG‑AAAATTTTTTCTTTGGGGCGGGGAAAAATTTCATCAACACTCCCGCGCCGGCAGACCCAACAATAATTACTAGCCCGGcc                                                >  1:289410/1‑148 (MQ=255)
                                                                                                          gCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCt                                         <  2:423728/151‑1 (MQ=255)
                                                                                                           cGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGC‑TCTGGAAGATGCCCC‑GA‑AAAATTTTTTCTTTTGCGCCGGTAAAATTTTCCTCAACTCTCCTGCGCCGGCAGAACCCACCTTAATTACTCGCCTTGCCGCTc                                            >  1:714735/1‑147 (MQ=255)
                                                                                                               aGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACGCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGc                                    <  1:477862/151‑1 (MQ=255)
                                                                                                                gTGATCGACATTCCGTTTTGTGTGGCGGGCGGGGTTTAGTG‑TCTGGCACATGCCCC‑GACAAATCTTTTCCTTTGG‑GTGGGTAAAATTTTCATCAACCCTCCCGCGCTGGCCGACACACCACTAATTACTCTCCCGGCCGATTGCctt                                       >  1:709510/1‑147 (MQ=255)
                                                                                                                   aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAGGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTg                                 <  1:857552/150‑1 (MQ=255)
                                                                                                                      gACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgt                             <  1:535145/151‑1 (MQ=255)
                                                                                                                       aCATTCCGTTTTGTGTGGCGGGGGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTTTTTTCTTTGGCGCGGGTAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtg                            >  2:143468/1‑151 (MQ=255)
                                                                                                                        cATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtgt                           <  2:894936/151‑1 (MQ=255)
                                                                                                                        cATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTC‑TCTGGGAGAAGCCGC‑GG‑AAAATTTTTTCTTTGGGGCGGGTAAAATTTTCATCAAATCTCCTGCGCGGGCAGAACCCACATTAATTAATCGCCTGGGCGATAGCTTTTGCGTGGAg                              >  1:665852/1‑148 (MQ=255)
                                                                                                                         aTTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTa                          <  2:912924/151‑1 (MQ=255)
                                                                                                                             cGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt                      <  1:551910/151‑1 (MQ=255)
                                                                                                                                     gtgGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGt               <  1:374337/150‑1 (MQ=255)
                                                                                                                                           ggTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATa        <  1:672161/151‑1 (MQ=255)
                                                                                                                                                 aTTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt  <  1:862453/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                         
GCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTC‑TCTGGAAGATGCCGC‑GA‑AAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  NC_000913/2095945‑2096240

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