Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,624,333 (T)5→6 coding (165/393 nt) ydeI ← hydrogen peroxide resistance OB fold protein, putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,624,3281.T98.6% 219.8 / ‑4.1 69coding (170/393 nt)ydeIhydrogen peroxide resistance OB fold protein, putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (33/35);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCCATCAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GG‑TTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAT  >  NC_000913/1624199‑1624458
                                                                                                                                   |                                                                                                                                  
cGTCCATCAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GG‑ttttttgcg                                                                                                                           <  2:1521724/139‑5 (MQ=255)
cGTCCATCAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GG‑ttttttgcg                                                                                                                           >  1:897230/1‑135 (MQ=255)
   ccATCAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGaaa                                                                                                                         >  2:1585483/1‑139 (MQ=255)
      tCAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCa                                                                                                                      >  1:981179/1‑139 (MQ=255)
       cAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCaa                                                                                                                     <  1:1653859/139‑1 (MQ=255)
               gCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAtt                                                                                                             >  1:895875/1‑139 (MQ=255)
                ccGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAAc                                                                                                                    >  2:629336/1‑131 (MQ=255)
                ccGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAAc                                                                                                                    <  1:629336/131‑1 (MQ=255)
                ccGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAttt                                                                                                            >  2:107259/1‑139 (MQ=255)
                      gAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                      <  2:513871/139‑3 (MQ=255)
                                 ttAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTc                                                                                          >  1:688548/1‑139 (MQ=255)
                                         aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                  <  2:981179/139‑1 (MQ=255)
                                         aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                  <  2:688548/139‑1 (MQ=255)
                                         aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                  <  2:897230/139‑1 (MQ=255)
                                           cGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                <  2:1619005/139‑1 (MQ=255)
                                              tctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                             <  2:976596/139‑1 (MQ=255)
                                              tctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                             <  2:895875/139‑1 (MQ=255)
                                              tctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                             <  2:273178/139‑1 (MQ=255)
                                               ctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                            >  2:2049978/1‑139 (MQ=255)
                                                  aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                >  1:142039/1‑132 (MQ=255)
                                                  aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                <  2:142039/132‑1 (MQ=255)
                                                          aaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCACTTTTTTTGCATCATCGGt                                                                                            >  1:310183/1‑113 (MQ=255)
                                                          aaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCACTTTTTTTGCATCATCGGt                                                                                            <  2:310183/113‑1 (MQ=255)
                                                              gTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggc                                                             <  1:1341416/139‑1 (MQ=255)
                                                                     tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAt                                                                              >  1:28222/1‑115 (MQ=255)
                                                                     tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAt                                                                              <  2:28222/115‑1 (MQ=255)
                                                                     tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                      <  2:771290/139‑1 (MQ=255)
                                                                     tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                      >  1:881647/1‑139 (MQ=255)
                                                                        tctcCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                      <  2:1329269/89‑3 (MQ=255)
                                                                        tctcCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                      >  1:1636264/1‑87 (MQ=255)
                                                                        tctcCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                      >  1:1329269/1‑87 (MQ=255)
                                                                        tctcCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                      <  2:1636264/89‑3 (MQ=255)
                                                                               tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                            <  2:1839468/139‑1 (MQ=255)
                                                                               tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                            <  2:313720/139‑1 (MQ=255)
                                                                               tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                            >  1:1839468/1‑139 (MQ=255)
                                                                                   agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACgg                                                                                                                  >  1:1679726/1‑66 (MQ=255)
                                                                                   agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACgg                                                                                                                  <  2:1679726/66‑1 (MQ=255)
                                                                                   agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCt                                                 <  1:1362950/130‑1 (MQ=255)
                                                                                   agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCt                                                 >  2:1362950/1‑130 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                <  1:1381236/80‑1 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                >  2:1381236/1‑80 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGAt                                                                                    >  2:649330/1‑92 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGAt                                                                                    <  1:649330/92‑1 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGtttt                                     >  1:1983472/1‑139 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGtttt                                     >  2:1920412/1‑139 (MQ=255)
                                                                                      aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGtttt                                     >  1:15430/1‑139 (MQ=255)
                                                                                       aTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTa                                    <  2:15430/139‑1 (MQ=255)
                                                                                              tACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGcc                             <  1:1920412/139‑1 (MQ=255)
                                                                                                ccGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCt                           >  1:328282/1‑139 (MQ=255)
                                                                                                   cgTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                  >  1:1815547/1‑81 (MQ=255)
                                                                                                   cgTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                  <  2:1815547/81‑1 (MQ=255)
                                                                                                   cgTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGt                        >  2:121922/1‑139 (MQ=255)
                                                                                                   cgTAGTGAAACCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGt                        >  1:1242378/1‑139 (MQ=255)
                                                                                                          aaaCCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGa                                                                <  1:1238443/92‑1 (MQ=255)
                                                                                                          aaaCCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGa                                                                >  2:1238443/1‑92 (MQ=255)
                                                                                                            aCCGAAGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCa               <  2:328282/139‑1 (MQ=255)
                                                                                                                aaGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAttttttt                                                                                                        >  1:1042410/1‑47 (MQ=38)
                                                                                                                aaGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCAttttttt                                                                                                        <  2:1042410/47‑1 (MQ=38)
                                                                                                                aaGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGc           >  1:381535/1‑139 (MQ=255)
                                                                                                                aaGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGc           >  1:644116/1‑139 (MQ=255)
                                                                                                                 aGCGCCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa          <  2:1242378/139‑1 (MQ=255)
                                                                                                                  gcgcCATCGTGCATAGT‑TTTTTGCGAAATCAACGGTCa                                                                                                               >  2:1250691/1‑38 (MQ=25)
                                                                                                                  gcgcCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCa                                                                                                               <  1:1250691/38‑1 (MQ=37)
                                                                                                                  gcgcCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                >  2:1164281/1‑52 (MQ=37)
                                                                                                                  gcgcCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                <  1:1164281/52‑1 (MQ=37)
                                                                                                                  gcgcCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCAt         >  2:697456/1‑139 (MQ=255)
                                                                                                                    gcCATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATaa       <  1:2049978/139‑1 (MQ=255)
                                                                                                                      cATCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAgg     >  2:1566429/1‑139 (MQ=255)
                                                                                                                       aTCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                            <  1:95101/99‑1 (MQ=255)
                                                                                                                       aTCGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                            >  2:95101/1‑99 (MQ=255)
                                                                                                                         cGTGCAT‑GGTTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt  <  2:881647/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                  
CGTCCATCAAAGACGGCCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCAT‑GG‑TTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAT  >  NC_000913/1624199‑1624458

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