BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,890,766	T→C	N165S (AAC→AGC)	b1808 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,890,766	0	T	C	99.6%	760.3 / ‑5.0	243	N165S (AAC→AGC)	b1808	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (132/110); total (132/111)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 255 total aligned reads displayed.

ACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAG  >  NC_000913/1890676‑1890855

accaccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGCCAGTTGTCGCTg                                                                                         >  2:6656324/1‑93 (MQ=255)

accaccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTg                                                                                         >  1:2969662/1‑93 (MQ=255)

  caccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGt                                                                                       <  2:2857268/93‑1 (MQ=255)

  caccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGt                                                                                       <  2:6946173/93‑1 (MQ=255)

   accacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTg                                                                                      >  2:3491104/1‑93 (MQ=255)

   accacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTg                                                                                      >  1:7730580/1‑93 (MQ=255)

    ccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGc                                                                                     >  2:4922927/1‑93 (MQ=255)

      ccaTCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTg                                                                                   <  1:7788341/92‑1 (MQ=255)

      acaTCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTg                                                                                   <  2:6373492/93‑1 (MQ=255)

       caTCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTgg                                                                                  <  2:7390940/93‑1 (MQ=255)

       caTCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTgg                                                                                  <  1:5582116/93‑1 (MQ=255)

         tCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGctggtgctggtg                                                                                >  2:1521008/1‑93 (MQ=255)

           ggCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGAc                                                                              >  2:3182253/1‑93 (MQ=255)

            gCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGAcc                                                                             <  2:7397871/93‑1 (MQ=255)

              gTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAg                                                                           <  1:5460551/93‑1 (MQ=255)

              gTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAg                                                                           <  2:274403/93‑1 (MQ=255)

               tCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGc                                                                          <  1:2016485/93‑1 (MQ=255)

               tCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGc                                                                          >  1:2493016/1‑93 (MQ=255)

               tCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGc                                                                          >  1:4266747/1‑93 (MQ=255)

               tCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGc                                                                          >  1:3268166/1‑93 (MQ=255)

                ccATCGCTTTTTTACGTGCTTTGGCCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGCTGCTGGTGACCAGCg                                                                         >  1:4480580/1‑93 (MQ=255)

                 cATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCgg                                                                        >  2:2038508/1‑93 (MQ=255)

                  aTCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGc                                                                       >  1:4023200/1‑93 (MQ=255)

                   tCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGcc                                                                      >  2:1709983/1‑93 (MQ=255)

                     gCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCAc                                                                    >  1:2713265/1‑93 (MQ=255)

                      cTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACg                                                                   >  1:7211318/1‑93 (MQ=255)

                       ttttttACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGc                                                                  >  2:1180240/1‑93 (MQ=255)

                          tttACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCtg                                                               >  2:5035298/1‑93 (MQ=255)

                           ttACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCtgt                                                              >  1:4676745/1‑93 (MQ=255)

                              cGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGaa                                                           <  1:2334515/93‑1 (MQ=255)

                               gTGCTTTGACCACAAAGCAATCTTTATCCATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAAt                                                          <  1:4973113/93‑1 (MQ=255)

                                  cTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCtt                                                       >  2:5992412/1‑93 (MQ=255)

                                   tttGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTc                                                      <  1:3868192/93‑1 (MQ=255)

                                   tttGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTc                                                      <  1:3222074/93‑1 (MQ=255)

                                    ttGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTcc                                                     <  1:6450193/93‑1 (MQ=255)

                                        ccACAAAGCAATCTTTATTCATCGGGCAGTCGCTGCCAAGATAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCa                                                 <  2:6853652/93‑1 (MQ=255)

                                         cacaAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCAc                                                >  1:1725625/1‑93 (MQ=255)

                                          acaAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCAcc                                               >  1:5223380/1‑92 (MQ=255)

                                          acaAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACg                                               >  2:5393527/1‑93 (MQ=255)

                                          acaAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACg                                               >  2:1639832/1‑93 (MQ=255)

                                          acaAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACg                                               >  1:1370820/1‑93 (MQ=255)

                                           caAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGc                                              >  2:1865593/1‑93 (MQ=255)

                                            aaaGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCt                                             <  2:2206016/93‑1 (MQ=255)

                                             aaGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTg                                            >  1:1793591/1‑93 (MQ=255)

                                             aaGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTg                                            >  1:5937480/1‑93 (MQ=255)

                                              aGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGa                                           <  2:6408249/93‑1 (MQ=255)

                                              aGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGa                                           >  2:2385083/1‑93 (MQ=255)

                                              aGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGa                                           <  1:4065403/93‑1 (MQ=255)

                                               gCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGAc                                          <  2:2618806/93‑1 (MQ=255)

                                                 aaTCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGc                                        >  1:3651661/1‑93 (MQ=255)

                                                  aTCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCa                                       >  1:4319021/1‑93 (MQ=255)

                                                   tCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAg                                      >  1:6561491/1‑93 (MQ=255)

                                                    ttttATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAgg                                     <  1:1482073/92‑1 (MQ=255)

                                                     tttATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAggg                                    >  1:6796797/1‑92 (MQ=255)

                                                     tttATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGt                                    <  1:6730845/93‑1 (MQ=255)

                                                       tataCATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACTCCTGTGAATCTTCCGCCACGCTGACGCAGGTGc                                  <  2:4364223/93‑1 (MQ=255)

                                                       tataCATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGc                                  <  1:6074702/93‑1 (MQ=255)

                                                       tataCATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGc                                  <  2:4739479/93‑1 (MQ=255)

                                                       tataCATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGACACGCCTGTGAATCTTCCGCCACGCTGACGCGGgcgc                                  >  2:447097/1‑93 (MQ=255)

                                                          aCATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGa                               <  2:5722151/93‑1 (MQ=255)

                                                           cATCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCGGGTGCTGat                              <  2:5312211/93‑1 (MQ=255)

                                                            aTCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGata                             <  2:7601400/93‑1 (MQ=255)

                                                            aTCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGata                             >  1:5820911/1‑93 (MQ=255)

                                                             tCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGatat                            >  1:7486836/1‑93 (MQ=255)

                                                             tCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGatat                            >  1:5046702/1‑93 (MQ=255)

                                                             tCGGGCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGatat                            >  1:3993937/1‑93 (MQ=255)

                                                               gggCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCa                          >  1:3769415/1‑93 (MQ=255)

                                                               gggCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCa                          >  1:159243/1‑93 (MQ=255)

                                                                 gCAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCAcc                        <  1:816982/93‑1 (MQ=255)

                                                                  cAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCGGGTGCTGATATCACCa                       <  1:7374796/93‑1 (MQ=255)

                                                                  cAGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCa                       <  2:6295319/93‑1 (MQ=255)

                                                                   aGTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCAt                      <  1:6937019/93‑1 (MQ=255)

                                                                    gTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATc                     >  1:892483/1‑93 (MQ=255)

                                                                    gTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATc                     >  2:2239500/1‑93 (MQ=255)

                                                                    gTCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATc                     <  1:3400078/93‑1 (MQ=255)

                                                                     tCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCg                    >  2:5073062/1‑93 (MQ=255)

                                                                     tCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCg                    <  1:1257520/93‑1 (MQ=255)

                                                                     tCGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGGGCTGATATCACCATCg                    >  1:310208/1‑93 (MQ=255)

                                                                      cGCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGa                   <  1:4432999/93‑1 (MQ=255)

                                                                       gCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGAc                  <  2:4572502/93‑1 (MQ=255)

                                                                       gCTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGAc                  >  2:7908461/1‑93 (MQ=255)

                                                                        cTGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACt                 >  2:7088606/1‑93 (MQ=255)

                                                                         tGCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTg                >  1:7193660/1‑93 (MQ=255)

                                                                          gagaaGACAGTGGGCGCTGGTGCTGGGAACCAGCGGCCACGCCTGTGAATCTTCCGCACCGCTGACGCAGGCGCTGATATCACCATCGACTGt               <  1:1408230/90‑1 (MQ=255)

                                                                          gCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGt               >  1:7514617/1‑93 (MQ=255)

                                                                          gCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGt               <  1:6220794/93‑1 (MQ=255)

                                                                          gCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGt               <  2:1302788/93‑1 (MQ=255)

                                                                          gCCAAGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGt               <  1:1987274/93‑1 (MQ=255)

                                                                              aGACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGa           <  2:3711553/93‑1 (MQ=255)

                                                                               gACAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt          <  1:217508/93‑1 (MQ=255)

                                                                                aCAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAtt         >  2:5631408/1‑93 (MQ=255)

                                                                                aCAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAtt         >  2:3436404/1‑93 (MQ=255)

                                                                                 cAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTa        <  1:454511/93‑1 (MQ=255)

                                                                                 cAGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTa        <  2:6615115/93‑1 (MQ=255)

                                                                                  aGTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAg       >  1:547984/1‑93 (MQ=255)

                                                                                   gTTGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGa      <  1:119812/93‑1 (MQ=255)

                                                                                    ttGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGAc     <  2:812578/93‑1 (MQ=255)

                                                                                    ttGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGAc     >  1:353308/1‑93 (MQ=255)

                                                                                    ttGTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGAc     >  1:6680771/1‑93 (MQ=255)

                                                                                      gTCGCTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCa   <  2:3842290/93‑1 (MQ=255)

ACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAG  >  NC_000913/1890676‑1890855

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