Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6610TC41.0% 81.4 / 27.5 63I44V (ATC→GTC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (27/9);  new base C (0/25);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 7.09e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705514‑3705810
                                                                                                                                                    |                                                                                                                                                      
cGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATccc                                                                                                                                                     <  1:620518/151‑1 (MQ=255)
        tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGtt                                                                                                                                                <  2:1176368/148‑1 (MQ=255)
         cccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑cgcc                                                                                                                                           >  2:349486/1‑151 (MQ=255)
          ccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑cgccc                                                                                                                                          <  1:976249/151‑1 (MQ=255)
            aaaCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑cgcccgc                                                                                                                                        >  2:498309/1‑151 (MQ=255)
                cGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATg                                                                                                                                    >  1:380156/1‑151 (MQ=255)
                      gCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGatcatc                                                                                                                              >  2:724885/1‑151 (MQ=255)
                       cGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGatcatc                                                                                                                              >  2:98414/1‑150 (MQ=255)
                          cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCg                                                                                                                          >  1:1155850/1‑151 (MQ=255)
                          cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCg                                                                                                                          >  2:530298/1‑151 (MQ=255)
                          cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCg                                                                                                                          >  1:308811/1‑151 (MQ=255)
                           ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCgg                                                                                                                         >  1:1132139/1‑151 (MQ=255)
                                aTTAAAATCCCTACACCCCCATACATAACCCCCCAAATGTAATGGAGATACTGCTCCCACCTCTTTTTATCTAACCCGATTC‑CACCCAGCCCCTGCCCGTCACCCTCGGGGGAGTCCCCACTTC‑CCCCCGCCATGATCATCACCGGATc                                                                                                                      <  1:125500/149‑1 (MQ=255)
                                 ttaaaaTGCCTAGATCCCCATGCAAACCCCCCAAAATGTAAGGGGGATCCTGCTGCCACATGGTTTTATCTAAGCCGGTTC‑CAGCAAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTC‑CGCCCGCCATGATCATCACCGGATc                                                                                                                      <  1:860923/145‑1 (MQ=255)
                                   gAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCg                                                                                                                 >  1:373491/1‑151 (MQ=255)
                                     aaTGCCGAGGCCGCCTGGCAAAACGCCCAAAATGTAAGGGATTTAGCGCTCCAACATGGTTTTATAT‑AGCCGGCTTCCAGCCGCCGGCTGGTCGTGGCGCGGGGGGGAGCCCCCACGTC‑CGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                               <  2:771307/151‑1 (MQ=255)
                                         cctatacccccATGCATACCCCCCAAAATGAAAGGGAGATCCTCCTCCCACATCGTTTTATCAAACCCGATTC‑CAGCCAGCAGCTGCGCGGGACGCTCTGGGGAGACCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                            <  1:887084/141‑1 (MQ=255)
                                         cctagtcccccATGCATACCCCCCAAAATTTAAGGGAGATCCTGCTCCCACTTCGTTTTATCAAACCCGAGTT‑CAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGATCCCACTTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGa                                                                                                              <  1:920294/139‑1 (MQ=255)
                                         ccTAGACCGCCATGTATACCCCCCAAAATAAAAGGGGGATACGGCCGCCACCGTGTTTTATCAGAGCCGATTC‑CACCCCGCCGCTGCGGGTGCCGCTGGGGGGAGA‑CCCCCGTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                            <  2:79178/150‑1 (MQ=255)
                                              acccccATGCAAAACCCCCAAAATGTAGTGGAGATCCTGCTGCCCCATCTGTTTACCAAGCCCGAGTT‑CAGCCAGCCGCGGCGCGGGCCGCTCGGGGGAGCCCCCGCGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                          <  2:184294/143‑1 (MQ=255)
                                              aTCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTgg                                                                                                      >  1:1115107/1‑151 (MQ=255)
                                                ggCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                                    >  2:465755/2‑151 (MQ=255)
                                                cGCCTTGCATACCCCCCAAAATGTAAGGGAAATACTGCTCCCACATCGTTTTATCTAAGCCGAGTT‑CAGCCAGCAGCGGCGGGTGCCGCTCGGGGGAGATCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGa                                                                                                     <  1:158358/150‑1 (MQ=255)
                                                cGCCATGCAAACCCCCCAAAATTTAATGAAGATCCTGCTCCCACATCGTTTTATCTAAGCGGAGTT‑CAGCCAGCGGCGGCGCGTGACGCTCGGGGGAGACCCCTCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                                    <  1:1003042/151‑1 (MQ=255)
                                                cGCCATCCATAACCCCCAAAATGTAGGGGAGATACTGCTCCCCCATGGTTTTATCTAAGCCGATTT‑CGGCCAGCGGGTGCGGGTGCCGCTCGGGGGAGCCCCCCCTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                                    <  1:346322/151‑1 (MQ=255)
                                                    ataaaaacccccaaaaaTTTAAGGGAGATCCGGCGGCCACATCGGTTTATCTACGCCGAGTT‑CAGCCCCCACCGGCGCGTGCCGCTGGGGGAAGCTCCCACGTC‑GGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                        <  2:149375/129‑1 (MQ=255)
                                                     tGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                                               >  2:523401/1‑151 (MQ=255)
                                                        ataccccccaAAATGTAATGGAGATACTGCTGCCCCGTCGTTTTATCTAAGCCGAGTT‑CAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaaa                                                                                                 <  1:779857/140‑1 (MQ=255)
                                                        aTACCGCCCAAAATGTAAGGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTT‑CAGCCGGCAGCTGCGCGTGCCGCTGGGGGGAGATCCCACTTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaa                                                                                            <  1:553147/151‑1 (MQ=255)
                                                            cccccaAAATTTAAGGGAGATCCTGCTGCCACATGGTTTTATCAAAGCCGAGTT‑CGGCCAGCAGCTGCGCGTGCCGCTGGGGGGACACCCCGCGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa                                                                                           <  1:807537/146‑1 (MQ=255)
                                                            cccccaAAATGAAAGGGATATACTGCCCCCACATCGTTTTATCTAAGCCGAGTC‑CAGCCGGCGGCGGCGCGGGACGCTCGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgt                                                                                        <  1:259307/149‑1 (MQ=255)
                                                                  aaaTGTAAGGGAGATTCTGCTGCCCCATCTTTTTATCTAACCCGAGTC‑CGGCCGGAGGGTGGGCGTGCCGCTGGGGGGAGATCCCAGGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                                    <  2:868494/149‑1 (MQ=255)
                                                                      gaaagggggatcgtgccgcCACATGGGTTTATAAAAGCGGATTC‑CGCCCGGCGGCGGCGGGGGCCGTTGGGGGGAGACCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                                                  <  2:863099/134‑1 (MQ=255)
                                                                          gggagatctgcctcccACGTCGTTTTACCTACCCGGAGTC‑CGCCCACCCGCTGCGCGTCACGCTGGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  2:590021/137‑1 (MQ=255)
                                                                            gagttagcgctgcCACTTCGTTTTATCAAAGCCGAGTT‑CCGCCAGCAGCGGCGGGGGGCGCTGGGGGGAGACCCCACGTG‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                               <  2:1057387/136‑1 (MQ=255)
                                                                              gaTCCTGCGCCCACATCGTTTTATCAAAGCCGAGTT‑CGCCCGGCAGCTGCGGGTGTCGCTGGGGGGAGACCCCTTTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                        <  1:851578/149‑1 (MQ=255)
                                                                              gaTACTGCTGCCACATCGGTTTATCTAACCCGATTT‑CAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGa                                                                       <  1:513531/150‑1 (MQ=255)
                                                                                tccgcctcccACGTCGTTTTATCAAAGCCGAGTC‑CAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                                     <  1:966462/142‑1 (MQ=255)
                                                                                  ctgctgCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAAc                                                                  >  2:942382/1‑151 (MQ=255)
                                                                                    gctgcCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACg                                                                 >  2:366000/1‑150 (MQ=255)
                                                                                     ctgcCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                               >  2:785567/1‑151 (MQ=255)
                                                                                           caTGGGTTTATCTAAGCGGTGTT‑CACCCGGCGGCTGCGGGTCCCGCTGGGGGGAGA‑CCCCCGTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                                         <  1:980530/151‑1 (MQ=255)
                                                                                              cGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  1:190390/1‑151 (MQ=255)
                                                                                               ggTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  2:874318/1‑150 (MQ=255)
                                                                                                 tttATCAAAGCGGTGTT‑CAGCCAGCCGCGGCGCGTGCCGCTCGGGGGAGCCCCCACTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAg                                                    <  1:863053/150‑1 (MQ=255)
                                                                                                    aTCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTccc                                                >  1:673885/1‑151 (MQ=255)
                                                                                                              aGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                       >  2:259118/1‑150 (MQ=255)
                                                                                                              aGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                       >  1:461655/1‑150 (MQ=255)
                                                                                                                tt‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgc                                       >  2:896219/1‑147 (MQ=255)
                                                                                                                 tccagccagcGGCTGGGCGGGCCGCTCGGGGGAGACCCCCCTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgag                                     <  1:221499/148‑1 (MQ=255)
                                                                                                                 tccagccagcCGCGGGGCGTGCCGCTGGGGGGATCCCCCACGTG‑CCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                      <  2:959227/147‑1 (MQ=255)
                                                                                                                   cagccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt                                  >  2:810365/1‑151 (MQ=255)
                                                                                                                   cacccgccggcTGCGCGTGACGCGCGGGACATCCCCCACTGG‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                      <  2:324173/138‑1 (MQ=255)
                                                                                                                      ccagcCGCGGGGCGTGCCGCTGGGGGGAGCCCCCCCTTC‑GGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATa                                 <  2:866491/149‑1 (MQ=255)
                                                                                                                       cagcggCTGCGCGTCACGCTGTGGGGACACCCAAGGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt                                  <  2:1167721/147‑1 (MQ=255)
                                                                                                                               gcggGTGCCGCTGTGGGGAGCCCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc                        <  1:610333/149‑1 (MQ=255)
                                                                                                                                 gggTGCCGCTGGGGGGAGCCCCCCCCTG‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACa                       <  2:997221/146‑1 (MQ=255)
                                                                                                                                          ggggggcAGCTCCCCCGTC‑CTCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAt             <  2:336016/146‑1 (MQ=255)
                                                                                                                                             ggggAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct        >  2:778996/1‑151 (MQ=255)
                                                                                                                                              gggATACCCCTCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCta       <  1:965703/151‑1 (MQ=255)
                                                                                                                                              gggAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct        >  1:122958/1‑150 (MQ=255)
                                                                                                                                                  gaTCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAg    >  2:287733/1‑150 (MQ=255)
                                                                                                                                                   aTCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTg  >  1:199941/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705514‑3705810

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