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

CGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705514‑3705810
                                                                                                                                                    |                                                                                                                                                      
cGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATccc                                                                                                                                                     <  1:620518/151‑1 (MQ=255)
        tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGtt                                                                                                                                                <  2:1176368/148‑1 (MQ=255)
         cccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgcc                                                                                                                                           >  2:349486/1‑151 (MQ=255)
          ccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgccc                                                                                                                                          <  1:976249/151‑1 (MQ=255)
            aaaCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgcccgc                                                                                                                                        >  2:498309/1‑151 (MQ=255)
                cGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATg                                                                                                                                    >  1:380156/1‑151 (MQ=255)
                      gCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcatc                                                                                                                              >  2:724885/1‑151 (MQ=255)
                       cGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcatc                                                                                                                              >  2:98414/1‑150 (MQ=255)
                          cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCg                                                                                                                          >  1:1155850/1‑151 (MQ=255)
                          cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCg                                                                                                                          >  2:530298/1‑151 (MQ=255)
                          cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCg                                                                                                                          >  1:308811/1‑151 (MQ=255)
                           ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                         >  1:1132139/1‑151 (MQ=255)
                                aTTAAAATCCCTACACCCCCATACATAACCCCCCAAATGTAATGGAGATACTGCTCCCACCTCTTTTTATCTAACCCGATTCCACCCAGCCCCTGCCCGTCACCCTCGGGGGAGTCCCCACTTCCCCCCGCCATGATCATCACCGGATc                                                                                                                      <  1:125500/149‑1 (MQ=255)
                                 ttaaaaTGCCTAGATCCCCATGCAAACCCCCCAAAATGTAAGGGGGATCCTGCTGCCACATGGTTTTATCTAAGCCGGTTCCAGCAAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATc                                                                                                                      <  1:860923/145‑1 (MQ=255)
                                   gAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCg                                                                                                                 >  1:373491/1‑151 (MQ=255)
                                     aaTGCCGAGGCCGCCTGGCAAAACGCCCAAAATGTAAGGGATTTAGCGCTCCAACATGGTTTTATATAGCCGGCTTCCAGCCGCCGGCTGGTCGTGGCGCGGGGGGGAGCCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                               <  2:771307/151‑1 (MQ=255)
                                         cctatacccccATGCATACCCCCCAAAATGAAAGGGAGATCCTCCTCCCACATCGTTTTATCAAACCCGATTCCAGCCAGCAGCTGCGCGGGACGCTCTGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                            <  1:887084/141‑1 (MQ=255)
                                         cctagtcccccATGCATACCCCCCAAAATTTAAGGGAGATCCTGCTCCCACTTCGTTTTATCAAACCCGAGTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGATCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGa                                                                                                              <  1:920294/139‑1 (MQ=255)
                                         ccTAGACCGCCATGTATACCCCCCAAAATAAAAGGGGGATACGGCCGCCACCGTGTTTTATCAGAGCCGATTCCACCCCGCCGCTGCGGGTGCCGCTGGGGGGAGACCCCCGTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                            <  2:79178/150‑1 (MQ=255)
                                              acccccATGCAAAACCCCCAAAATGTAGTGGAGATCCTGCTGCCCCATCTGTTTACCAAGCCCGAGTTCAGCCAGCCGCGGCGCGGGCCGCTCGGGGGAGCCCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                          <  2:184294/143‑1 (MQ=255)
                                              aTCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTgg                                                                                                      >  1:1115107/1‑151 (MQ=255)
                                                ggCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                                    >  2:465755/2‑151 (MQ=255)
                                                cGCCTTGCATACCCCCCAAAATGTAAGGGAAATACTGCTCCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGGGTGCCGCTCGGGGGAGATCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGa                                                                                                     <  1:158358/150‑1 (MQ=255)
                                                cGCCATGCAAACCCCCCAAAATTTAATGAAGATCCTGCTCCCACATCGTTTTATCTAAGCGGAGTTCAGCCAGCGGCGGCGCGTGACGCTCGGGGGAGACCCCTCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                                    <  1:1003042/151‑1 (MQ=255)
                                                cGCCATCCATAACCCCCAAAATGTAGGGGAGATACTGCTCCCCCATGGTTTTATCTAAGCCGATTTCGGCCAGCGGGTGCGGGTGCCGCTCGGGGGAGCCCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                                    <  1:346322/151‑1 (MQ=255)
                                                    ataaaaacccccaaaaaTTTAAGGGAGATCCGGCGGCCACATCGGTTTATCTACGCCGAGTTCAGCCCCCACCGGCGCGTGCCGCTGGGGGAAGCTCCCACGTCGGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                        <  2:149375/129‑1 (MQ=255)
                                                     tGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                                               >  2:523401/1‑151 (MQ=255)
                                                        ataccccccaAAATGTAATGGAGATACTGCTGCCCCGTCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaaa                                                                                                 <  1:779857/140‑1 (MQ=255)
                                                        aTACCGCCCAAAATGTAAGGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCGGCAGCTGCGCGTGCCGCTGGGGGGAGATCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaa                                                                                            <  1:553147/151‑1 (MQ=255)
                                                            cccccaAAATTTAAGGGAGATCCTGCTGCCACATGGTTTTATCAAAGCCGAGTTCGGCCAGCAGCTGCGCGTGCCGCTGGGGGGACACCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa                                                                                           <  1:807537/146‑1 (MQ=255)
                                                            cccccaAAATGAAAGGGATATACTGCCCCCACATCGTTTTATCTAAGCCGAGTCCAGCCGGCGGCGGCGCGGGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgt                                                                                        <  1:259307/149‑1 (MQ=255)
                                                                  aaaTGTAAGGGAGATTCTGCTGCCCCATCTTTTTATCTAACCCGAGTCCGGCCGGAGGGTGGGCGTGCCGCTGGGGGGAGATCCCAGGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                                    <  2:868494/149‑1 (MQ=255)
                                                                      gaaagggggatcgtgccgcCACATGGGTTTATAAAAGCGGATTCCGCCCGGCGGCGGCGGGGGCCGTTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                                                  <  2:863099/134‑1 (MQ=255)
                                                                          gggagatctgcctcccACGTCGTTTTACCTACCCGGAGTCCGCCCACCCGCTGCGCGTCACGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  2:590021/137‑1 (MQ=255)
                                                                            gagttagcgctgcCACTTCGTTTTATCAAAGCCGAGTTCCGCCAGCAGCGGCGGGGGGCGCTGGGGGGAGACCCCACGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                               <  2:1057387/136‑1 (MQ=255)
                                                                              gaTCCTGCGCCCACATCGTTTTATCAAAGCCGAGTTCGCCCGGCAGCTGCGGGTGTCGCTGGGGGGAGACCCCTTTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                        <  1:851578/149‑1 (MQ=255)
                                                                              gaTACTGCTGCCACATCGGTTTATCTAACCCGATTTCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGa                                                                       <  1:513531/150‑1 (MQ=255)
                                                                                tccgcctcccACGTCGTTTTATCAAAGCCGAGTCCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                                     <  1:966462/142‑1 (MQ=255)
                                                                                  ctgctgCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAAc                                                                  >  2:942382/1‑151 (MQ=255)
                                                                                    gctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACg                                                                 >  2:366000/1‑150 (MQ=255)
                                                                                     ctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                               >  2:785567/1‑151 (MQ=255)
                                                                                           caTGGGTTTATCTAAGCGGTGTTCACCCGGCGGCTGCGGGTCCCGCTGGGGGGAGACCCCCGTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                                         <  1:980530/151‑1 (MQ=255)
                                                                                              cGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  1:190390/1‑151 (MQ=255)
                                                                                               ggTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  2:874318/1‑150 (MQ=255)
                                                                                                 tttATCAAAGCGGTGTTCAGCCAGCCGCGGCGCGTGCCGCTCGGGGGAGCCCCCACTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAg                                                    <  1:863053/150‑1 (MQ=255)
                                                                                                    aTCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTccc                                                >  1:673885/1‑151 (MQ=255)
                                                                                                              aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                       >  2:259118/1‑150 (MQ=255)
                                                                                                              aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                       >  1:461655/1‑150 (MQ=255)
                                                                                                                ttCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgc                                       >  2:896219/1‑147 (MQ=255)
                                                                                                                 tccagccagcGGCTGGGCGGGCCGCTCGGGGGAGACCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgag                                     <  1:221499/148‑1 (MQ=255)
                                                                                                                 tccagccagcCGCGGGGCGTGCCGCTGGGGGGATCCCCCACGTGCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                      <  2:959227/147‑1 (MQ=255)
                                                                                                                   cagccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt                                  >  2:810365/1‑151 (MQ=255)
                                                                                                                   cacccgccggcTGCGCGTGACGCGCGGGACATCCCCCACTGGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                      <  2:324173/138‑1 (MQ=255)
                                                                                                                      ccagcCGCGGGGCGTGCCGCTGGGGGGAGCCCCCCCTTCGGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATa                                 <  2:866491/149‑1 (MQ=255)
                                                                                                                       cagcggCTGCGCGTCACGCTGTGGGGACACCCAAGGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt                                  <  2:1167721/147‑1 (MQ=255)
                                                                                                                               gcggGTGCCGCTGTGGGGAGCCCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc                        <  1:610333/149‑1 (MQ=255)
                                                                                                                                 gggTGCCGCTGGGGGGAGCCCCCCCCTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACa                       <  2:997221/146‑1 (MQ=255)
                                                                                                                                          ggggggcAGCTCCCCCGTCCTCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAt             <  2:336016/146‑1 (MQ=255)
                                                                                                                                             ggggAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct        >  2:778996/1‑151 (MQ=255)
                                                                                                                                              gggATACCCCTCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCta       <  1:965703/151‑1 (MQ=255)
                                                                                                                                              gggAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct        >  1:122958/1‑150 (MQ=255)
                                                                                                                                                  gaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAg    >  2:287733/1‑150 (MQ=255)
                                                                                                                                                   aTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTg  >  1:199941/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  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: