BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,494,457	C→T	A149A (GCC→GCT)	nirB →	nitrite reductase catalytic subunit NirB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,494,457	0	C	T	100.0%	947.5 / NA	230	A149A (GCC→GCT)	nirB	nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (113/117); total (113/117)

Only 100 of 232 total aligned reads displayed.

CCAATCAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGG  >  NC_000913/3494359‑3494554

ccAATCAAAGGTTCTGATGCTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTgt                                                                                                 >  2:4561963/1‑101 (MQ=255)

   aTCAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCtgttgt                                                                                              >  1:4561966/1‑101 (MQ=255)

    tCAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTgttgtt                                                                                             <  1:4561478/101‑1 (MQ=255)

    tCAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTgttgtt                                                                                             <  2:4561402/101‑1 (MQ=255)

     cAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGttgttg                                                                                            >  1:4561968/1‑101 (MQ=255)

     cAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGttgttg                                                                                            <  2:4561687/101‑1 (MQ=255)

      aaaGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTtgg                                                                                           >  2:4561970/1‑101 (MQ=255)

       aaGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTtggt                                                                                          >  1:4561973/1‑101 (MQ=255)

       aaGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTtggt                                                                                          >  2:4561974/1‑101 (MQ=255)

          gTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTggc                                                                                       <  2:4561722/101‑1 (MQ=255)

             cTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTggcggc                                                                                    >  1:4561978/1‑101 (MQ=255)

             cTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTggcggc                                                                                    >  2:4561980/1‑101 (MQ=255)

              tGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTtggtgg                                                                                        >  2:4561982/1‑96 (MQ=255)

              tGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGcc                                                                                   <  2:4561739/101‑1 (MQ=255)

               gATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCt                                                                                  <  1:4561556/101‑1 (MQ=255)

                 tACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGt                                                                                >  1:4561985/1‑101 (MQ=255)

                  aCTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGtt                                                                               >  1:4561987/1‑101 (MQ=255)

                   cTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTa                                                                              >  1:4561990/1‑101 (MQ=255)

                    tCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAg                                                                             >  2:4561991/1‑101 (MQ=255)

                    tCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAg                                                                             <  1:4561569/101‑1 (MQ=255)

                     cAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAgg                                                                            >  1:4561993/1‑101 (MQ=255)

                       ggACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTc                                                                          <  2:4561765/101‑1 (MQ=255)

                       ggACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTc                                                                          >  1:4561995/1‑101 (MQ=255)

                        gACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCt                                                                         <  1:4561535/101‑1 (MQ=255)

                         aCTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTg                                                                        <  2:4561553/101‑1 (MQ=255)

                            gCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGaa                                                                     >  1:4561998/1‑101 (MQ=255)

                            gCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGaa                                                                     >  1:4562001/1‑101 (MQ=255)

                             cTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAg                                                                    <  1:4561461/101‑1 (MQ=255)

                             cTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAg                                                                    <  1:4561647/101‑1 (MQ=255)

                              tttGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGc                                                                   >  2:4562004/1‑101 (MQ=255)

                               ttGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCt                                                                         >  1:4562010/1‑94 (MQ=255)

                               ttGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGcc                                                                  >  1:4562006/1‑101 (MQ=255)

                               ttGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGACTGTTAGGTCTGGAAGcc                                                                  >  2:4562007/1‑101 (MQ=255)

                                  tCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCa                                                               <  2:4561606/101‑1 (MQ=255)

                                  tCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCa                                                               <  1:4561956/101‑1 (MQ=255)

                                   cTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAg                                                              <  1:4561721/101‑1 (MQ=255)

                                   cTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAg                                                              >  2:4562013/1‑101 (MQ=255)

                                    tATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAgg                                                             <  1:4561682/101‑1 (MQ=255)

                                    tATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAgg                                                             <  1:4561665/101‑1 (MQ=255)

                                     aTCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGgc                                                            <  2:4561568/101‑1 (MQ=255)

                                        gCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGgcgcg                                                         >  1:4562017/1‑101 (MQ=255)

                                         cACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTggcg                                                                                      >  2:4562020/1‑71 (MQ=255)

                                         cACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGgcgcgc                                                        <  1:4561497/101‑1 (MQ=255)

                                          aCTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCt                                                       >  2:4562023/1‑101 (MQ=255)

                                          aCTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCt                                                       >  1:4562022/1‑101 (MQ=255)

                                            tATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGa                                                     >  2:4562025/1‑101 (MQ=255)

                                            tATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGa                                                     <  2:4561476/101‑1 (MQ=255)

                                            tATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGa                                                     <  1:4561576/101‑1 (MQ=255)

                                              ttGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGaaa                                                   <  1:4561784/101‑1 (MQ=255)

                                               tGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGaaaa                                                  <  2:4561949/101‑1 (MQ=255)

                                               tGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGaaaa                                                  >  1:4562029/1‑101 (MQ=255)

                                                 aaGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAAc                                                <  2:4561588/101‑1 (MQ=255)

                                                   gtccTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAg                                                              <  1:4561615/83‑1 (MQ=255)

                                                   gACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACtt                                              <  1:4561523/101‑1 (MQ=255)

                                                    aCCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTa                                             <  1:4561620/101‑1 (MQ=255)

                                                    aCCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTa                                             >  1:4562032/1‑101 (MQ=255)

                                                    aCCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTa                                             >  2:4562033/1‑101 (MQ=255)

                                                      cTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAgg                                           <  2:4561771/101‑1 (MQ=255)

                                                      cTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAgg                                           <  1:4561842/101‑1 (MQ=255)

                                                       tCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGt                                          >  1:4562036/1‑101 (MQ=255)

                                                       tCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGt                                          >  2:4562037/1‑101 (MQ=255)

                                                           cGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTg                                      >  2:4562040/1‑101 (MQ=255)

                                                            gCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGgcgc                                                          <  2:4561583/80‑1 (MQ=255)

                                                            gCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGa                                     >  1:4562042/1‑101 (MQ=255)

                                                            gCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGa                                     <  2:4561752/101‑1 (MQ=255)

                                                             ccATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGaa                                    <  1:4561650/101‑1 (MQ=255)

                                                              cATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGaaa                                   <  1:4561693/101‑1 (MQ=255)

                                                              cATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGaaa                                   >  1:4562045/1‑101 (MQ=255)

                                                               aTTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAAc                                  <  2:4561708/101‑1 (MQ=255)

                                                                 tGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAAccc                                >  1:4562049/1‑101 (MQ=255)

                                                                 tGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAAccc                                >  2:4562048/1‑101 (MQ=255)

                                                                  gAATCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCa                               >  2:4562051/1‑101 (MQ=255)

                                                                   aaTCCTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCAc                              <  1:4561767/101‑1 (MQ=255)

                                                                      ccTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGtt                           <  2:4561713/101‑1 (MQ=255)

                                                                       cTGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTa                          <  1:4561704/101‑1 (MQ=255)

                                                                        tGCGCCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTAt                         <  1:4561803/101‑1 (MQ=255)

                                                                         gcgcCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATc                        >  2:4562055/1‑101 (MQ=255)

                                                                          cgcCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCg                       <  2:4561683/101‑1 (MQ=255)

                                                                           gcCCGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACg                             >  2:4562058/1‑94 (MQ=255)

                                                                             ccGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAAt                    <  1:4561904/101‑1 (MQ=255)

                                                                              cGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAAtt                   >  2:4562060/1‑101 (MQ=255)

                                                                              cGTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAAtt                   <  1:4561511/101‑1 (MQ=255)

                                                                               gTCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAAttt                  <  1:4561717/101‑1 (MQ=255)

                                                                                tCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTg                 <  2:4561452/101‑1 (MQ=255)

                                                                                tCGCAGCAAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTg                 <  2:4561745/101‑1 (MQ=255)

                                                                                   cagcaAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGccc              >  2:4562062/1‑101 (MQ=255)

                                                                                    agcaAACGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGcccc             >  2:4562065/1‑101 (MQ=255)

                                                                                      caAACGCGGTGCTGTTGTTGGTGGCGGACTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTa           >  1:4562067/1‑101 (MQ=255)

                                                                                       aaaCGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTAt          <  1:4561594/101‑1 (MQ=255)

                                                                                       aaaCGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTAt          >  2:4562068/1‑101 (MQ=255)

                                                                                        aaCGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATg         <  1:4561625/101‑1 (MQ=255)

                                                                                         aCGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCt            <  2:4561699/97‑1 (MQ=255)

                                                                                         aCGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGc        <  2:4561663/101‑1 (MQ=255)

                                                                                         aCGCGGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGc        >  1:4562070/1‑101 (MQ=255)

                                                                                           gcgGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTg      >  2:4562073/1‑101 (MQ=255)

                                                                                            cgGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGa     <  2:4561559/101‑1 (MQ=255)

                                                                                            cgGTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGa     >  1:4562074/1‑101 (MQ=255)

                                                                                             ggTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGAt    <  1:4561746/101‑1 (MQ=255)

                                                                                              gTGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATg   >  1:4562076/1‑101 (MQ=255)

                                                                                               tGCTGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATgg  <  1:4561728/101‑1 (MQ=255)

CCAATCAAAGGTTCTGATACTCAGGACTGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGG  >  NC_000913/3494359‑3494554

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