Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,763,869 (T)7→8 coding (1285/1617 nt) yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,763,8621.T95.3% 214.7 / 1.6 64coding (1292/1617 nt)yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base T (29/32);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑T‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCAC  >  NC_000913/2763732‑2763988
                                                                                                                                   |                                                                                                                                 
tagATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑t‑tt‑tt‑ttt                                                                                                                        <  1:896321/139‑2 (MQ=255)
     tAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTgcgcg                                                                                                                    <  1:652180/139‑1 (MQ=255)
      agagTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTgcgcgc                                                                                                                   >  2:1155939/1‑139 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTgcgcgc                                                                                                                   <  2:467015/137‑1 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTgcgcgc                                                                                                                   >  1:467015/1‑137 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑t‑tt‑tt‑tttgc                                                                                                                      >  2:1126293/1‑130 (MQ=255)
        agTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑t‑tt‑tt‑tttgc                                                                                                                      <  1:1126293/133‑4 (MQ=255)
          tCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTa                                                                                                               >  1:353309/1‑139 (MQ=255)
             tAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCg                                                                                                            >  2:253986/1‑139 (MQ=255)
                cccACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTTAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCTAGGA‑TATT‑TT‑TTGCGCGCGGTATCGATg                                                                                                         <  2:816721/139‑1 (MQ=255)
                  cACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCa                                                                                                       >  1:169543/1‑139 (MQ=255)
                  cACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGAATCGATGCa                                                                                                       >  1:776327/1‑139 (MQ=255)
                          gaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGaa                                                                                                >  1:739876/1‑138 (MQ=255)
                          gaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGaa                                                                                                <  2:739876/138‑1 (MQ=255)
                           aacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTgcgcg                                                                                                                    >  1:882924/1‑117 (MQ=255)
                           aacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTAt                                                                                                              <  1:487148/123‑1 (MQ=255)
                           aacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTAt                                                                                                              >  2:487148/1‑123 (MQ=255)
                           aacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑T‑TTGTT‑TTgcgcg                                                                                                                    <  2:882924/117‑1 (MQ=255)
                            acaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATg                                                                                             <  1:551535/139‑1 (MQ=255)
                             caAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                            <  1:253986/139‑1 (MQ=255)
                               aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCCATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          <  1:1227105/139‑1 (MQ=255)
                               aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          <  1:991167/139‑1 (MQ=255)
                               aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          <  2:900721/139‑1 (MQ=255)
                                aCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCt                                                                                         >  1:1184944/1‑139 (MQ=255)
                                   tGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGt                                                                                      <  2:555652/139‑1 (MQ=255)
                                     aaaCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTAc                                                                                    <  2:391256/139‑1 (MQ=255)
                                           tCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                            >  1:526802/1‑125 (MQ=255)
                                           tCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGt                                                                                            <  2:526802/125‑1 (MQ=255)
                                                     cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTAtct                                                                    >  2:348927/1‑139 (MQ=255)
                                                     cATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGAATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTAtct                                                                    >  1:941987/1‑139 (MQ=255)
                                                         gAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCgg                                                                <  2:634670/139‑1 (MQ=255)
                                                               ttGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCt                                                          >  1:639916/1‑139 (MQ=255)
                                                                 ggCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          <  1:328686/105‑1 (MQ=255)
                                                                 ggCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          >  2:328686/1‑105 (MQ=255)
                                                                    tGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACgg                                                     <  2:1087977/139‑1 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                    <  2:1087898/139‑1 (MQ=255)
                                                                     gCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                    <  2:1184944/139‑1 (MQ=255)
                                                                      ccATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                   <  2:754761/139‑1 (MQ=255)
                                                                                tttCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAatat                                         <  1:818224/139‑1 (MQ=255)
                                                                                        aTGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCAcc                                 >  1:1316361/1‑139 (MQ=255)
                                                                                          gTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCaa                               <  2:639916/139‑1 (MQ=255)
                                                                                           tCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑T‑TT‑TTATTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAAc                              <  1:934090/139‑1 (MQ=255)
                                                                                                              tATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          >  1:480082/1‑60 (MQ=255)
                                                                                                              tATCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                                                                          <  2:480082/60‑1 (MQ=255)
                                                                                                               aTCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAAt                                                 >  1:309087/1‑100 (MQ=255)
                                                                                                               aTCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAAt                                                 <  2:309087/100‑1 (MQ=255)
                                                                                                               aTCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAg          >  2:138559/1‑139 (MQ=255)
                                                                                                                tCGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGc         >  2:1070845/1‑139 (MQ=255)
                                                                                                                 cGTCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCt        >  1:375859/1‑139 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                    >  2:438979/1‑93 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGc                                                    <  1:438979/93‑1 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCaaa           <  1:513547/134‑1 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCaaa           >  2:513547/1‑134 (MQ=255)
                                                                                                                   tCTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGa      >  1:60894/1‑139 (MQ=255)
                                                                                                                    cTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGNAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt     <  2:702822/139‑1 (MQ=255)
                                                                                                                    cTTTAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGAt     <  2:941987/139‑1 (MQ=255)
                                                                                                                     tttAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc    <  1:663853/139‑1 (MQ=255)
                                                                                                                     tttAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc    <  2:375859/139‑1 (MQ=255)
                                                                                                                       tAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                                                                                 >  2:399006/1‑60 (MQ=255)
                                                                                                                       tAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                                                                                 <  1:399006/60‑1 (MQ=255)
                                                                                                                       tAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                   <  1:911384/90‑1 (MQ=255)
                                                                                                                       tAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                   >  2:911384/1‑90 (MQ=255)
                                                                                                                       tAATGTCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCAc  >  2:485902/1‑139 (MQ=255)
                                                                                                                            tCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                   >  2:1284583/1‑85 (MQ=255)
                                                                                                                            tCGAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCa                                                   <  1:1284583/85‑1 (MQ=255)
                                                                                                                              gAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc    <  2:231533/130‑1 (MQ=255)
                                                                                                                              gAGGATT‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATc    >  1:231533/1‑130 (MQ=255)
                                                                                                                                   |                                                                                                                                 
TAGATTAGAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑T‑TT‑TT‑TTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCTGACGGCAATGTTAATATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCAC  >  NC_000913/2763732‑2763988

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