BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,663,474	A→T	15.4%	A11A (GCA→GCT)	suhB →	inositol monophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,663,474	0	A	T	15.4%	262.5 / 38.7	104	A11A (GCA→GCT)	suhB	inositol monophosphatase
Reads supporting (aligned to +/- strand): ref base A (45/43); new base T (8/8); total (53/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 116 total aligned reads displayed.

TGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGA  >  NC_000913/2663378‑2663572

tGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTgcgcgctgcg                                                                                                 <  2:1281824/100‑1 (MQ=255)

tGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcg                                                                                                 >  2:2363234/1‑100 (MQ=255)

  ccGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcg                                                                                               >  1:294787/1‑100 (MQ=255)

      tctcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                           <  1:2787984/100‑1 (MQ=255)

        tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                         <  2:2155911/100‑1 (MQ=255)

        tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGc                                                                                         >  2:217089/1‑100 (MQ=255)

        tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGc                                                                                         >  2:1297529/1‑100 (MQ=255)

         cGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGGACATCGCCGTGCGCGCAGCGCGCAAGGCg                                                                                        >  2:2862612/1‑100 (MQ=255)

          gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                       <  1:1105399/100‑1 (MQ=255)

           ctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                      <  2:2327492/100‑3 (MQ=39)

           ctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                      <  1:1976626/100‑1 (MQ=255)

           ctgctATACTCTGAGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                      >  2:1976637/1‑100 (MQ=255)

            tgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGt                                                                                     <  2:632798/100‑1 (MQ=255)

             gctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgt                                                                                    <  1:2673958/100‑5 (MQ=38)

              ctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCGAGTGAGAGAGACATATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                      <  2:705646/97‑3 (MQ=37)

               tagaCTCTGCGCCGTTTGCCCGTTCTTGAACATCCAGTGAGAGAGACTGATGCATCCGATGCTGAGCATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaat                                                                                  <  2:118122/97‑1 (MQ=255)

                ataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttca                                                                                 <  1:2173940/100‑8 (MQ=37)

                ataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatt                                                                                 >  2:2000876/1‑100 (MQ=255)

                 taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttcag                                                                                <  1:1025516/100‑9 (MQ=37)

                 taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttcag                                                                                <  1:202294/100‑9 (MQ=37)

                 taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttcag                                                                                <  1:1117470/100‑9 (MQ=37)

                 taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattt                                                                                <  1:1571740/100‑1 (MQ=255)

                   ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                              <  2:15299/100‑1 (MQ=255)

                   ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                              >  1:335048/1‑100 (MQ=255)

                    tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                             <  2:1941860/100‑1 (MQ=255)

                      tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                           <  1:1039806/100‑1 (MQ=255)

                      tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                           <  1:1297524/100‑1 (MQ=255)

                       gcgcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGc                                                                          <  2:2607160/100‑1 (MQ=255)

                         gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCa                                                                        >  1:1664354/1‑100 (MQ=255)

                         gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCa                                                                        >  1:437212/1‑100 (MQ=255)

                           cGTTTTCCCGTTCTTTAACATCCAGTGAGCGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaa                                                                       <  1:2230202/99‑1 (MQ=255)

                              tttCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCAGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAc                                                                   >  2:1724966/1‑100 (MQ=255)

                                tCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                 <  2:1032817/100‑1 (MQ=255)

                                 cccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                >  1:1203087/1‑100 (MQ=255)

                                  ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                               >  1:2341219/1‑100 (MQ=255)

                                  ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                               <  2:2642117/100‑1 (MQ=255)

                                  ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCACGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                               >  1:1535777/1‑100 (MQ=255)

                                       cTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTCCCAAAAACTATGaaaac                                                          >  1:1177135/1‑98 (MQ=255)

                                        tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACCTCGCCGTGCGCACAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                         <  1:2044395/100‑1 (MQ=255)

                                        tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                         >  2:1541371/1‑100 (MQ=255)

                                        tttAACATCCAGAGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                         >  2:1920430/1‑100 (MQ=255)

                                           aaCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCgg                                                      <  2:2438469/100‑1 (MQ=255)

                                            aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGa                                                     <  1:910684/100‑1 (MQ=255)

                                             cATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCCAGGCGGGTAATTTAATTGCCCAAAACTATGAAACCCCCGAc                                                    >  1:902451/1‑100 (MQ=255)

                                             cATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGAc                                                    <  1:1753021/100‑1 (MQ=255)

                                                ccAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCt                                                 <  1:2755932/100‑1 (MQ=255)

                                                 cAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                >  1:445652/1‑100 (MQ=255)

                                                 cAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                >  1:2551380/1‑100 (MQ=255)

                                                 cAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                <  1:76935/100‑1 (MQ=255)

                                                  aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGt                                               <  2:2356881/100‑1 (MQ=255)

                                                     gagagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGa                                            >  1:809177/1‑100 (MQ=255)

                                                        agagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGc                                         <  1:431296/100‑1 (MQ=255)

                                                        agagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGc                                         <  1:1945259/100‑1 (MQ=255)

                                                         gagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGCCGCTGTAGAAGCg                                        >  1:872975/1‑100 (MQ=255)

                                                           gaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                      <  1:623485/100‑1 (MQ=255)

                                                            aCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCACAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                       >  1:1723135/1‑98 (MQ=255)

                                                            aCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGc                                     <  2:566605/100‑1 (MQ=255)

                                                             ccGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGcc                                    >  2:1397908/1‑100 (MQ=255)

                                                               gATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAg                                  >  2:1265557/1‑100 (MQ=255)

                                                               gATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAg                                  >  2:1543926/1‑100 (MQ=255)

                                                               gATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAg                                  >  1:192923/1‑100 (MQ=255)

                                                                 tGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaa                                >  1:2534506/1‑100 (MQ=255)

                                                                 tGCATCCGACGCGGAATACCGCTGTGCGCGCAGCGCGCGAGGGGGGTAACCGAATGGCGAAACACGATGAAAACCCGGACACTGTACAAGCGTGCCa                                   <  1:490401/97‑1 (MQ=255)

                                                                   cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAg                              >  1:274972/1‑100 (MQ=255)

                                                                     tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGCCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                            >  1:1553958/1‑100 (MQ=255)

                                                                      ccGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                           >  2:2576794/1‑100 (MQ=255)

                                                                      ccGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                           <  2:2411330/100‑1 (MQ=255)

                                                                       cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                          <  2:2321911/100‑1 (MQ=255)

                                                                       cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                          >  2:25715/1‑100 (MQ=255)

                                                                       cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                          >  1:1787964/1‑100 (MQ=255)

                                                                        gaggCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                          <  1:2461429/96‑1 (MQ=255)

                                                                        gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCAGGTAATTTAATTGCCAAAAACTCTGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGt                         >  1:2397550/1‑100 (MQ=255)

                                                                         aTGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTa                        <  2:1248041/100‑1 (MQ=255)

                                                                           gCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAc                      <  2:2112743/100‑1 (MQ=255)

                                                                           gCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAc                      >  1:794753/1‑100 (MQ=255)

                                                                             tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                    >  2:704565/1‑100 (MQ=255)

                                                                             ngAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGTCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                    >  2:1615073/2‑100 (MQ=255)

                                                                             aattacccGCCTTGCGCGCTGCGCGCAAGGCGGGTAATTTAAATGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                    >  2:1972652/8‑100 (MQ=38)

                                                                              gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                    >  1:2177153/1‑99 (MQ=255)

                                                                                aCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAttt                 <  2:100840/100‑1 (MQ=255)

                                                                                  cccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg               >  2:671277/3‑100 (MQ=38)

                                                                                  cccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg               >  2:804326/3‑100 (MQ=38)

                                                                                  cccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg               >  2:2336537/3‑100 (MQ=38)

                                                                                  aTCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg               <  1:800616/100‑1 (MQ=255)

                                                                                   tCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt              <  2:2363621/100‑1 (MQ=255)

                                                                                   ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt              >  2:780556/2‑100 (MQ=39)

                                                                                    cGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg             >  1:615744/1‑100 (MQ=255)

                                                                                    cGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg             <  2:2341234/100‑1 (MQ=255)

                                                                                     gCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa            >  2:2306048/1‑100 (MQ=255)

                                                                                      ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAc           <  2:2821664/100‑1 (MQ=255)

                                                                                      ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAc           <  2:1553964/100‑1 (MQ=255)

                                                                                      ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAc           <  2:726627/100‑1 (MQ=255)

                                                                                       cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAcc          >  2:2092099/1‑100 (MQ=255)

                                                                                       cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAcc          <  2:2868257/100‑1 (MQ=255)

                                                                                            gcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGt     >  2:2059042/1‑100 (MQ=255)

                                                                                             cgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa    <  2:313957/100‑1 (MQ=255)

                                                                                              gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAg   <  2:1450995/100‑1 (MQ=255)

                                                                                              gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAg   >  2:1932943/1‑100 (MQ=255)

                                                                                              gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAg   <  1:1874531/100‑1 (MQ=255)

                                                                                               cagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGa  >  2:1228563/1‑100 (MQ=255)

TGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGA  >  NC_000913/2663378‑2663572

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