Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,243,627 Δ1 bp coding (448/612 nt) emtA → lytic murein transglycosylase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,243,6240G.98.4% 254.8 / ‑2.4 62coding (445/612 nt)emtAlytic murein transglycosylase E
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (30/31);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGACCAGCGAGCTGAAAAATCCAGAGCGTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATA  >  NC_000913/1243487‑1243756
                                                                                                                                         |                                                                                                                                    
cgacCAGCGAGCTGAAAAATCCAGAGCGTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACggg                                                                                                                                    <  2:1019812/140‑1 (MQ=255)
                 aaTCCAGAGCGTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACgg                                                                                                                  <  1:645739/140‑1 (MQ=255)
                      agagCGTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACttt                                                                                                             <  2:495977/140‑1 (MQ=255)
                        agCGTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTtct                                                                                                           >  1:639352/1‑140 (MQ=255)
                           gTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGt                                                                                                        <  1:604854/140‑1 (MQ=255)
                           gTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGt                                                                                                        <  1:883638/140‑1 (MQ=255)
                                                   aCCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCa                                                                                >  1:406334/1‑140 (MQ=255)
                                                        aaTATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATc                                                                           <  2:957975/140‑1 (MQ=255)
                                                         atatTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCa                                                                          <  2:281287/140‑1 (MQ=255)
                                                          tatTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCaa                                                                         >  1:1102858/1‑140 (MQ=255)
                                                             tCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGa                                                                      >  1:1399280/1‑140 (MQ=255)
                                                             tCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGa                                                                      >  1:452547/1‑140 (MQ=255)
                                                              cTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGa                                                                      >  2:177041/1‑140 (MQ=255)
                                                               tGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGAtt                                                                    >  2:1222187/1‑140 (MQ=255)
                                                                  aaaCCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTgg                                                                 >  2:1217880/1‑140 (MQ=255)
                                                                  aaaCCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTgg                                                                 >  2:586940/1‑140 (MQ=255)
                                                                   aaCCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGa                                                                <  1:586940/140‑1 (MQ=255)
                                                                    aCCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGc                                                                                 >  2:147761/1‑122 (MQ=255)
                                                                    aCCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGc                                                                                 <  1:147761/122‑1 (MQ=255)
                                                                    aCCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGAt                                                               >  2:333327/1‑140 (MQ=255)
                                                                        gCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTg                                                           >  2:293983/1‑140 (MQ=255)
                                                                          ccGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGAc                                                         <  1:1160613/140‑1 (MQ=255)
                                                                             cTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAg                                                      <  2:1102858/140‑1 (MQ=255)
                                                                                gcaggcaTTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTc                                                   <  2:639352/140‑1 (MQ=255)
                                                                                        tGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAaca                                           >  2:709826/1‑140 (MQ=255)
                                                                                         gAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAacac                                          <  2:1399280/140‑1 (MQ=255)
                                                                                           aGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCt                                                 <  2:132558/131‑1 (MQ=255)
                                                                                           aGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCt                                                 >  1:132558/1‑131 (MQ=255)
                                                                                           aGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGt                                        <  1:286858/140‑1 (MQ=255)
                                                                                               ccGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAgcg                                    >  1:163392/1‑140 (MQ=255)
                                                                                               ccGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAgcg                                    >  1:1405664/1‑140 (MQ=255)
                                                                                               ccGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAgcg                                    <  1:1217880/140‑1 (MQ=255)
                                                                                                       aCTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCa                            >  1:451581/1‑140 (MQ=255)
                                                                                                       aCTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCa                            >  1:1378569/1‑140 (MQ=255)
                                                                                                         tGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTTACGAGTTCCTCGAACACGTAGCGCGAAATCAcc                          >  2:1144136/1‑140 (MQ=255)
                                                                                                         tGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCAcc                          >  1:1144285/1‑140 (MQ=255)
                                                                                                         tGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCAcc                          >  2:62288/1‑140 (MQ=255)
                                                                                                         tGCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCAca                          >  2:999496/1‑139 (MQ=255)
                                                                                                          gCAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCAccc                         >  1:1012390/1‑140 (MQ=255)
                                                                                                           cAATATGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCt                        <  2:1386469/140‑1 (MQ=255)
                                                                                                             atatGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTgc                      <  2:1012390/140‑1 (MQ=255)
                                                                                                             atatGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTgc                      <  2:1144285/140‑1 (MQ=255)
                                                                                                             atatGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTgc                      <  1:37330/140‑1 (MQ=255)
                                                                                                                tGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGcgcc                   >  1:1079392/1‑140 (MQ=255)
                                                                                                                tGCGCTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGcgcc                   >  2:919238/1‑140 (MQ=255)
                                                                                                                  cgcTGGTGGTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGcgccgc                 >  1:195248/1‑140 (MQ=255)
                                                                                                                     tggtggTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAgg              >  2:1100704/1‑140 (MQ=255)
                                                                                                                     tggtggTGTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAgg              <  1:1222187/140‑1 (MQ=255)
                                                                                                                       gtggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCGGACCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCt            <  1:1254287/140‑1 (MQ=255)
                                                                                                                       gtggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCt            <  2:452547/140‑1 (MQ=255)
                                                                                                                       gtggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCt            >  1:141855/1‑140 (MQ=255)
                                                                                                                       gtggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACACTGCGCCGCAGGCt            <  1:999496/140‑1 (MQ=255)
                                                                                                                        tggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTc           <  2:406334/140‑1 (MQ=255)
                                                                                                                        tggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTc           >  1:713860/1‑140 (MQ=255)
                                                                                                                        tggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTc           <  2:163392/140‑1 (MQ=255)
                                                                                                                        tggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTc           <  2:713860/140‑1 (MQ=255)
                                                                                                                        tggtgTCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTc           <  1:1100704/140‑1 (MQ=255)
                                                                                                                             tCATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCgcgc      <  1:293983/140‑1 (MQ=255)
                                                                                                                              cATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCt     <  2:195248/140‑1 (MQ=255)
                                                                                                                              cATACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCt     <  2:1079392/140‑1 (MQ=255)
                                                                                                                                tACGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCtat   <  2:451581/140‑1 (MQ=255)
                                                                                                                                 aCGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCcgcg       >  1:557745/1‑135 (MQ=255)
                                                                                                                                 aCGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCcgcg       <  2:557745/135‑1 (MQ=255)
                                                                                                                                 aCGCTAACGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCtata  >  1:111539/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CGACCAGCGAGCTGAAAAATCCAGAGCGTAATATTTCAATGGGGGCGGCTTACCTGAATATTCTGGAAACCGGCCCGCTGGCAGGCATTGAAGATCCGAAGGTACTGCAATATGCGCTGGTGGTGTCATACGCTAACGGGGCAGGTGCGCTGCTACGGACTTTCTCGTCAGATCGGAAAAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATA  >  NC_000913/1243487‑1243756

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