Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 725,836 Δ1 bp 5.7% coding (1201/2049 nt) kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913725,8360T. 5.7% 232.0 / 6.3 70coding (1201/2049 nt)kdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base . (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGATATCTTTCAGCGCAATAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGTTAAAACGCTGCTTGGCGAGGAT  >  NC_000913/725690‑725977
                                                                                                                                                  |                                                                                                                                             
gACGATATCTTTCAGCGCAATAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATg                                                                                                                                              <  3:156076/148‑1 (MQ=255)
 aCGATATCTTTCAGCGCAATAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGc                                                                                                                                            <  4:293170/149‑1 (MQ=255)
       tCTTTCAGCGCAATAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGa                                                                                                                                       <  7:262128/148‑1 (MQ=255)
                   aTAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACg                                                                                                                          <  2:369894/149‑1 (MQ=255)
                    tAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTa                                                                                                                            >  7:342284/1‑146 (MQ=255)
                    tAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTa                                                                                                                            <  8:342284/146‑1 (MQ=255)
                          ccAGCACACGAGAACCTTCCACCACCAGCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATg                                                                                                                   >  2:66908/1‑149 (MQ=255)
                               acacGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGtt                                                                                                              >  7:250348/1‑149 (MQ=255)
                                 acGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACg                                                                                                                          >  7:300509/1‑135 (MQ=255)
                                 acGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACg                                                                                                                          <  8:300509/135‑1 (MQ=255)
                                 acGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGt                                                                                                            >  8:251318/1‑149 (MQ=255)
                                   gagaACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCgg                                                                                                                <  3:209229/143‑1 (MQ=255)
                                   gagaACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCgg                                                                                                                >  4:209229/1‑143 (MQ=255)
                                   gagaACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                                                                                                           >  1:172362/1‑148 (MQ=255)
                                   gagaACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTGAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc                                                                                                           >  7:188515/1‑148 (MQ=255)
                                     gacccTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGAt                                                                                                        <  3:327161/146‑1 (MQ=255)
                                          ttCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCttt                                                                                                                             >  2:273755/1‑123 (MQ=255)
                                          ttCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCttt                                                                                                                             <  1:273755/123‑1 (MQ=255)
                                            ccaccaccaccaGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGa                                                                                                   <  6:3271/147‑1 (MQ=255)
                                                caccaccaGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGc                                                                                             <  6:351434/149‑1 (MQ=255)
                                                   caccaGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGc                                                                                             >  5:215317/1‑146 (MQ=255)
                                                   caccaGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGc                                                                                             <  6:215317/146‑1 (MQ=255)
                                                            cgTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCttt                                                                                 <  8:188515/149‑1 (MQ=255)
                                                             gTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGAt                                                                                                        >  7:134774/1‑125 (MQ=255)
                                                             gTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGAt                                                                                                        <  8:134774/125‑1 (MQ=255)
                                                                  tCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAg                                                                           >  2:370054/1‑149 (MQ=255)
                                                                    ccTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGt                                                                          >  5:92285/1‑148 (MQ=255)
                                                                    ccTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGt                                                                          >  3:255971/1‑148 (MQ=255)
                                                                    ccTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGt                                                                          >  4:131311/1‑148 (MQ=255)
                                                                       gACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAAc                                                                      >  3:262548/1‑149 (MQ=255)
                                                                        aCGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACg                                                                     <  4:255971/149‑1 (MQ=255)
                                                                         cgcgCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATg                                                                                                                                              >  2:336422/1‑75 (MQ=255)
                                                                         cgcgCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCTGGCTTTGCGCAGTAAACgg                                                                    <  1:66908/149‑1 (MQ=255)
                                                                           cgcAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTa                                                                  <  1:370054/149‑1 (MQ=255)
                                                                           cgcAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTa                                                                  <  4:262548/149‑1 (MQ=255)
                                                                                   gATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACgg                                                                    <  8:199810/139‑1 (MQ=255)
                                                                                   gATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACgg                                                                    >  7:199810/1‑139 (MQ=255)
                                                                                   gATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTg                                                          >  7:382896/1‑149 (MQ=255)
                                                                                   gATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTg                                                          >  2:49046/1‑149 (MQ=255)
                                                                                    aTCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATc                                                                                                 >  7:388416/1‑109 (MQ=255)
                                                                                    aTCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTgg                                                         <  8:382896/149‑1 (MQ=255)
                                                                                    aTCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTgg                                                         >  3:137008/1‑149 (MQ=255)
                                                                                     tCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTgg                                                         >  1:190045/1‑148 (MQ=255)
                                                                                       gACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGACCCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCAt                                                      >  6:212629/1‑149 (MQ=255)
                                                                                        aCTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATg                                                     <  2:172362/149‑1 (MQ=255)
                                                                                            tttGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAg                                                  <  7:251318/148‑1 (MQ=255)
                                                                                             ttGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGc                                                 <  6:92285/148‑1 (MQ=255)
                                                                                                 tCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTg                                            >  1:76252/1‑149 (MQ=255)
                                                                                                   aaCATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGca                                          >  1:70902/1‑149 (MQ=255)
                                                                                                      aTCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTAGTCGATGTTGATCCCGCTCATCCGGCATAGCGCAGAAAACGGTACAAAGGTGGCATGGAGCGACTGcaca                                        >  8:335865/1‑148 (MQ=255)
                                                                                                               gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTg                                                          <  2:21669/121‑1 (MQ=255)
                                                                                                               gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTg                                                          >  1:21669/1‑121 (MQ=255)
                                                                                                               gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACAtcgcgctcgc                              <  7:75665/149‑1 (MQ=255)
                                                                                                               gAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACAtcgcgctcgc                              <  3:131311/149‑1 (MQ=255)
                                                                                                                 aGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACAtcgcgctcgcgc                            <  8:250348/149‑1 (MQ=255)
                                                                                                                  gTGACCACCGTTAGCCTCAACATGGCGACGAAGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCttt                                                                                 <  1:161825/94‑1 (MQ=255)
                                                                                                                  gTGACCACCGTTAGCCTCAACATGGCGACGAAGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCttt                                                                                 >  2:161825/1‑94 (MQ=255)
                                                                                                                  gTGACCACCGTTAGCCTCAACATGGCGACGAAGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCttt                                                                                 >  2:156529/1‑94 (MQ=255)
                                                                                                                  gTGACCACCGTTAGCCTCAACATGGCGACGAAGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCttt                                                                                 <  1:156529/94‑1 (MQ=255)
                                                                                                                      ccaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCa                           <  4:291979/145‑1 (MQ=255)
                                                                                                                      ccaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCa                           >  3:291979/1‑145 (MQ=255)
                                                                                                                       caccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGtt                       >  6:222206/1‑148 (MQ=255)
                                                                                                                       caccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGtt                       >  5:305668/1‑148 (MQ=255)
                                                                                                                       caccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGtt                       >  8:3681/1‑148 (MQ=255)
                                                                                                                                 cTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGTTAAAACgctgct            >  2:228840/1‑149 (MQ=255)
                                                                                                                                  tCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGAc                                              <  7:223738/114‑1 (MQ=255)
                                                                                                                                  tCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGAc                                              >  8:223738/1‑114 (MQ=255)
                                                                                                                                          gCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGTTAAAACGCTGCTTGGCGAGGa   <  2:70902/149‑1 (MQ=255)
                                                                                                                                          gCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGTTAAAACGCTGCTTGGCGAGGa   <  1:49046/149‑1 (MQ=255)
                                                                                                                                            gacgaATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCTGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGTTAAAACGCTGCTTGGCGAGGAt  >  6:120249/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GACGATATCTTTCAGCGCAATAACGCCCAGCACACGAGAACCTTCCACCACCACCAGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTCGATGTTGATCCCGCTCATCCGGCTTTGCGCAGTAAACGGTACAAAGGTGGCATGGAGCGACTGCACATCGCGCTCGCGCAGGTTAAAACGCTGCTTGGCGAGGAT  >  NC_000913/725690‑725977

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