BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,115,116	+C	coding (2013/2028 nt)	b2974 ←	putative endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,115,115	1	.	C	98.8%	535.1 / ‑3.5	165	K672E (AAA→GAA)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (69/93); minor base A (2/0); total (71/94)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 180 total aligned reads displayed.

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

gTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                                          >  2:2672253/1‑93 (MQ=255)

 tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑c                                                                                        >  1:976730/1‑93 (MQ=255)

 tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑c                                                                                        <  1:2936799/93‑1 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       <  2:1311393/93‑2 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       <  1:2046163/93‑2 (MQ=255)

    gCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctg                                                                                     <  2:4555925/93‑4 (MQ=255)

     cTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctgc                                                                                    <  1:991882/93‑5 (MQ=255)

       ggTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCgg                                                                                   >  1:3138615/1‑93 (MQ=255)

       ggTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCgg                                                                                   >  1:3893703/1‑93 (MQ=255)

        gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGa                                                                                  <  1:134507/93‑1 (MQ=255)

        gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGa                                                                                  <  2:4360420/93‑1 (MQ=255)

        gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGa                                                                                  <  1:2518277/93‑1 (MQ=255)

         tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGAt                                                                                 <  1:3501349/93‑1 (MQ=255)

         tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATAATCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGAt                                                                                 <  1:802098/93‑1 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  1:759506/1‑93 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  1:3402707/1‑93 (MQ=255)

            cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCga                                                                              <  2:2403572/93‑1 (MQ=255)

            cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATGCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCga                                                                              <  1:1518198/93‑1 (MQ=255)

              gCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgaga                                                                            <  1:4089951/93‑1 (MQ=255)

               cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAt                                                                           <  2:4715215/93‑1 (MQ=255)

               cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAt                                                                           <  2:3427796/93‑1 (MQ=255)

                 cGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCg                                                                         <  1:1434749/93‑1 (MQ=255)

                  gCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGc                                                                        <  2:752423/93‑1 (MQ=255)

                  gCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGc                                                                        >  1:3434780/1‑93 (MQ=255)

                    aGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa                                                                      >  2:568598/1‑93 (MQ=255)

                    aGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa                                                                      <  2:3289288/93‑1 (MQ=255)

                     gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCAt                                                                     >  1:4461199/1‑93 (MQ=255)

                     gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCAt                                                                     <  1:502445/93‑1 (MQ=255)

                     gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCAt                                                                     <  2:412630/93‑1 (MQ=255)

                      cGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATa                                                                    <  2:1928509/93‑1 (MQ=255)

                      cGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATa                                                                    >  2:3846951/1‑93 (MQ=255)

                        ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                                  >  1:4022461/1‑93 (MQ=255)

                        ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                                  >  2:4182495/1‑93 (MQ=255)

                         gATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCt                                                                 <  2:4631973/93‑1 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                >  2:2763602/1‑93 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                <  2:4619106/93‑1 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                <  1:3246436/93‑1 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                >  1:3151976/1‑93 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                <  1:13062/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               >  2:1825415/1‑93 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               >  2:4162803/1‑93 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  2:2868068/93‑1 (MQ=255)

                                tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  1:1531529/1‑93 (MQ=255)

                                 aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                                         <  1:1640843/93‑1 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       >  2:2107289/1‑93 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       <  2:605788/93‑1 (MQ=255)

                                    gATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                                      <  1:363675/93‑1 (MQ=255)

                                     aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     <  2:1047600/93‑1 (MQ=255)

                                      tCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCa                                                    >  1:885342/1‑93 (MQ=255)

                                         ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                                 <  2:3616089/93‑1 (MQ=255)

                                          tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                <  1:3092448/93‑1 (MQ=255)

                                          tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                <  1:3092447/93‑1 (MQ=255)

                                           aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               <  1:2080030/93‑1 (MQ=255)

                                            cccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACagag                                              <  2:943676/93‑1 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             >  2:1561791/1‑93 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             >  1:1567553/1‑93 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             >  1:4743197/1‑93 (MQ=255)

                                              cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            <  2:1304832/93‑1 (MQ=255)

                                              cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            <  1:3466314/93‑1 (MQ=255)

                                                ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGc                                          >  1:2174537/1‑93 (MQ=255)

                                                 cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         <  1:2794841/93‑1 (MQ=255)

                                                 cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         >  2:4100069/1‑93 (MQ=255)

                                                     gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCa                                     >  1:370877/1‑93 (MQ=255)

                                                        acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  2:1215371/1‑93 (MQ=255)

                                                         caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 >  1:3751194/1‑93 (MQ=255)

                                                         caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 <  1:1897434/93‑1 (MQ=255)

                                                         caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 <  2:293297/93‑1 (MQ=255)

                                                          atatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATc                                >  2:1454799/1‑93 (MQ=255)

                                                           tatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               >  2:1096357/1‑93 (MQ=255)

                                                             ttCAAAGCCGTTAAGGATCCTCAGCGGCATTTCCTGCGGATCGAGATCGCTATGGGTAGCTAACGGCTAAACAGAAGTGCCCGCAAAAtacag                             >  1:2314516/1‑89 (MQ=255)

                                                             ttCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCgg                             >  1:2314515/1‑93 (MQ=255)

                                                              tCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGa                            <  2:1715151/93‑1 (MQ=255)

                                                              tCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGa                            <  1:2137977/93‑1 (MQ=255)

                                                               cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt                           <  2:4600865/93‑1 (MQ=255)

                                                               cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt                           <  2:3060393/93‑1 (MQ=255)

                                                                  aGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        >  2:2151786/1‑93 (MQ=255)

                                                                    ccGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg                      <  2:3071561/93‑1 (MQ=255)

                                                                    ccGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg                      <  2:3327250/93‑1 (MQ=255)

                                                                     cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     >  2:1946505/1‑93 (MQ=255)

                                                                      gTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    >  1:2420992/1‑93 (MQ=255)

                                                                      gTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    >  1:3346993/1‑93 (MQ=255)

                                                                      gTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    >  2:4309551/1‑93 (MQ=255)

                                                                        tAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  >  1:1044103/1‑93 (MQ=255)

                                                                         aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  1:418456/93‑1 (MQ=255)

                                                                         aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 >  1:1025715/1‑93 (MQ=255)

                                                                          aGGATCAGCAGCGGCATTTACTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAAc                >  2:1090790/1‑93 (MQ=255)

                                                                            gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGc              <  2:867733/93‑1 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  2:3336183/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  2:4607864/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  2:4514098/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  2:1177541/1‑93 (MQ=255)

                                                                               cagcagCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgc           <  1:1059543/93‑1 (MQ=255)

                                                                                agcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcg          >  1:2774942/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         <  1:3255224/93‑1 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         >  2:2186231/1‑93 (MQ=255)

                                                                                   agcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTg       >  1:4494138/1‑93 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:2666994/93‑1 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:2666995/93‑1 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  2:3748143/1‑93 (MQ=255)

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  1:241081/93‑1 (MQ=255)

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`