Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,854 A→T 6.9% intergenic (‑141/‑142) yefM ← / → hisL YefM antitoxin of the YoeB‑YefM toxin‑antitoxin pair and DNA binding transcriptional repressor/his operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,8540AT 6.9% 163.6 / 4.9 59intergenic (‑141/‑142)yefM/hisLYefM antitoxin of the YoeB‑YefM toxin‑antitoxin pair and DNA binding transcriptional repressor/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base T (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-02

ATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTA  >  NC_000913/2089712‑2089996
                                                                                                                                              |                                                                                                                                              
attttaaTTTCCTCTTTTGTACAGTGCATTTCACAAAAAAGAGCGTGAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACa                                                                                                                                                <  1:748673/138‑1 (MQ=255)
aTGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACa                                                                                                                                                >  2:577180/1‑143 (MQ=255)
  gTCAATCTCCTCTTTTGGACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACGGt                                                                                                                                              >  1:1026584/1‑143 (MQ=255)
     aaTCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCct                                                                                                                                           >  2:504680/1‑143 (MQ=255)
     aaTCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCATCACTATTTTAGATGCCGTCAGCAGAGAGGACAGTGct                                                                                                                                           >  2:527190/1‑143 (MQ=255)
          cctcTTTTGTACAGCTCATTGTACAATCATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAg                                                                                                                                      <  1:959765/143‑1 (MQ=255)
            tctTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTACGAGGACAGTCCTCTTAgac                                                                                                                                    >  1:1034608/1‑141 (MQ=255)
            tctTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGcc                                                                                                                                    >  2:1029106/1‑143 (MQ=255)
            tctTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAAGTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATAGTTGGTGCCATCAGCTACAAGGACAGTCCTCGTAGcc                                                                                                                                    >  2:915134/1‑143 (MQ=255)
              ttttGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCGTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGctga                                                                                                                                  <  1:814805/143‑4 (MQ=255)
                    aCAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAATTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctccat                                                                                                                            >  2:576335/1‑141 (MQ=255)
                            ttGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCt                                                                                                                    >  1:964710/1‑143 (MQ=255)
                            ttGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCt                                                                                                                    >  1:1071767/1‑143 (MQ=255)
                             tGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtc                                                                                                                   >  2:489992/1‑143 (MQ=255)
                                 cAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattc                                                                                                               >  1:591312/1‑143 (MQ=255)
                                 cAATGATGAGCGTTAATTAACTAGTGATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattc                                                                                                               <  2:1063992/143‑1 (MQ=255)
                                  aaTGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAACATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                                                              <  1:1048301/143‑1 (MQ=255)
                                    tgatgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattcatg                                                                                                            >  2:478787/1‑142 (MQ=255)
                                       tgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTaaa                                                                                                         <  2:1047053/143‑1 (MQ=255)
                                                attaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCAt                                                                                                <  2:42287/143‑1 (MQ=255)
                                                      tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACACATCCATTCCCAt                                                                                          >  2:366386/1‑143 (MQ=255)
                                                       atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAACCAAATCCATTCCCAtc                                                                                         >  2:887330/1‑142 (MQ=255)
                                                        tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaa                                                                                        >  1:531181/1‑143 (MQ=255)
                                                                ttaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAAtctata                                                                                >  2:834822/1‑143 (MQ=255)
                                                                 cagcttGTAGATCAAGGTATTGTCACTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGTACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaa                                                                               <  1:80470/139‑1 (MQ=255)
                                                                    tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaaa                                                                            >  2:878775/1‑143 (MQ=255)
                                                                      tGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAg                                                                          >  2:667054/1‑143 (MQ=255)
                                                                        tAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAACATATATAAAAAAGcc                                                                        >  1:463349/1‑143 (MQ=255)
                                                                         aGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGccc                                                                       >  1:832389/1‑143 (MQ=255)
                                                                          gATCAAGGTATTGTCAGTGAGACGAAAATTCAGGCACCCCTATGTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAAGCCATTGCCATAAAATATATAAAAAAGCCCt                                                                      <  2:1070229/143‑1 (MQ=255)
                                                                                ggTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCATGCttt                                                                >  1:967599/1‑143 (MQ=255)
                                                                                   aTTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTa                                                             <  1:369586/143‑1 (MQ=255)
                                                                                     tggcAGAGAGACGAAAACCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAAc                                                           <  1:438925/140‑1 (MQ=255)
                                                                                     tGTCAGTGAGACGAGTTGCCATGCGTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAAc                                                           <  2:1009690/143‑1 (MQ=255)
                                                                                         aGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCTCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGa                                                       >  1:512512/1‑143 (MQ=255)
                                                                                           tGAGACGAAAACCCCGGCTCCGCGATCTTTTCTGCCAGCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaaa                                                     <  2:267204/143‑1 (MQ=255)
                                                                                               aCGAAAATCCAGGCTGCGCTATTTTTGGTCCCATCAGCTAAGAGGACAGTCCTCTTAACCCCCTCCATTCCCCGCTCATTCATTAAACAAATCCATTGCCAAAACAACTATAAAAAAGCACCTGCTTTCTAACGACAAACTgg                                                 >  2:737063/1‑143 (MQ=255)
                                                                                                 gAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAGATATGTCCAACAGTCCTTGCGCTCTAAGGTGACAgaatgc                                               >  1:32333/1‑138 (MQ=255)
                                                                                                  aaaaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt                                              <  2:156734/143‑1 (MQ=255)
                                                                                                    aaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAACATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAg                                            >  2:830944/1‑143 (MQ=255)
                                                                                                    aaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAg                                            <  1:346202/143‑1 (MQ=255)
                                                                                                     aTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAgg                                           <  1:343070/143‑1 (MQ=255)
                                                                                                      tCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAACTATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTCGGt                                          >  1:343423/1‑143 (MQ=255)
                                                                                                         aGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCGTTCATTAAACAAATCCATTGCCATAAAATCTATAAAAAAGCCCTTGATTTCTAACGTGAAAGTGGCTCAGGTTaa                                       >  2:197635/1‑143 (MQ=255)
                                                                                                           gCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTaaaa                                     <  2:561462/143‑1 (MQ=255)
                                                                                                               cGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACa                                 <  1:870505/143‑1 (MQ=255)
                                                                                                               cGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACATGAAAGTTGTTTAGGTTAAAAGACa                                 >  1:846481/1‑143 (MQ=255)
                                                                                                                gCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTCCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACAt                                >  1:559764/1‑143 (MQ=255)
                                                                                                                     ttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCAGTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGtt                           <  2:105001/143‑1 (MQ=255)
                                                                                                                     ttttGGCGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCCTTCCCCGCTCACTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGtt                           <  2:707995/143‑1 (MQ=255)
                                                                                                                               aTCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACAt                 >  1:570744/1‑143 (MQ=255)
                                                                                                                               aTCAGATAAGAGTACACTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAACCCAATCCATTGCCATAAAATATATAAAAAAGCCCTTGTTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACAt                 <  1:359004/143‑1 (MQ=255)
                                                                                                                                  aGCTATGAGGACAGTCCTCTTAGCCCCCTCCTTTCCGCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAATGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTca              <  2:430032/143‑1 (MQ=255)
                                                                                                                                   gCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTcac             <  1:878775/143‑1 (MQ=255)
                                                                                                                                    cTAAGAGGCCAGGGTGCTTAGCCCCCCCCTTTCACCGCTTATTCATTAAACAAATCCATTGCCATAAAAGATATAAAAAAGCGCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTcaca            <  2:1067333/143‑1 (MQ=255)
                                                                                                                                    cTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCAATCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTcaca            <  2:12558/143‑1 (MQ=255)
                                                                                                                                     tAAGAGGGCTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACag           >  1:1029735/1‑143 (MQ=255)
                                                                                                                                     tAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTAATTAAACAAAACCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACag           >  1:758554/1‑143 (MQ=255)
                                                                                                                                     tAAGAGGACACTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACag           <  1:320031/143‑1 (MQ=255)
                                                                                                                                           gACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACtt     <  1:912336/143‑1 (MQ=255)
                                                                                                                                             cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACtttt   <  1:489992/143‑1 (MQ=255)
                                                                                                                                              aGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATACAATATATAAAACAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCATAGAGAcactca  >  2:653119/1‑138 (MQ=255)
                                                                                                                                              aGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTa  <  2:593079/143‑1 (MQ=255)
                                                                                                                                              aGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAAGTTAAAAGACATCAGTTGAATAAACATTCAAAGAGACTTTTa  >  2:27116/1‑143 (MQ=255)
                                                                                                                                              aGTCCTCATAGCCCCCTCCTTTCCGCGCTCAGTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTa  <  2:47738/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTA  >  NC_000913/2089712‑2089996

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: