Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,243,335 C→T 6.4% V52V (GTC→GTT)  emtA → lytic murein endotransglycosylase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,243,3350CT 6.4% 154.8 / 2.7 62V52V (GTC→GTT) emtAlytic murein endotransglycosylase E
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base T (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAGTGGTGAGC  >  NC_000913/1243188‑1243483
                                                                                                                                                   |                                                                                                                                                    
gagaTGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGATAGTTCCGGTGCAACGTGCGATGCACTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCg                                                                                                                                                     >  2:762091/1‑149 (MQ=255)
  gaTGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGACAGCCGGTGCAGCCTGGGGCGTCGAt                                                                                                                                                   <  2:375013/149‑1 (MQ=255)
  gaTGGTGTGCCTTTTTGATTGTGTTAATAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGAt                                                                                                                                                   <  1:785439/149‑1 (MQ=255)
          gCCTTTTTGATTGTGTGATTAGAGGGCTGTTCATCAACGCATGACTATACGAACCCGCCGTCGAACTAGAAAGTTCAGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCATAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAAtt                                                                                                                                           <  2:173851/149‑1 (MQ=255)
             tttttGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCCGCCTGGGGCGTCGATCCACAATTGAt                                                                                                                                        >  1:100085/1‑149 (MQ=255)
             tttttGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGAt                                                                                                                                        >  2:808441/1‑149 (MQ=255)
              ttttGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATc                                                                                                                                       >  2:263682/1‑149 (MQ=255)
                       tgtTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCACTAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGAtt                                                                                                                              >  1:101372/1‑149 (MQ=255)
                         ttattaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAATGTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTAt                                                                                                                            <  2:668560/149‑1 (MQ=255)
                         ttattaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTTCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTAt                                                                                                                            >  1:67013/1‑149 (MQ=255)
                         ttattaGAGTGTTGTTCATCAAAGCATGACTATACGAACCAGCCGTGGAACGCGAAAGTTACGGTGCAACGTGCGATGCAGTGTATGCCAATAAGCCAGAAAGCAGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTAt                                                                                                                            <  2:149272/149‑1 (MQ=255)
                           attaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCg                                                                                                                          >  1:318884/1‑149 (MQ=255)
                                     tgttCATCAAAGCATGACTATAAGAACCCGGCGTTGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATc                                                                                                                <  2:750397/149‑1 (MQ=255)
                                      gttCATCAAAGCATGACTATACGAACCCGCCGTCGAACGCGAAACTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCAAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAAATGATCACGGCGATTATCGCTATCTAATcc                                                                                                               >  2:752486/1‑148 (MQ=255)
                                       ttCATCAAAGCATGACTATACGAGCCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATACGCCAGAAAGCCGGTGCAGCCTGGGGGGTCGATCCTCATGTGCTGAGGGGGATTTTCGGTATCGGGTCgg                                                                                                              >  2:619766/1‑149 (MQ=255)
                                       ttCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATACGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCCCAATTGATCACGGCGATTATCGCTATCGAATCgg                                                                                                              >  2:867001/1‑149 (MQ=255)
                                         catccAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGGTCCGGTGCAACGGGCGAAGCAGTGGATGCCAATAAGCCACAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGggt                                                                                                            <  1:681395/149‑1 (MQ=255)
                                          atcaAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGCTCCACAATTGATCACGGCGATTATCGCTATCGATTCGggtg                                                                                                           >  1:509718/1‑149 (MQ=255)
                                            caAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGggtggt                                                                                                         >  2:847021/1‑149 (MQ=255)
                                               aGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAAt                                                                                                      <  1:976445/149‑1 (MQ=255)
                                                gCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATc                                                                                                     <  1:560518/149‑1 (MQ=255)
                                                 cATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATcc                                                                                                    <  2:1016364/149‑1 (MQ=255)
                                                      cTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAAcg                                                                                               >  1:658420/1‑149 (MQ=255)
                                                          aCGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGt                                                                                           <  1:702072/149‑1 (MQ=255)
                                                            gAACCCGCAGTGGAACTCTAAAGTTCCGGTTCAACGTTCGATGGAGTGGATGCCAATTAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCTCTATCGATTCGGGTGGTAATCCGAACGCGGTGa                                                                                         <  1:737188/149‑1 (MQ=255)
                                                            ccaCCCGCCGGGGCACGCGAATGTTCCGGTGCAACGGGCGATGGAGTGGATGCCAATACGCCCGAAAGCCGGTGCTGGCGGGGGCGTCGATCCACAATTGCTCCCGGCGTTTATCGCTATCGAATCGGGGGGTAATCCGAAAGCGGTGa                                                                                         <  2:256628/147‑1 (MQ=255)
                                                                      tGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCTGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTAACGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGaa                                                                               >  1:7354/1‑149 (MQ=255)
                                                                      tGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAATCCGGTGAAGCCTGGGGCGTCGATCCACAATTGATAACGGCGATTATCGCTATCGAATCGGGTGGTAATGCGAACGCGGTGAGTAAATCGaa                                                                               <  1:736792/149‑1 (MQ=255)
                                                                      tGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGaa                                                                               >  2:130689/1‑149 (MQ=255)
                                                                        gAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATg                                                                             >  1:854821/1‑149 (MQ=255)
                                                                           cgcgAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCa                                                                          <  1:808441/149‑1 (MQ=255)
                                                                            gcgAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATAATCGCTATCGAATCGGGTGGTAATCAGAACGCGGAGTGTAAATCGAATCCCAt                                                                         >  2:109937/1‑149 (MQ=255)
                                                                             cgAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCATCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCTAATGCCAtt                                                                        >  1:307066/1‑149 (MQ=255)
                                                                              gAAAGTTCCGGTGCAACGTGCGATGCAGGGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTg                                                                       <  2:610993/149‑1 (MQ=255)
                                                                                aaGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGt                                                                     <  1:919711/149‑1 (MQ=255)
                                                                                aaGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGt                                                                     <  2:854821/149‑1 (MQ=255)
                                                                                         tGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGATAGCCGGTGCAGCATGGGGCGTCGATCCACAATTGATCACGGCGAATATCGCTATCGAATCGGGTGGTAATACGATCGCTGTGAGTTAATCGAATGCCATTGGTTTGATGCAg                                                            >  1:259495/1‑149 (MQ=255)
                                                                                           cAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGtt                                                          <  2:318884/149‑1 (MQ=255)
                                                                                               gTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTaaaa                                                      <  2:183185/149‑1 (MQ=255)
                                                                                                 gCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGc                                                    >  1:450637/1‑149 (MQ=255)
                                                                                                       cAGTGTATGCCAATAAGCCAGAAAGCCGGGGCAGCCTGGTGCGTCGATCCACAATTGATCACGGCGAATATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCCAATGCCATTGGTTTGATCCAGTTAAAAGCTTCAAc                                              <  1:859819/149‑1 (MQ=255)
                                                                                                        aGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAAcc                                             >  1:931854/1‑149 (MQ=255)
                                                                                                          tGGATGCCAATACGCCAGAAAGCCGGTGCAGCCTGGGTCGTCGATCCAAAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCTTTGGTTTGATGCAGTTAAAAGCTTCAACCTc                                           <  2:934117/149‑1 (MQ=255)
                                                                                                                ccAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACg                                     >  2:830448/1‑149 (MQ=255)
                                                                                                                 cAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGt                                    >  2:186543/1‑149 (MQ=255)
                                                                                                                 aaaTAAGCCAGATAGCCGGTGCAGCCTGGGGCGTCGATCGACAATTGATCACGGCGTTTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGt                                    <  1:41413/148‑1 (MQ=255)
                                                                                                                  aaTAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTg                                   >  1:904653/1‑149 (MQ=255)
                                                                                                                         cAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAAACCGAACGCGGTGAGTAAATCGAATGCCATTGGATTGAAGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTa                            <  1:847021/149‑1 (MQ=255)
                                                                                                                           tataGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGTAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATc                          <  1:4899/146‑1 (MQ=255)
                                                                                                                              aGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGGGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGcc                       <  2:237789/149‑1 (MQ=255)
                                                                                                                               gCCGGTGCAGCCTGGGGCGTTGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGCCGTGATGTTTATCGCCg                      >  2:214687/1‑149 (MQ=255)
                                                                                                                               gCCGGTGCAGCCTGGGGCGTTGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTATGATGCAGTTAAAACCTTCAACCTCCGGACGTGATGTTTATCGCCg                      >  2:352939/1‑149 (MQ=255)
                                                                                                                                   gggccgCCTGGGGCGTCGATCCACAATTGATCACGTCGATTATCCCTAGCGACTCCGGTGGTAACCCGAACGCGGTGTGTAACTCGATTGCCTTTGGTTAGATGCAGTTAAAAGCTTCAACCTCCGGACGTTATGTTTATCGCCGTATg                  <  1:971020/144‑1 (MQ=255)
                                                                                                                                   gTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATg                  >  1:321944/1‑149 (MQ=255)
                                                                                                                                     gCAGCCTGGGGCGTCGATCCACAATTGATCACGGCCATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATggg                >  1:522599/1‑149 (MQ=255)
                                                                                                                                      cAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGc               <  2:904653/149‑1 (MQ=255)
                                                                                                                                        gCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTg             <  1:241186/149‑1 (MQ=255)
                                                                                                                                          cTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAAACGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGa           <  1:876400/149‑1 (MQ=255)
                                                                                                                                             gggCGTTGATCCACAATTGATCACGGCGATTATCGCTTTCTAATCGGGGGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGGTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAgtg        <  1:352939/149‑1 (MQ=255)
                                                                                                                                             gggCGTTGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAgtg        <  1:214687/149‑1 (MQ=255)
                                                                                                                                              ggCGTCGCTCCACAATTGATCTCGGCGCTTCTCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAACAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAgtgg       <  2:693396/149‑1 (MQ=255)
                                                                                                                                                cGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGGGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAgtggtg     <  1:867001/149‑1 (MQ=255)
                                                                                                                                                   cGATCCACAATTGATCACGGCGATTTTCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAGTGGTGAGc  <  2:100085/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTGGTTTGATGCAGTTAAAAGCTTCAACCTCCGGACGTGATGTTTATCGCCGTATGGGCTGGAGTGGTGAGC  >  NC_000913/1243188‑1243483

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