Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6610TC43.5% 215.4 / 42.6 138I44V (ATC→GTC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (67/11);  new base C (0/60);  total (67/71)
Fisher's exact test for biased strand distribution p-value = 4.12e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

GCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705511‑3705810
                                                                                                                                                          |                                                                                                                                                     
gCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGAt                                                                                                                                                       >  3:101103/1‑151 (MQ=255)
 cACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATc                                                                                                                                                      <  2:1373631/151‑1 (MQ=255)
         ctTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTc                                                                                                                                              >  2:967449/1‑151 (MQ=255)
         ctTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGAAACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGt                                                                                                                                                <  2:1289346/149‑1 (MQ=255)
           tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgc                                                                                                                                            >  2:1086912/1‑151 (MQ=255)
             ccAAACCGGGATGCGGCTTTTCAATGAAATGCCTAGATTGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCTCATCGGTTTATCTAAGCAGAGTTCAGCCAGCCGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgccc                                                                                                                                          <  4:332997/151‑1 (MQ=255)
                 aCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCa                                                                                                                                      <  2:642651/151‑1 (MQ=255)
                   cGGGATGCGGCTTTTCATAAAAATGCCTAGACCGCCATGAATAACCCCCAAAATTAAAGGGAAATCCTGCTGCCACATCGGTTTATCTAGGCCGATTTCAGCCAGCAGCTGGGCGTGACGCTCGGGGGAACCCCCCCGTTCGCCCGCCATg                                                                                                                                    <  4:843049/151‑1 (MQ=255)
                      gATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatc                                                                                                                                 >  2:817126/1‑151 (MQ=255)
                       aTGCGGCTTTTCATTAAAATGCCTAGCCCGCCATGCAAACCGCCCAAAATGTAAGGGAGATCCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCGGGGGAGATCCCACGTCCGCCCGCCATGatc                                                                                                                                 <  3:515756/150‑1 (MQ=255)
                        ggCGGTTTTTCATTAAATTGCCTAGTCCGACATGCATAACCCCCAAATTGTAAGGGAGATCCGGCTGCCCAATCGTTTTATATAACCCGAGTCCGGCCGGCGGCGGCGGGGGACGCTCTGGGGAGACCCCCCGTTCGCCCGCCATGatca                                                                                                                                <  3:324648/149‑1 (MQ=255)
                          cGGTTTTTCATTGAAATGCCTAGATCGCCATGCATACCCCCCAAATTGTAAGGGAGATCCGCCTGCCACATCGTTTTATCTAACCCGAGTCCGGCCAGCAGCGGCGCGGGCCGCTGGGGGAGACCCCCCCGTCCCCCCGCCATGAtcatca                                                                                                                             <  4:1279710/151‑1 (MQ=255)
                              ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCAGGTTCGCCCGCCATGATCATCACCgg                                                                                                                         >  2:1460821/1‑151 (MQ=255)
                               tttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                         >  4:862742/1‑150 (MQ=255)
                                ttCATTAAAATCCCTAGACCGCCATGAAAAACCCCCAAAATGTAAGGGAGATCCTGCTCCCACATGGTTTTATCAAAGCCGATTTCAGCCAGCAGCTGCGCGTGCCGCTGGGGGGAGATCCCCCGTTCGCCCGCCATGATCATCACCGGAt                                                                                                                       <  1:16223/151‑1 (MQ=255)
                                       aaaTCCCTAGACCGCCATCAAAACCGCCCAAAATGTAATGGAGATCCTGCTGCCCCTTGGTTTTATCAAAGCCGAGTTCAGCCGGCAGCTGCGCGTGCCGCTCGGGGGATACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCgg                                                                                                                <  4:62461/151‑1 (MQ=255)
                                       aaaTCCCTAAATCGCCATGCAAACCCCCCCAAATTAAATGAAGATCCTGCTCCCACATCGTTTTACCAAGCCCATTTTCAGCCGGCAGCTGCGCGTGCCGCTCTGGGGAGACCCCACTTCCCCCCGCCATGATCATCACCGGATCGCCCgg                                                                                                                <  2:635467/151‑1 (MQ=255)
                                          tGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGAt                                                                                                             >  1:425440/1‑151 (MQ=255)
                                            ccTAGACCCCCATCCAAAACCCCCCAAATAAAATGGAGATCCTGCTCCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCGGCTGCGCGGGCCGCTCGGGGGAGACCCCACTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCa                                                                                                           <  4:1350152/151‑1 (MQ=255)
                                             cTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                          >  4:55534/1‑151 (MQ=255)
                                               aGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                        >  4:1404514/1‑151 (MQ=255)
                                               aGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                        >  4:642981/1‑151 (MQ=255)
                                                      cATGCATACCCCCCAAAATGTAAGGGAGATACGGCTCCCACATCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTGTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaaa                                                                                                 <  1:101103/151‑1 (MQ=255)
                                                       atccATACCCCCCAAAATGTAAGGGAGATCCTGCTCCCACTTGGTTTTATCAAAGCCGATTCCAGCCCGCGGGTGCGGGGGCCGCTCGGGGGAGACCCCACGGTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAg                                                                                                <  4:391977/148‑1 (MQ=255)
                                                        tGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                                               >  2:691092/1‑151 (MQ=255)
                                                         gCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                                               >  3:106114/1‑150 (MQ=255)
                                                         aaaaaccccccaAAATGTAATGGAGATCCGGCGGCCACAGCGTTTTATCAAAGCCGAGTCCAGCCAGCAGCGGCGGGTGCCGCTCGGGGGAGACCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGc                                                                                              <  3:207436/143‑1 (MQ=255)
                                                          cATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGc                                                                                              >  4:565032/1‑150 (MQ=255)
                                                             aaCGCCCAAAATGTAATGGAGATACTGCTCCCACATGGTTTTATCTAAGCCGATTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAt                                                                                          <  2:906239/151‑1 (MQ=255)
                                                              ccccccaAAATGAAGTGGAGATCCGGCTGCCACATGGTTTTATCAACCCCGGGTTCAGCCGGCAGCGGCGGGTGCCGGTCGGGGGAGACCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                         <  3:522045/148‑1 (MQ=255)
                                                              ccGCCCAAAATGTAATGGAGATACTGCGGCCACGTGGTTTTCTCTACGCGGAGTGCAGCAAGCGGCGGCGCGCGCCGCTGGGGGGAGATCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                         <  3:19437/150‑1 (MQ=255)
                                                              accccccAAATGTAATGGAGATACTGCTCCCACATCGTTTTATCAAACCCGAGTTCACCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                         <  1:378170/148‑1 (MQ=255)
                                                                cccccAATGGTAGTGGAGTTCCTGCTCCCACTCCGTTTTTTTTACGCGGCGTCCGCCCGGAAGCTGCGGGTACCGCTGGGGGGAAACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                       <  4:1372362/150‑1 (MQ=255)
                                                                cccccAAATGAAGTGGAAATCCGGCTCCCGCATTGTTTTATCAAAGCGATTTACAGCCGCCGGGTGCGGGTGCCGCTGGGGGGATCCCCCCCTTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                       <  4:1313202/150‑1 (MQ=255)
                                                                 ccccAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                                      >  3:157772/1‑151 (MQ=255)
                                                                 cccaAAATGAAAGGGATTTTCGGCCGCCACAGCGTTTTATTAACGCCGTGTTCACCCGGCAGCTGCGCGGGGCCCCCGGGGGAGCCCCCGGGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                                      <  2:1280342/151‑1 (MQ=255)
                                                                   ccAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                       >  2:108952/1‑148 (MQ=255)
                                                                   caaattgaaagggAGTTTCTGCTGCCACTTCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                                    <  3:303866/140‑1 (MQ=255)
                                                                   caaaaTGAAAGGGAGACACTGCGGCCAAATGGTTTTACCGAACCCGGTTCCACCCACTGACTGCGCGTAACGTATTGGGGAACCCCCAAGTTCACCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                                    <  2:395885/149‑1 (MQ=255)
                                                                    cAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                                    >  4:1222955/1‑150 (MQ=255)
                                                                     aaaTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                                                  >  3:89619/1‑151 (MQ=255)
                                                                       aTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                >  3:565608/1‑151 (MQ=255)
                                                                       aTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                >  3:531855/1‑151 (MQ=255)
                                                                         gaaagggAGATACTGCTCCCACATCGTTTTATCAAAGCCGAGTCCAGCCAGCAGCTGCGCGGGCCGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                              <  2:883819/146‑1 (MQ=255)
                                                                         gTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                              >  3:40059/1‑151 (MQ=255)
                                                                         gTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                              >  3:139000/1‑151 (MQ=255)
                                                                         gTAAGGGAGATCCTGCTGCCACATCGTTTTATCTAAGCGGAGTCCAGCCAGCGGCGGCGCGTGACGCTCGGGGGAGACCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                              <  2:646298/151‑1 (MQ=255)
                                                                         gTAAGGGAGATCCTCCTGCCACTTCGGTTTATAAACGCCGATTCCACCCGCCACCTGCGCGGGCCGCTCTGGGGACACCCCGCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                               <  4:900598/150‑1 (MQ=255)
                                                                           aaTGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAAt                                                                            >  4:1108235/1‑151 (MQ=255)
                                                                            aTGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAAt                                                                            >  4:928278/1‑150 (MQ=255)
                                                                             gggAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCTGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  2:1041617/150‑1 (MQ=255)
                                                                              ggggagacgccccccACAGCGTTTTATCAAAGCGGAGTTCCACCACCAGGTGCGCGTCCCCCTCGGGGGAGACCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  3:506967/137‑1 (MQ=255)
                                                                              ggAGATCCTGCTCCCACATCGTTTTATCTAACCCGATTCCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  1:11739/150‑1 (MQ=255)
                                                                              ggAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGACGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  1:349496/150‑1 (MQ=255)
                                                                              ggAGATACTGCTGCCACATCGTTTTATCTAAGCCAAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  1:524169/150‑1 (MQ=255)
                                                                              ggAGATACTGCTCCCACATCGTTTTATCTAAGCAGATTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCACCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          <  1:190778/150‑1 (MQ=255)
                                                                              gagttttctgctcCCGCAGCGTTTTATCGAAGCCGAGTTCGCCCAGCGGCTGGGCGTGCCGCTCTGGGGAGACCCCGCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAcc                                                                         <  2:689602/144‑1 (MQ=255)
                                                                               gagaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                        >  2:937990/1‑151 (MQ=255)
                                                                                agaTCCTGCCGCCACATCGTTTTATCAACGCCGATTCCAGCCAGCAGCGGCGGGGGCCGCGCGGGGGAGACCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                        <  3:131474/150‑1 (MQ=255)
                                                                                agaTCCTCCTGCCACACCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTGTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                         <  1:230208/133‑1 (MQ=255)
                                                                                agaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                        >  4:657652/1‑150 (MQ=255)
                                                                                  atcccgctccCCCATGGTTTTATCTAAGCCGCTTCCCCCAACCGTCTGCGCGTGACGCTTGGGGGCGCCCCCCCGCTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                                     <  2:944618/146‑1 (MQ=255)
                                                                                   tACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                    >  1:212430/1‑151 (MQ=255)
                                                                                   tACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                                     >  4:891655/1‑150 (MQ=255)
                                                                                    ccggcgcccACATCGTTTTATCTAACCCGAGTTCGGCCAGCAGCTGCGGGGGCCGCTCGGGGGAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                    <  2:1318961/143‑1 (MQ=255)
                                                                                    aCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                    >  4:732943/1‑150 (MQ=255)
                                                                                    aCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                    >  4:540238/1‑150 (MQ=255)
                                                                                    aCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                    >  4:301515/1‑150 (MQ=255)
                                                                                     ctgctgCCACATCGTTTTATCTAAGCCGAGTTCAGCCACCAGCTGCGCGGGCCGCTGGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAAc                                                                  <  3:453012/151‑1 (MQ=255)
                                                                                          gccgcgtCGTTTTATCTAAGCCGAGTTCAGCCGCCGGCTCCGCGGACCGCTGGGGGAGGCCCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTggg                                                             <  2:485680/145‑1 (MQ=255)
                                                                                           ccACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                            >  4:710507/1‑151 (MQ=255)
                                                                                            cccttCGTTTTACCAACGCCGATTTCAGCCGGCAGCGGCGCGTGACCCTCGGGGGAGCCCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGa                                                           <  3:535692/147‑1 (MQ=255)
                                                                                              caTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                                         <  2:1145330/151‑1 (MQ=255)
                                                                                                 ggttttATCAAAGCGGAGTCCACCCAGCAGGTGCGCGGGCCGCGGGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      <  4:1386637/148‑1 (MQ=255)
                                                                                                 cGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  2:1167484/1‑151 (MQ=255)
                                                                                                 cGGTTTACCTAAGCAGAGTTCAGCCAGCCCCGGCGCGTGCCGCGCTGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      <  2:1427673/151‑1 (MQ=255)
                                                                                                  ggTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  2:248077/1‑150 (MQ=255)
                                                                                                  ggTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGa                                                     >  1:27463/1‑151 (MQ=255)
                                                                                                    tttATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGt                                                   >  2:536361/1‑151 (MQ=255)
                                                                                                     ttATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTc                                                  >  4:721806/1‑151 (MQ=255)
                                                                                                      tATCTAGCCCGAGTCCAGCCGGCAGGGGCGGGTGCCGCTGGGGGGAGACCCCAGGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTc                                                  <  4:670624/150‑1 (MQ=255)
                                                                                                          tAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCCTCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGCTGACGAGTCCCaa                                              >  4:933770/1‑150 (MQ=255)
                                                                                                           aaGCCGAGTTCAGCCAGCAGCTGCGGGGGCCGCTCTGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                            <  2:831897/151‑1 (MQ=255)
                                                                                                           aaGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                            >  2:721676/1‑151 (MQ=255)
                                                                                                             gCCGATTTCAGCCAGCAGCTGCGCGTGACGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                            <  1:338596/149‑1 (MQ=255)
                                                                                                             gCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGt                                          >  3:552463/1‑151 (MQ=255)
                                                                                                               ggAGTCCACCGAGCGGCTGCGCGTGCCGCTTGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCg                                        <  2:227555/150‑1 (MQ=255)
                                                                                                                         gccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaa                                >  3:512696/1‑151 (MQ=255)
                                                                                                                          ccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaa                                >  2:371630/1‑150 (MQ=255)
                                                                                                                             gcggcTGCGCGTGACCCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACt                             <  3:506769/147‑1 (MQ=255)
                                                                                                                              cagcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGc                           >  2:378101/1‑151 (MQ=255)
                                                                                                                                  cTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc                        <  1:513441/151‑1 (MQ=255)
                                                                                                                                      gcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATaa                    >  1:185607/1‑151 (MQ=255)
                                                                                                                                        gTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAcc                  >  4:292517/1‑151 (MQ=255)
                                                                                                                                         ggccGCTGTGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGCACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAcc                  <  1:399640/147‑1 (MQ=255)
                                                                                                                                            cGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCgg               >  2:170207/1‑150 (MQ=255)
                                                                                                                                                ctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtct          >  4:619124/1‑151 (MQ=255)
                                                                                                                                                ctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtct          >  4:903103/1‑151 (MQ=255)
                                                                                                                                                 tGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtctc         >  2:1220463/1‑151 (MQ=255)
                                                                                                                                                         ctCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGt   <  3:286663/149‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                     
GCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705511‑3705810

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