BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,180,239	A→C	K166Q (AAA→CAA)	ygiB →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,180,239	0	A	C	97.6%	555.6 / 0.3	165	K166Q (AAA→CAA)	ygiB	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base A (2/2); new base C (84/77); total (86/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 169 total aligned reads displayed.

TCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTAAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG  >  NC_000913/3180142‑3180336

tCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTcaa                                                                                                 <  1:1673802/100‑1 (MQ=255)

 cTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTcaaa                                                                                                >  2:3342363/1‑100 (MQ=255)

 cTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTcaaa                                                                                                >  2:1212633/1‑100 (MQ=255)

 cTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTcaaa                                                                                                >  1:2702875/1‑100 (MQ=255)

   gATGGGCGGCGGGGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACAGACGCGACGGGGCAAAAc                                                                                              <  1:3406007/100‑1 (MQ=255)

    aTGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACt                                                                                             >  1:2931029/1‑100 (MQ=255)

    aTGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACt                                                                                             >  1:2746744/1‑100 (MQ=255)

      gggCGGCGGGCCGGGAGTTGCACAGCGGCCGCTGTTCGCCGCGAGAAACCCAGCCAGTCCGGCTTACGGTGAATATACCGACGCGACGGGTCAAAGCTAt                                                                                           <  2:1330293/100‑1 (MQ=255)

      gggCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTAt                                                                                           >  1:1198341/1‑100 (MQ=255)

       ggcggcggcGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATg                                                                                          <  2:3369368/100‑1 (MQ=255)

       ggcggcggcGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATg                                                                                          <  2:2890131/100‑1 (MQ=255)

       ggcggcggcGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATg                                                                                          <  1:640068/100‑1 (MQ=255)

        gcggcggcgCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATgg                                                                                         >  1:2954872/1‑100 (MQ=255)

        gcggcggcgCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATgg                                                                                         <  1:926935/100‑1 (MQ=255)

         cggcggcgCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGgc                                                                                        <  1:607782/100‑1 (MQ=255)

               cgcgGGATTTGCACAGCAGCCGCAGAAAAAAAAAAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcc                                                                                  <  2:1591338/100‑1 (MQ=255)

                 cgGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagccca                                                                                >  1:71159/1‑100 (MQ=255)

                  gggATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCGGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcccag                                                                               <  1:404340/100‑1 (MQ=255)

                  gggATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcccag                                                                               >  1:678590/1‑100 (MQ=255)

                  gggATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcccag                                                                               >  2:287702/1‑100 (MQ=255)

                   ggATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcccagc                                                                              >  1:2346744/1‑100 (MQ=255)

                    gATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcccagcc                                                                             >  1:297698/1‑100 (MQ=255)

                    gATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGcagcccagcc                                                                             <  2:2644448/100‑1 (MQ=255)

                     aTTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCccagcca                                                                            >  1:1052833/1‑100 (MQ=255)

                       ttGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAgg                                                                          <  1:4034002/100‑1 (MQ=255)

                        tGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGc                                                                         <  1:756677/100‑1 (MQ=255)

                         gCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGcc                                                                        <  1:3778810/100‑1 (MQ=255)

                         gCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGcc                                                                        >  2:774586/1‑100 (MQ=255)

                           acaGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGc                                                                      >  1:2594087/1‑100 (MQ=255)

                              gcagcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCAcc                                                                   >  1:58178/1‑100 (MQ=255)

                              gcagcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCAcc                                                                   <  2:1796834/100‑1 (MQ=255)

                              gcagcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGCCAAAACTATGGCGCACCCCAGCCAGGCCGCAcc                                                                   >  2:2312764/1‑100 (MQ=255)

                               cagcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCa                                                                  <  1:99898/100‑1 (MQ=255)

                               cagcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCa                                                                  >  1:55768/1‑100 (MQ=255)

                                 gcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGCCGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATg                                                                <  2:3451071/100‑1 (MQ=255)

                                 gcCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATg                                                                <  1:3592518/100‑1 (MQ=255)

                                  ccGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGa                                                               >  1:3391311/1‑100 (MQ=255)

                                   cGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGAc                                                              <  1:2994747/100‑1 (MQ=255)

                                   cGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGAc                                                              <  1:872422/100‑1 (MQ=255)

                                     cTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCg                                                            >  2:2073981/1‑100 (MQ=255)

                                          ctcctcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCg                                                       >  1:2886335/1‑100 (MQ=255)

                                          ctcctcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCg                                                       >  2:617858/1‑100 (MQ=255)

                                          ctcctcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCg                                                       >  1:1959887/1‑100 (MQ=255)

                                            atctgAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGaa                                                     <  1:4004329/96‑1 (MQ=255)

                                             ctcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGaan                                                    <  1:3240483/100‑2 (MQ=255)

                                             ctcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAg                                                    >  1:107678/1‑100 (MQ=255)

                                             ctcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAg                                                    <  1:2479295/100‑1 (MQ=255)

                                              tcGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGa                                                   >  2:2579355/1‑100 (MQ=255)

                                               cGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGAc                                                  >  1:2788760/1‑100 (MQ=255)

                                               cGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACACGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGAc                                                  >  1:2404884/1‑100 (MQ=255)

                                                gAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGGCAGGCCGCACCATGACCGTACCGAAGACg                                                 >  2:420824/1‑100 (MQ=255)

                                                gAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGaaagat                                                 >  2:3921441/1‑96 (MQ=255)

                                                gAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACg                                                 >  1:3169773/1‑100 (MQ=255)

                                                  aaaaCCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGc                                               <  1:2455745/100‑1 (MQ=255)

                                                    aaCCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCaa                                             >  2:2637702/1‑100 (MQ=255)

                                                    aaCCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCaa                                             >  2:2637277/1‑100 (MQ=255)

                                                     aCCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAAt                                            >  1:3211306/1‑100 (MQ=255)

                                                     aCCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAAt                                            <  1:752187/100‑1 (MQ=255)

                                                       ccagccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATgg                                          >  1:3641739/1‑100 (MQ=255)

                                                       ccagccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATgg                                          >  1:541094/1‑100 (MQ=255)

                                                        cagccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGc                                         <  1:2370376/100‑1 (MQ=255)

                                                        cagccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGc                                         <  2:2418094/100‑1 (MQ=255)

                                                        cagccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGc                                         <  1:2358453/100‑1 (MQ=255)

                                                         agccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCa                                        <  2:57176/100‑1 (MQ=255)

                                                          gccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCAc                                       <  2:1406124/100‑1 (MQ=255)

                                                           ccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCAcc                                      <  1:2110299/100‑1 (MQ=255)

                                                           ccagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCAcc                                      <  1:2011229/100‑1 (MQ=255)

                                                            cagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCa                                     >  1:3587856/1‑100 (MQ=255)

                                                            cagTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCa                                     <  1:2296597/100‑1 (MQ=255)

                                                            cagTCCGGCTGACGGGAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAAGGGCACCa                                     >  2:3030804/1‑100 (MQ=255)

                                                             agTCCGGCTTACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCaa                                    >  1:1757164/1‑100 (MQ=255)

                                                                     ttACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCg                            <  2:2136258/100‑1 (MQ=255)

                                                                      tACGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGa                           <  1:2139489/100‑1 (MQ=255)

                                                                      tACGGTAAATATACCGACGCGACGGGTAAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGCAGACGGCTCTGGCACCAAAACCGGCGa                           >  2:3575887/1‑100 (MQ=255)

                                                                       aCGGTAAATATACCGACGCGACGGGTCAAAACGATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGac                          >  2:3597319/1‑100 (MQ=255)

                                                                        cGGTACATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCCGACAGGCCGCACCATGACCGTACCGAAGACGGCACTGGCACCAAAACCGGCGacc                         >  2:1102440/1‑100 (MQ=255)

                                                                        cGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGacc                         >  2:1636895/1‑100 (MQ=255)

                                                                        cGGTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGacc                         <  1:2254675/100‑1 (MQ=255)

                                                                         ggTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGacca                        >  1:71926/1‑100 (MQ=255)

                                                                          gTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGaccac                       >  1:588976/1‑100 (MQ=255)

                                                                          gTAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGaccac                       <  1:2649355/100‑1 (MQ=255)

                                                                           tAAATATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGaccacc                      <  1:1498109/100‑1 (MQ=255)

                                                                             aaTATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACcaccac                    >  1:1454436/1‑100 (MQ=255)

                                                                             aaTATACCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACcaccac                    >  1:1538555/1‑100 (MQ=255)

                                                                               tataCCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTa                  >  1:1430038/1‑100 (MQ=255)

                                                                                ataCCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTAc                 <  2:1031089/100‑1 (MQ=255)

                                                                                 taCCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTAcc                >  2:2130725/1‑100 (MQ=255)

                                                                                  aCCGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCg               <  2:1280760/100‑1 (MQ=255)

                                                                                   ccGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGt              >  1:1642974/1‑100 (MQ=255)

                                                                                   ccGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGt              <  2:1326245/100‑1 (MQ=255)

                                                                                   ccGACGCGACGGGTCAAAACTATGGCGCAGCCCAGCCAGGCAGCACCATAACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGt              >  1:3874886/1‑100 (MQ=255)

                                                                                    cgacgcgacgGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGtt             >  1:1432017/1‑100 (MQ=255)

                                                                                    cgacgcgacgGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGtt             >  1:1486597/1‑100 (MQ=255)

                                                                                     gacgcgacgGGTCAAAACTATGGCGCAGCCCAGCCAAGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTa            >  2:3750121/1‑100 (MQ=255)

                                                                                       cgcgacgGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTAcc          <  1:3533983/100‑1 (MQ=255)

                                                                                       cgcgacgGGGCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTAcc          <  1:2561994/100‑1 (MQ=255)

                                                                                         cgacgGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCNGCGACCACCACTACCGTTACCCg        >  2:699319/1‑100 (MQ=255)

                                                                                          gacgGGTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGt       >  2:2323100/1‑100 (MQ=255)

                                                                                             gggTCAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGc    <  2:3951922/100‑1 (MQ=255)

                                                                                              ggTCAAAACTATGGCGCAGCCCAGCGAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGATACCCGTGGCg   >  2:3191327/1‑100 (MQ=255)

TCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCAGCCAGTCCGGCTTACGGTAAATATACCGACGCGACGGGTAAAAACTATGGCGCAGCCCAGCCAGGCCGCACCATGACCGTACCGAAGACGGCAATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG  >  NC_000913/3180142‑3180336

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: `‑`