BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	159,762	A→G	R117R (CGT→CGC)	yadB ←	putative tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	159,762	0	A	G	88.8%	807.7 / 73.9	305	R117R (CGT→CGC)	yadB	putative tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (17/17); new base G (126/145); total (143/162)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 315 total aligned reads displayed.

GGATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCA  >  NC_000913/159670‑159854

ggATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATa                                                                                              <  1:3861988/93‑1 (MQ=255)

   tgctgcTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATgcgc                                                                                           >  1:252362/1‑93 (MQ=255)

   tgctgcTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATAcgc                                                                                           >  2:3773005/1‑93 (MQ=255)

     ctgctgGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATgcgcgc                                                                                         >  1:661175/1‑93 (MQ=255)

      tgctgGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCa                                                                                        <  2:4049610/93‑1 (MQ=255)

      tgctgGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCa                                                                                        >  2:6436601/1‑93 (MQ=255)

         tgGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCAcgc                                                                                     >  1:4100705/1‑93 (MQ=255)

          ggCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCAcgcg                                                                                    >  2:2582136/1‑93 (MQ=255)

            cGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTa                                                                                  <  1:5905162/93‑1 (MQ=255)

            cGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTa                                                                                  <  2:5414290/93‑1 (MQ=255)

             ggATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTAc                                                                                 >  2:6217970/1‑93 (MQ=255)

             ggATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTAc                                                                                 <  2:7450174/93‑1 (MQ=255)

              gATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACa                                                                                <  2:1640501/93‑1 (MQ=255)

              gATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACa                                                                                <  2:7206847/93‑1 (MQ=255)

                 aCGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGt                                                                             <  1:4270120/93‑1 (MQ=255)

                  cGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTg                                                                            <  2:6697939/93‑1 (MQ=255)

                  cGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTg                                                                            <  2:2933282/93‑1 (MQ=255)

                      cTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGt                                                                        <  1:1932438/93‑1 (MQ=255)

                       tGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGta                                                                       >  1:4958357/1‑93 (MQ=255)

                         cggcgTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGtaat                                                                     <  1:97444/93‑1 (MQ=255)

                         cggcgTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGtaat                                                                     <  1:5784636/93‑1 (MQ=255)

                           gcgTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGtaataa                                                                   >  1:3222949/1‑93 (MQ=255)

                             gTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACt                                                                 >  1:1684143/1‑93 (MQ=255)

                              ttGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACtt                                                                <  2:4346643/93‑1 (MQ=255)

                              ttGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACtt                                                                >  2:5592382/1‑93 (MQ=255)

                               tGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTa                                                               <  2:457610/93‑1 (MQ=255)

                                 tCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGt                                                             <  2:2312118/93‑1 (MQ=255)

                                 tCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGt                                                             <  1:2052392/93‑1 (MQ=255)

                                    ggTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCt                                                          >  1:3581991/1‑93 (MQ=255)

                                    ggTCCATGATACAACACCCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTcc                                                           <  1:6153490/93‑1 (MQ=255)

                                     gTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCtt                                                         <  1:5741377/93‑1 (MQ=255)

                                     gTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCtt                                                         <  1:5526433/93‑1 (MQ=255)

                                     gTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCtt                                                         <  1:1424503/93‑1 (MQ=255)

                                       ccATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGt                                                       <  2:6782011/93‑1 (MQ=255)

                                          tgatgCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCa                                                    <  1:337406/93‑1 (MQ=255)

                                           gatgCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCAt                                                   <  2:5991195/93‑1 (MQ=255)

                                             tgCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGt                                                 >  2:440406/1‑93 (MQ=255)

                                             tgCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGt                                                 <  2:1570228/93‑1 (MQ=255)

                                                aaCA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAAc                                              >  1:5155587/1‑93 (MQ=255)

                                                aaCA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAAc                                              <  2:2990311/93‑1 (MQ=255)

                                                aaCA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAAc                                              >  1:7117536/1‑93 (MQ=255)

                                                aaCA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAAc                                              >  1:6790442/1‑93 (MQ=255)

                                                 aca‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAAcc                                             <  1:2427242/93‑1 (MQ=255)

                                                  ca‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCa                                            >  1:1790222/1‑93 (MQ=255)

                                                  ca‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCa                                            <  1:156078/93‑1 (MQ=255)

                                                   a‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAg                                           <  2:4283529/93‑1 (MQ=255)

                                                   a‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAg                                           >  2:5889040/1‑93 (MQ=255)

                                                   a‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAg                                           >  1:6775979/1‑93 (MQ=255)

                                                   a‑cCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAg                                           <  1:5226176/93‑1 (MQ=255)

                                                      ccGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGc                                         <  1:920749/93‑1 (MQ=255)

                                                        ggCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGa                                       <  1:2087417/93‑1 (MQ=255)

                                                         gCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAg                                      >  2:2779996/1‑93 (MQ=255)

                                                         gCAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAg                                      >  2:2693456/1‑93 (MQ=255)

                                                          cAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGt                                     <  1:5374397/93‑1 (MQ=255)

                                                          cAATGACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGt                                     >  2:1388700/1‑93 (MQ=255)

                                                              gACCGTCGTAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCtt                                 <  1:2742066/93‑1 (MQ=255)

                                                                 cgtcgtAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCAc                              <  2:3919011/93‑1 (MQ=255)

                                                                  gtcgtAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACg                             >  2:4263466/1‑93 (MQ=255)

                                                                   tcgtAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGa                            >  1:2644391/1‑93 (MQ=255)

                                                                    cgtAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGAt                           >  1:4404813/1‑93 (MQ=255)

                                                                     gtAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATa                          <  2:6539707/93‑1 (MQ=255)

                                                                      tAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAg                         >  2:2122361/1‑93 (MQ=255)

                                                                      tAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAg                         >  1:6944672/1‑93 (MQ=255)

                                                                      tAAATACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAg                         <  2:706292/93‑1 (MQ=255)

                                                                       aaaTACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAgg                        <  2:6024924/93‑1 (MQ=255)

                                                                        aaTACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGc                       >  1:1459007/1‑93 (MQ=255)

                                                                        aaTACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGc                       >  1:1000833/1‑93 (MQ=255)

                                                                        aaTACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGc                       <  1:3764398/93‑1 (MQ=255)

                                                                         aTACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcg                      <  2:3727757/93‑1 (MQ=255)

                                                                          tACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgt                     >  2:1953469/1‑93 (MQ=255)

                                                                          tACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgt                     >  2:5322379/1‑93 (MQ=255)

                                                                          tACCGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgt                     >  2:4173326/1‑93 (MQ=255)

                                                                            ccGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgtcg                   >  1:2242293/1‑93 (MQ=255)

                                                                            ccGCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgtcg                   <  1:1228165/93‑1 (MQ=255)

                                                                              gCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTg                 >  1:6086088/1‑93 (MQ=255)

                                                                              gCCAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTg                 >  1:6706366/1‑93 (MQ=255)

                                                                              gCCAATGCTTTGAATGCGCGCACGCGCACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTg                 >  1:800739/1‑93 (MQ=255)

                                                                               ccAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGa                <  2:6909482/93‑1 (MQ=255)

                                                                               ccAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGa                <  1:5082633/93‑1 (MQ=255)

                                                                                cAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGAc               >  2:4946925/1‑93 (MQ=255)

                                                                                cAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGAc               >  1:3114133/1‑93 (MQ=255)

                                                                                cAATGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGAc               >  1:371106/1‑93 (MQ=255)

                                                                                 aaTGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACg              >  1:6511605/1‑93 (MQ=255)

                                                                                  aTGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGt             <  1:289447/93‑1 (MQ=255)

                                                                                   tGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGtt            >  1:2830244/1‑93 (MQ=255)

                                                                                   tGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGtt            <  2:4912192/93‑1 (MQ=255)

                                                                                   tGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGtt            >  1:3674465/1‑93 (MQ=255)

                                                                                   tGCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGtt            <  1:4610816/93‑1 (MQ=255)

                                                                                    gCTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTg           >  2:1868886/1‑93 (MQ=255)

                                                                                     cTTTGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGc          <  2:2074862/93‑1 (MQ=255)

                                                                                       ttGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGa        <  2:814973/93‑1 (MQ=255)

                                                                                       ttGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGa        <  1:7358751/93‑1 (MQ=255)

                                                                                        tGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGAt       <  2:5796182/93‑1 (MQ=255)

                                                                                        tGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGAt       <  2:5052733/93‑1 (MQ=255)

                                                                                        tGAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGAt       <  2:3295456/93‑1 (MQ=255)

                                                                                        tGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGAt       <  2:147173/93‑1 (MQ=255)

                                                                                         gAATGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGAtt      <  1:234691/93‑1 (MQ=255)

                                                                                          aaTGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTg     >  1:891587/1‑93 (MQ=255)

                                                                                          aaTGCGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTg     <  1:3013881/93‑1 (MQ=255)

                                                                                            tACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGcc   >  2:3472947/1‑93 (MQ=255)

GGATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACA‑CCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCA  >  NC_000913/159670‑159854

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