Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6610TC44.3% 123.3 / 23.1 79I44V (ATC→GTC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (41/3);  new base C (0/35);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 2.59e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705513‑3705810
                                                                                                                                                    |                                                                                                                                                     
aCGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATcc                                                                                                                                                     <  1:18997/151‑1 (MQ=255)
aCGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATcc                                                                                                                                                     >  2:550863/1‑151 (MQ=255)
  gAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATccc                                                                                                                                                    >  1:482694/1‑150 (MQ=255)
   aaCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCAc                                                                                                                                                  <  1:127846/151‑1 (MQ=255)
    aCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACg                                                                                                                                                 >  2:938598/1‑151 (MQ=255)
                ccGGGATGCGTTTTTTCATTAAAATCCCTAGACCGCCATACATACCGCCCAAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAACCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTCTGGGGAGACCCCACGTTCCCCCGCCAt                                                                                                                                     <  2:35246/151‑1 (MQ=255)
                   ggATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGat                                                                                                                                  >  2:306421/1‑151 (MQ=255)
                          gCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCAcc                                                                                                                           >  1:355267/1‑151 (MQ=255)
                             tttAATTAAAAGGCCAAGTCCGCCATGCATACCCCCCAAAATGAAGGGGAGATCCTGCTGCCACATCGTTTTACCTACGCCGAGTACGGCCGGCATCTGCGCGGGCCGCTCGGGGGAGTCCCCCCGTTCGCCCGCCATGATCATCACCg                                                                                                                          <  2:212319/149‑1 (MQ=255)
                                 aTTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGc                                                                                                                    >  2:619223/1‑151 (MQ=255)
                                     aaaTGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCgg                                                                                                                >  2:157970/1‑151 (MQ=255)
                                     aaaTCCCTAGATCCCCATGCATAACCCCCAAAATGTAAGGGAGATCCTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGCCCCCCCGTTCGCCCGCCATGATCATCACCGGATCGccc                                                                                                                  <  1:737131/149‑1 (MQ=255)
                                     aaaTCCCTAGACCGCATTCCATACCCCCCAAAATGAAAGGGAGATACTGCCGCCCCCTCGTTTAATAAAAGCCGGGTCCAGCCAGCAGCTGCGCGGGCCGCTGGGGGAAGGCCCCACGGTCGCCCGCCATGATCATCACCGGa                                                                                                                        <  2:753308/143‑1 (MQ=255)
                                      aaTGCCTGGATCGCCTTGCAAACCGCCCAAAATGAAAGGGAGATCCTGCGGCCACATGGTTTTATCTAAGCCGAGTTCGGCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                               <  2:167313/151‑1 (MQ=255)
                                       aTGCCTAGACCCCCATGCATACCGCCCAAAATGTAATGGAGATACTGCTGCCACATGGTTTTATCTAAGCCGATTTCAGCCAGCAGCGGCGCGTGACGCTCTGGGGAGCCCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGa                                                                                                              <  1:762634/151‑1 (MQ=255)
                                        tGCCTAGATCGCCATGCATAACCCCCCAAAGGTAGGGGAGATCCTGCGGCCCCTGCGTTTTATCAAGGCCGAGTTCAGCCAGCGGCGGCGCGTGCCGCTCGGGGGAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGa                                                                                                              <  2:224741/150‑1 (MQ=255)
                                         gCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                            >  2:58785/1‑151 (MQ=255)
                                          ccTAGCTCCCCATGCATAACCCCCAAAATGTAAGGGAGATACTGCTGCCACATCGGTTTATCAAAGCCGATTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                            <  1:298046/150‑1 (MQ=255)
                                           cTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                          >  1:959213/1‑151 (MQ=255)
                                            tAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtg                                                                                                         >  1:723733/1‑151 (MQ=255)
                                             aGACCGCCATGCATAACCCCCCAAATTAAATGGGGATCCTGCTGCCACATCGTTTTATCTAGGCCGAGTTCACCCAGCAGCTGCGGGGGACGCTTGGGGGAGTCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtg                                                                                                         <  1:996998/150‑1 (MQ=255)
                                                   ccATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACa                                                                                                   >  2:391863/1‑150 (MQ=255)
                                                   ccATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACa                                                                                                   >  1:917858/1‑150 (MQ=255)
                                                   ccATGCAAACCCCCCAAAAGGAAAGGGAGATCCTGCTGCCCCATGGTTTTATCTAGGCCGAGTCCAGCCAGCGGCTGCGCGGGCCGCTCGGGGGAGATCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtg                                                                                                         <  2:724780/144‑1 (MQ=255)
                                                     aTGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa                                                                                                  >  1:482258/1‑149 (MQ=255)
                                                     aTGCATAACCCCCAAAATGTAAGGGAGATCCTGCTGCCACATCGTTTTATCTACCCGGAGTTCACCCGGCGGGTGCGCGGGCCGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAg                                                                                                <  2:202277/151‑1 (MQ=255)
                                                           accccccaAAATGTAATGGAGATCCTGCCCCCCCTTCGTTTTATCTAACCCGATTCCAGCCCGCAGCGGGGCGCGCCCCCGGGGGGAAACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGc                                                                                              <  2:962868/143‑1 (MQ=255)
                                                           aaCGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAt                                                                                          >  2:975845/1‑151 (MQ=255)
                                                            aCGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                         >  2:166546/1‑151 (MQ=255)
                                                             ccccccAAATGTAAGGGAGATCCTGCTCCCACTTCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                         <  1:924107/148‑1 (MQ=255)
                                                               cccaAAATGTAAGGGAGATACTGCCCCCACCTGGTTTTATCAAAGCCGAGTCCAGCCAGCAGCGGCGCGGGACGCTGGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                                      <  2:490230/151‑1 (MQ=255)
                                                                cccAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                       >  2:322707/1‑149 (MQ=255)
                                                                    aaTGTAAGGGAGATCCTGCTCCCACATCGGTTTATCTAAGCGTAGTTCAGCCAGCAGCTGCGGGTGCCCCTCTGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                                 <  1:714846/151‑1 (MQ=255)
                                                                     aTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                >  2:6418/1‑151 (MQ=255)
                                                                     aTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                >  2:982472/1‑151 (MQ=255)
                                                                       gTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                               >  2:959586/1‑150 (MQ=255)
                                                                        tAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                              >  2:846541/1‑150 (MQ=255)
                                                                         aagggAGATTCTGCGCCCACATGGTTTTATAAAACCCGATTTCAGCCAGCAGCTGCGGGTGACGCTCGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTaa                                                                             <  2:744685/147‑1 (MQ=255)
                                                                         aagggAGATACTGCTGCCACATCGTTTTATCTAAGCCGATTTCAGCCAGCAGCTGCGCGGGCCGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                <  1:978197/144‑1 (MQ=255)
                                                                         aaggaagaTCCTGCCCCCACATGGTTTTATCGAACCCGAGTCCGCCCGGCAGCTGCGGGGGCCGCTTGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                               <  2:691563/143‑1 (MQ=255)
                                                                         aaTGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAAt                                                                            >  1:271865/1‑151 (MQ=255)
                                                                           tGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                          >  2:439553/1‑151 (MQ=255)
                                                                            ggAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAcc                                                                         >  1:564578/1‑151 (MQ=255)
                                                                             gagaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                        >  2:535458/1‑151 (MQ=255)
                                                                                aTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                      >  2:668908/1‑150 (MQ=255)
                                                                                  aCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                                    >  2:77773/1‑150 (MQ=255)
                                                                                   ctgctgCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAAc                                                                  >  1:908687/1‑151 (MQ=255)
                                                                                         ccACATCGTTTTATCAAAGCCGAGTTCACCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                               <  1:585266/148‑1 (MQ=255)
                                                                                          cacaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                            >  1:851321/1‑150 (MQ=255)
                                                                                           acaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                          >  1:434525/1‑151 (MQ=255)
                                                                                            caTCGTTTTATCAAAGCCGAGTTCACCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                          <  1:298203/150‑1 (MQ=255)
                                                                                                gttttATCAAAGCCGATTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                          <  1:28578/144‑1 (MQ=255)
                                                                                                ggTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                                      >  2:618721/1‑150 (MQ=255)
                                                                                                  tttATCAAAGCCGAGTGCAGCCAGCAGCTGCGCGTGCCGCTCTGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                          <  2:196738/144‑1 (MQ=255)
                                                                                                     aTCTAAGCCGATTCCAGCCAGCGGCTGCGGGTGACGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc                                                 <  1:811641/150‑1 (MQ=255)
                                                                                                        tAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCaa                                              >  2:327959/1‑150 (MQ=255)
                                                                                                          aGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACg                                           >  2:27574/1‑151 (MQ=255)
                                                                                                            ccGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                                         >  1:887234/1‑151 (MQ=255)
                                                                                                            ccGAGTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGCCCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                            <  1:797063/148‑1 (MQ=255)
                                                                                                             cGAGTTCAGCCAGCGGCTGCGCGTGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                            <  1:798854/147‑1 (MQ=255)
                                                                                                               tgttccccccgcagcGGCGGGGGCCGCTGGGGGGAGACCCCAGGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCaaa                                             <  2:722560/134‑1 (MQ=255)
                                                                                                               aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                       >  2:146864/1‑150 (MQ=255)
                                                                                                                 ttCAGCCAGCAGCTGCGCGTGCCGCTCTGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                      <  1:825844/149‑1 (MQ=255)
                                                                                                                 ttCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCg                                        <  1:2596/147‑1 (MQ=255)
                                                                                                                  tCACCCAGCAGCGGCGCGTGCCGCGCTGGGGAAACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgaga                                    <  1:51751/150‑1 (MQ=255)
                                                                                                                     gccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATa                                 >  2:52350/1‑150 (MQ=255)
                                                                                                                          cagcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGc                           >  1:247646/1‑151 (MQ=255)
                                                                                                                               gcggGTGACGCTGTGGGAAGACCCCTCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACAt                      <  1:567964/151‑1 (MQ=255)
                                                                                                                               gcgcGGGCCGCTGGGGGGAACCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                                         <  2:865914/132‑1 (MQ=255)
                                                                                                                                 gcgtgccgcgggggggATACCCCCTGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACt                             <  2:724608/130‑1 (MQ=255)
                                                                                                                                 gcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATa                     >  2:601478/1‑150 (MQ=255)
                                                                                                                                 gcgTGACGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGc                           <  1:576096/144‑1 (MQ=255)
                                                                                                                                   gggcccctctGGGGAACCCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATaa                    <  1:619289/143‑1 (MQ=255)
                                                                                                                                         ctctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAtt            >  1:639586/1‑151 (MQ=255)
                                                                                                                                            tGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtctc         >  1:30571/1‑151 (MQ=255)
                                                                                                                                             ggggAGCCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct        <  2:959088/151‑1 (MQ=255)
                                                                                                                                             ggggAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtctc         <  2:905720/150‑1 (MQ=255)
                                                                                                                                                gagaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtat      >  1:284796/1‑150 (MQ=255)
                                                                                                                                                   aTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTg  >  1:200936/1‑151 (MQ=255)
                                                                                                                                                    tCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTg  >  1:612038/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTG  >  NC_000913/3705513‑3705810

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