BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,011,749	A→G	A369A (GCA→GCG)	metE →	tetrahydropteroyltriglutamate methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,011,749	0	A	G	100.0%	980.8 / NA	313	A369A (GCA→GCG)	metE	tetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (150/163); total (150/163)

Only 100 of 342 total aligned reads displayed.

TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838

 gCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGct                                                                                         <  1:1857413/93‑1 (MQ=255)

  cAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctc                                                                                        <  1:4364837/93‑1 (MQ=255)

   aGAAGTGAAAAGCTGGTTTGCCTTCGCTCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctct                                                                                       <  2:3214488/93‑1 (MQ=255)

    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCCGAACAGTGGTGACACGGCGGCTCTg                                                                                      <  1:6361820/93‑1 (MQ=255)

     aaGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTgg                                                                                     <  2:4215166/93‑1 (MQ=255)

      aGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGc                                                                                    >  1:3296231/1‑93 (MQ=255)

       gTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCa                                                                                   <  2:100705/93‑1 (MQ=255)

          aaaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCagag                                                                                <  2:7324413/93‑1 (MQ=255)

          aaaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCagag                                                                                <  2:5186796/93‑1 (MQ=255)

              gCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGa                                                                            <  2:4887124/93‑1 (MQ=255)

               cTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAg                                                                           <  1:148539/93‑1 (MQ=255)

               cTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAg                                                                           <  2:2976559/93‑1 (MQ=255)

                   tttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑cc                                                                      >  2:3130980/1‑93 (MQ=255)

                      gCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCg                                                                   >  2:1063673/1‑93 (MQ=255)

                      gCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCg                                                                   >  1:3791753/1‑93 (MQ=255)

                       ccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGa                                                                  <  2:7574976/93‑1 (MQ=255)

                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGAtt                                                                >  1:5351967/1‑93 (MQ=255)

                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGAtt                                                                >  1:4497145/1‑93 (MQ=255)

                          tCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTc                                                               <  2:5808333/93‑1 (MQ=255)

                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCa                                                              <  1:6277323/93‑1 (MQ=255)

                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAg                                                             >  1:6734063/1‑93 (MQ=255)

                              ccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGc                                                           >  1:287284/1‑93 (MQ=255)

                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGcc                                                          >  1:355075/1‑92 (MQ=255)

                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCa                                                          <  1:1443866/93‑1 (MQ=255)

                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCa                                                          <  1:2437389/93‑1 (MQ=255)

                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCa                                                          <  2:7199287/93‑1 (MQ=255)

                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAc                                                         >  2:769630/1‑93 (MQ=255)

                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAc                                                         >  1:6474819/1‑93 (MQ=255)

                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAc                                                         >  1:73619/1‑93 (MQ=255)

                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAc                                                         >  1:1564499/1‑93 (MQ=255)

                                  cAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAcgt                                                       >  2:1519071/1‑93 (MQ=255)

                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCCcgtc                                                      >  2:1278373/1‑93 (MQ=255)

                                    aaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAcgtcg                                                     >  2:2735750/1‑93 (MQ=255)

                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCAcgtcgc                                                    >  1:800632/1‑92 (MQ=255)

                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACgtcgtc                                                   >  2:3378761/1‑93 (MQ=255)

                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACgtcgtc                                                   <  1:4051326/93‑1 (MQ=255)

                                       aTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCa                                                  <  1:7432285/93‑1 (MQ=255)

                                       aTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCa                                                  >  1:3010325/1‑93 (MQ=255)

                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑GCCCGGTTCAGGCACGTCGTCAc                                                 >  2:1315688/1‑93 (MQ=255)

                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCAc                                                 >  2:6377837/1‑93 (MQ=255)

                                         gCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCAct                                                >  2:6446160/1‑93 (MQ=255)

                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCActc                                               <  2:6193847/93‑1 (MQ=255)

                                           cATGAACTGGCACTGCTGCGCGATGCGCTGCACAGTGGTGACCCGGCGGCTCTGGCAGAGTGGAGCG‑GCCCGATTCAGGCACGTCGTCAcgct                                              >  2:546784/1‑93 (MQ=255)

                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTAc                                            <  2:4641834/93‑1 (MQ=255)

                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTAc                                            >  1:5104253/1‑93 (MQ=255)

                                               aaCTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTAccc                                          <  2:923856/93‑1 (MQ=255)

                                                aCTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGCGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCcg                                         >  1:3398024/1‑93 (MQ=255)

                                                 cTGGCACTGCTGCGCCATGCGCTGAACAGTGGTGTCACGGCGTCTCTGGCAGAGAAGAGCG‑CCTCGCTGCATGCACGTCGTCACTCCACCcgc                                        >  2:3083803/1‑93 (MQ=255)

                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       >  2:790586/1‑93 (MQ=255)

                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   >  2:4400177/1‑93 (MQ=255)

                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   <  1:4789881/93‑1 (MQ=255)

                                                        tgctgcGCGATGCCCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATa                                 <  1:6579131/93‑1 (MQ=255)

                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATa                                 <  1:2751587/93‑1 (MQ=255)

                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACAt                                  <  1:5555523/91‑1 (MQ=255)

                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                >  2:5345253/1‑93 (MQ=255)

                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                <  2:1826643/93‑1 (MQ=255)

                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                >  1:1196336/1‑93 (MQ=255)

                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               <  1:2518383/93‑1 (MQ=255)

                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               >  2:7391102/1‑93 (MQ=255)

                                                           tgcGCGATGCGCTGGACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATc                              >  2:6058892/1‑93 (MQ=255)

                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             <  1:5248916/93‑1 (MQ=255)

                                                             cgcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcg                            <  1:4573046/93‑1 (MQ=255)

                                                             cgcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcg                            <  2:4028886/93‑1 (MQ=255)

                                                              gcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcgg                           >  2:2032235/1‑93 (MQ=255)

                                                               cgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggc                          <  2:3521038/93‑1 (MQ=255)

                                                               cgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggc                          <  1:7062408/93‑1 (MQ=255)

                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCTGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         <  1:4166630/93‑1 (MQ=255)

                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         <  2:648345/93‑1 (MQ=255)

                                                                 aTGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcgg                        <  2:381810/93‑1 (MQ=255)

                                                                  tGCTCTGAACAGTGGTGACACGGCGGCTCTTGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGt                       <  2:5665754/93‑1 (MQ=255)

                                                                  tGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGt                       <  2:4600221/93‑1 (MQ=255)

                                                                  tGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGt                       <  2:4724407/93‑1 (MQ=255)

                                                                   gcgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGCCGTCACTCTACCCGCGTACATAATCCGGCGGTa                      >  1:5595587/1‑93 (MQ=255)

                                                                     gcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGa                    >  2:4782476/1‑93 (MQ=255)

                                                                     gcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGa                    >  1:4302938/1‑93 (MQ=255)

                                                                     gcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGa                    >  2:6880474/1‑93 (MQ=255)

                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  <  1:1117293/93‑1 (MQ=255)

                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 >  2:3640332/1‑93 (MQ=255)

                                                                         aaCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAg                <  1:7209045/93‑1 (MQ=255)

                                                                          aCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGc               <  1:3117512/93‑1 (MQ=255)

                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCg              <  2:2399283/93‑1 (MQ=255)

                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  2:4528775/93‑1 (MQ=255)

                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  2:6617947/93‑1 (MQ=255)

                                                                              tggtgTCACGGCGGCTCTGGCAGAGTGGAGCG‑CCCTGATTCAGGCACGGCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCt           <  1:1343218/93‑1 (MQ=255)

                                                                              tggtgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCt           <  1:4144418/93‑1 (MQ=255)

                                                                              tggtgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCt           <  2:894393/93‑1 (MQ=255)

                                                                              tggtgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCt           >  1:7738148/1‑93 (MQ=255)

                                                                                 tgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        <  2:3978442/93‑1 (MQ=255)

                                                                                 tgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        >  1:6076505/1‑93 (MQ=255)

                                                                                 tgACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        <  1:2190864/93‑1 (MQ=255)

                                                                                 ggACACGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        <  2:83441/92‑1 (MQ=255)

                                                                                   acacGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcgg      >  1:1002212/1‑93 (MQ=255)

                                                                                    cacGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcggc     <  1:5760588/93‑1 (MQ=255)

                                                                                     acGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGgcggcg    >  2:7052854/1‑93 (MQ=255)

                                                                                     acGGCGGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGgcggcg    <  2:2585630/93‑1 (MQ=255)

                                                                                      cggcggCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGa   >  1:5909301/1‑93 (MQ=255)

                                                                                      cggcggCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGa   >  1:6913028/1‑93 (MQ=255)

                                                                                      cggcggCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGa   <  2:6006133/93‑1 (MQ=255)

                                                                                       ggcggcTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAt  <  2:2305896/93‑1 (MQ=255)

                                                                                       ggcggcTCTGGCAGAGCGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATCATCCGGCGGTAGAAAAGCGTCTGGCGGCGAt  >  1:4941617/1‑93 (MQ=255)

TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCG‑CCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838

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