BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,751,176	G→A	L227F (CTT→TTT)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,751,176	0	G	A	99.5%	758.7 / ‑5.4	204	L227F (CTT→TTT)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (109/94); total (110/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 208 total aligned reads displayed.

TACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCTTTCA  >  NC_000913/1751086‑1751265

tACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAAtt                                                                                         >  2:3929894/1‑93 (MQ=255)

tACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAAtt                                                                                         >  2:5065783/1‑93 (MQ=255)

tACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAAtt                                                                                         >  1:300619/1‑93 (MQ=255)

 aCGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTc                                                                                        <  2:2657778/93‑1 (MQ=255)

 aCGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTc                                                                                        >  1:3884378/1‑93 (MQ=255)

   gTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAg                                                                                      <  1:5352426/93‑1 (MQ=255)

   gTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAg                                                                                      >  1:3310406/1‑93 (MQ=255)

     gCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGTCATGGTTGTTATTCGCACCAATAAATTCAGtc                                                                                    >  1:675868/1‑93 (MQ=255)

     gCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtc                                                                                    <  1:1201753/93‑1 (MQ=255)

     gCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtc                                                                                    <  2:238108/93‑1 (MQ=255)

        gTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatc                                                                                 <  1:5145847/93‑1 (MQ=255)

        gTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGtcatc                                                                                 <  2:5902290/93‑1 (MQ=255)

            aGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGt                                                                             >  1:4893189/1‑93 (MQ=255)

             gCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTa                                                                            <  2:4179484/93‑1 (MQ=255)

             gCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTa                                                                            >  2:5349193/1‑93 (MQ=255)

             gCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTAATTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTa                                                                            <  1:3610802/93‑1 (MQ=255)

               gTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAAt                                                                          >  2:1937729/1‑93 (MQ=255)

               gTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAAt                                                                          <  2:3534604/93‑1 (MQ=255)

               gTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAAt                                                                          <  2:692669/93‑1 (MQ=255)

               gTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAAt                                                                          <  1:409546/93‑1 (MQ=255)

                tGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATc                                                                         <  1:5117116/93‑1 (MQ=255)

                  aCTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCAt                                                                       <  2:4067406/93‑1 (MQ=255)

                   cTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCAtt                                                                      >  1:6288842/1‑93 (MQ=255)

                   cTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCAtt                                                                      >  2:3390816/1‑93 (MQ=255)

                   cTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCAtt                                                                      >  2:4691698/1‑93 (MQ=255)

                   cTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCAtt                                                                      <  1:1898576/93‑1 (MQ=255)

                    ttGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTTGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTg                                                                     >  2:2523810/1‑93 (MQ=255)

                    ttGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTg                                                                     <  2:6001209/93‑1 (MQ=255)

                    ttGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTg                                                                     >  2:1758038/1‑93 (MQ=255)

                     tGGTATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTga                                                                    <  2:5514988/93‑1 (MQ=255)

                     tGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTga                                                                    <  1:3570365/93‑1 (MQ=255)

                     tGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTga                                                                    >  2:6182898/1‑93 (MQ=255)

                        gATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGAc                                                                 <  1:5997211/93‑1 (MQ=255)

                        gATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGAc                                                                 <  1:2644181/93‑1 (MQ=255)

                         aTCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACg                                                                >  2:2674226/1‑93 (MQ=255)

                          tCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGc                                                               >  2:6281326/1‑93 (MQ=255)

                            tGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCtt                                                             >  1:2037217/1‑93 (MQ=255)

                            tGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGAtt                                                             >  1:6027886/1‑93 (MQ=255)

                             gAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTc                                                            <  2:3432236/93‑1 (MQ=255)

                              aaTACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCa                                                           >  2:5130358/1‑93 (MQ=255)

                                tACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCAtc                                                         >  1:4755513/1‑93 (MQ=255)

                                tACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCAtc                                                         >  2:445538/1‑93 (MQ=255)

                                   tctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTcc                                                      >  1:1854375/1‑93 (MQ=255)

                                   tctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTcc                                                      <  1:4246967/93‑1 (MQ=255)

                                   tctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTcc                                                      <  1:21954/93‑1 (MQ=255)

                                    ctGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCt                                                     <  2:3857898/93‑1 (MQ=255)

                                     tGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTg                                                    >  1:3104224/1‑93 (MQ=255)

                                       cccACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGAc                                                  <  2:3293656/93‑1 (MQ=255)

                                          aCGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGAt                                               <  1:1001842/93‑1 (MQ=255)

                                           cGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATc                                              <  1:4579273/93‑1 (MQ=255)

                                            gATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCg                                             >  1:1453522/1‑93 (MQ=255)

                                            gATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCg                                             >  2:1436367/1‑93 (MQ=255)

                                            gATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCg                                             >  2:5750700/1‑93 (MQ=255)

                                             aTTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCgg                                            >  1:5325453/1‑93 (MQ=255)

                                             aTTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCgg                                            >  2:2326798/1‑93 (MQ=255)

                                                gtgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCaa                                         >  2:3310674/1‑93 (MQ=255)

                                                 tgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAAt                                        <  1:1109190/93‑1 (MQ=255)

                                                 tgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAAt                                        <  1:6112900/93‑1 (MQ=255)

                                                 tgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAAt                                        >  1:4632967/1‑93 (MQ=255)

                                                 tgGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAAt                                        <  2:6239546/93‑1 (MQ=255)

                                                  ggCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATg                                       >  2:2146970/1‑93 (MQ=255)

                                                   gCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGt                                      <  1:826978/93‑1 (MQ=255)

                                                   gCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGt                                      >  2:6079582/1‑93 (MQ=255)

                                                   gCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGt                                      >  1:3509972/1‑93 (MQ=255)

                                                   gCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGt                                      >  1:4500837/1‑93 (MQ=255)

                                                   gCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGt                                      <  1:934687/93‑1 (MQ=255)

                                                     gTACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTg                                    <  2:4509803/93‑1 (MQ=255)

                                                      tACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGa                                   >  2:1154955/1‑93 (MQ=255)

                                                      tACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGa                                   >  1:607877/1‑93 (MQ=255)

                                                      tACTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGa                                   >  1:1930945/1‑93 (MQ=255)

                                                       aCTTGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAc                                  >  2:2119091/1‑93 (MQ=255)

                                                         ttGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAccc                                <  1:1811228/93‑1 (MQ=255)

                                                         ttGGCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGAccc                                >  1:971913/1‑93 (MQ=255)

                                                           ggCACGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCt                              <  1:561851/93‑1 (MQ=255)

                                                              aCGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                           >  2:2980028/1‑93 (MQ=255)

                                                              aCGACATGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTc                           >  1:3680525/1‑93 (MQ=255)

                                                                   aTGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTccc                      >  2:578878/1‑93 (MQ=255)

                                                                    tGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCt                     >  1:2754581/1‑93 (MQ=255)

                                                                    tGGTTGTTATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCt                     >  1:5478948/1‑93 (MQ=255)

                                                                       ttgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCa                  <  1:5021381/93‑1 (MQ=255)

                                                                       ttgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCa                  >  1:3743074/1‑93 (MQ=255)

                                                                       ttgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCa                  <  2:5732083/93‑1 (MQ=255)

                                                                        tgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCaa                 <  2:1367240/93‑1 (MQ=255)

                                                                        tgttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCaa                 <  2:562140/93‑1 (MQ=255)

                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTTAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc                <  2:2049551/93‑1 (MQ=255)

                                                                         gttATTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAAc                <  1:2454408/93‑1 (MQ=255)

                                                                            aTTCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGc             >  2:2238439/1‑93 (MQ=255)

                                                                             ttCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCa            >  2:5171685/1‑93 (MQ=255)

                                                                              tCGCACCAATAAATTCAGTCATCATGTAATCATTGATACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCaa           <  1:3115781/93‑1 (MQ=255)

                                                                              tCGCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCaa           <  2:2613074/93‑1 (MQ=255)

                                                                                gCACCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATc         <  2:4788240/93‑1 (MQ=255)

                                                                                  aCCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGc       >  2:1280764/1‑93 (MQ=255)

                                                                                  aCCAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGc       >  1:1493334/1‑93 (MQ=255)

                                                                                   ccAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCt      >  2:5859844/1‑93 (MQ=255)

                                                                                   ccAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCt      >  2:706308/1‑93 (MQ=255)

                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt     >  1:2871566/1‑93 (MQ=255)

                                                                                    cAATAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCtt     >  2:1573765/1‑93 (MQ=255)

                                                                                     aaTAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCttt    >  1:4213401/1‑93 (MQ=255)

                                                                                     aaTAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCttt    >  2:4332099/1‑93 (MQ=255)

                                                                                       tAAATTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCTTTCa  >  1:5747/1‑93 (MQ=255)

TACGTGCTGTCCAGCGTGACTTGGGATCTGAATACTCTGCCCACGATTGTGGCGTACTTGGCACGACATGGTTGTTATTCGCACCAATAAGTTCAGTCATCATGTAATCATTGAGACGCTTCATCTCCTGACGATCGGCAATGTTGACCCCTTTCGTCCCTTCAACCGCAATCGCTTTCA  >  NC_000913/1751086‑1751265

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