BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,126,892	T→C	H202H (CAT→CAC)	glcC →	GlcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,126,892	0	T	C	98.9%	598.5 / ‑2.5	187	H202H (CAT→CAC)	glcC	GlcC
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (107/78); total (109/78)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 197 total aligned reads displayed.

ATTGACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCATGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGAATCT  >  NC_000913/3126802‑3126981

aTTGACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACgc                                                                                         <  2:3712760/93‑1 (MQ=255)

aTTGACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACgc                                                                                         >  1:345936/1‑93 (MQ=255)

  tGACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACgcgc                                                                                       <  1:4419212/93‑1 (MQ=255)

  tGACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACgcgc                                                                                       >  2:1302040/1‑93 (MQ=255)

   gACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCAcgcgcg                                                                                      >  2:2260249/1‑93 (MQ=255)

    aCCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCgg                                                                                     >  2:3927521/1‑93 (MQ=255)

      cGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGAt                                                                                   <  2:1088921/93‑1 (MQ=255)

       gATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATc                                                                                  <  1:2403082/93‑1 (MQ=255)

        aTCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCt                                                                                 >  1:564127/1‑93 (MQ=255)

        aTCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCt                                                                                 <  2:4488392/93‑1 (MQ=255)

        aTCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCt                                                                                 >  1:590819/1‑93 (MQ=255)

        aTCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCt                                                                                 >  2:192185/1‑93 (MQ=255)

         tCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTa                                                                                >  2:1548135/1‑93 (MQ=255)

           tGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACa                                                                              >  2:4568411/1‑93 (MQ=255)

           tGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATAGATCGCCAGCATCCGCGGATATACa                                                                              >  1:3916501/1‑93 (MQ=255)

            gATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACaa                                                                             <  1:4373993/93‑1 (MQ=255)

            gATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACaa                                                                             <  1:1480900/93‑1 (MQ=255)

             ttGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAAc                                                                            <  1:893771/92‑1 (MQ=255)

             aTGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAAc                                                                            <  1:1570951/93‑1 (MQ=255)

              tGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAAcg                                                                           >  1:3415081/1‑93 (MQ=255)

               gTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAAcgc                                                                          >  2:3404004/1‑93 (MQ=255)

               gTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCAACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAAcgc                                                                          >  2:2916613/1‑93 (MQ=255)

                 taaaTTCATTGTTTG‑CGGGGTAAATAATCCATTCCCATTACCGAAGGAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGgt                                                                       <  2:2881515/91‑1 (MQ=255)

                  tAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGgt                                                                       <  2:631742/93‑1 (MQ=255)

                  tAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGgt                                                                       <  1:2222221/93‑1 (MQ=255)

                    aTTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGgtgt                                                                     >  2:3030500/1‑93 (MQ=255)

                      tCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTg                                                                   >  2:1317743/1‑93 (MQ=255)

                       cAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGc                                                                  <  1:3117157/93‑1 (MQ=255)

                        aGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCa                                                                 >  2:1361768/1‑93 (MQ=255)

                        aGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCa                                                                 >  2:553949/1‑93 (MQ=255)

                        aGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCa                                                                 >  1:1997602/1‑93 (MQ=255)

                         gtgtTTGCCAGCGTTAATAATCTCTAACATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAg                                                                >  2:432791/1‑93 (MQ=255)

                         gtgtTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAg                                                                <  2:3553497/93‑1 (MQ=255)

                          gggttGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGc                                                               <  1:3231463/90‑1 (MQ=255)

                           gtTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCg                                                              <  2:1272206/93‑1 (MQ=255)

                           gtTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCg                                                              >  2:3194593/1‑93 (MQ=255)

                             ttGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGc                                                            >  2:1997581/1‑93 (MQ=255)

                             ttGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGc                                                            <  1:102110/93‑1 (MQ=255)

                             ttGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGc                                                            <  1:3084119/93‑1 (MQ=255)

                             ttGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGc                                                            <  1:1046305/93‑1 (MQ=255)

                              tGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCt                                                           <  2:4423525/93‑1 (MQ=255)

                                ccAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGc                                                         >  1:4602272/1‑93 (MQ=255)

                                 cAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGCGGTGCAGCGGCTGcc                                                        >  1:4208755/1‑93 (MQ=255)

                                    cGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCa                                                     <  1:4473207/93‑1 (MQ=255)

                                     gTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCAc                                                    <  1:509235/93‑1 (MQ=255)

                                     gTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCAc                                                    <  2:385752/93‑1 (MQ=255)

                                     gTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCAc                                                    <  2:4682636/93‑1 (MQ=255)

                                     gTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCAc                                                    <  1:3841159/93‑1 (MQ=255)

                                      tAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACg                                                   >  1:2933992/1‑93 (MQ=255)

                                       aaaTAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGt                                                  >  2:2690374/1‑93 (MQ=255)

                                          taatCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGc                                               >  1:4049490/1‑93 (MQ=255)

                                            atCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGccc                                             >  1:3569473/1‑93 (MQ=255)

                                             tctcTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCa                                            >  1:3933937/1‑93 (MQ=255)

                                             tctcTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCa                                            >  1:2067116/1‑93 (MQ=255)

                                             tctcTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCa                                            >  2:3493981/1‑93 (MQ=255)

                                             tctcTACCATCGAACACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCa                                            <  1:4070594/93‑1 (MQ=255)

                                              ctctACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAg                                           >  2:3310023/1‑93 (MQ=255)

                                                 tACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAgcgc                                        >  2:2090999/1‑93 (MQ=255)

                                                   ccATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCCgcgcgc                                      >  2:4165180/1‑93 (MQ=255)

                                                    caccGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCgca                                     <  1:3814984/90‑1 (MQ=255)

                                                    cATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCgca                                     <  1:430372/93‑1 (MQ=255)

                                                    cATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCgca                                     >  2:1520899/1‑93 (MQ=255)

                                                     aTCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCgcag                                    <  1:520518/93‑1 (MQ=255)

                                                        gACCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCAc                                 <  2:2811506/93‑1 (MQ=255)

                                                         aCCACAGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCAcg                                >  2:3107389/1‑93 (MQ=255)

                                                            acaGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGa                             >  1:2373559/1‑93 (MQ=255)

                                                              aGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCTGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATc                           <  1:463303/93‑1 (MQ=255)

                                                              aGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATc                           <  1:2089938/93‑1 (MQ=255)

                                                              aGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATc                           >  1:4238004/1‑93 (MQ=255)

                                                              aGCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATc                           <  2:1637727/93‑1 (MQ=255)

                                                               gCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCa                          <  2:2442030/93‑1 (MQ=255)

                                                               gCAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCa                          <  2:2265848/93‑1 (MQ=255)

                                                                cAAAAAAAGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCAt                         >  2:4133321/1‑93 (MQ=255)

                                                                  aaaaaaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCAtgt                       <  2:2775094/93‑1 (MQ=255)

                                                                  aaaaaaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCAtgt                       <  1:4196849/93‑1 (MQ=255)

                                                                  aaaaaaACAGATCGATCCCCAGCACGCGCGTATCTACAACGCGGTGTTCCAGCGGCTGCCTCACGTCGCCCAGCGCGCAGCACGGGATCAtgt                       <  2:3591505/93‑1 (MQ=255)

                                                                    aaaaGCAGATCGATCGGCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGc                     >  2:4506647/1‑93 (MQ=255)

                                                                    aaaaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGc                     <  2:480077/93‑1 (MQ=255)

                                                                     aaaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCg                    <  2:879221/93‑1 (MQ=255)

                                                                     aaaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCg                    <  1:2728930/93‑1 (MQ=255)

                                                                     aaaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCg                    >  2:3378902/1‑93 (MQ=255)

                                                                      aaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCgg                   >  2:2370684/1‑93 (MQ=255)

                                                                      aaGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCgg                   >  2:2007387/1‑93 (MQ=255)

                                                                       aGCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGa                  <  1:4681809/93‑1 (MQ=255)

                                                                        gCAGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGAc                 >  1:3037740/1‑93 (MQ=255)

                                                                          aGATCGATCGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCg               >  2:2021761/1‑93 (MQ=255)

                                                                           gatcgatcGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGt              >  1:2760798/1‑93 (MQ=255)

                                                                           gatcgatcGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGt              >  1:1004208/1‑93 (MQ=255)

                                                                           gatcgatcGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGt              <  2:216962/93‑1 (MQ=255)

                                                                           gatcgatcGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGt              <  1:3881348/93‑1 (MQ=255)

                                                                           gatcgatcGCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGt              >  1:333126/1‑93 (MQ=255)

                                                                            atcgatcgCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTg             >  1:2590411/1‑93 (MQ=255)

                                                                               gatcgCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGaaa          <  1:1816450/93‑1 (MQ=255)

                                                                                   gCCAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGa      <  1:1337158/93‑1 (MQ=255)

                                                                                    ccAGCACGCGCGGATCTCCAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGaa     <  2:4206717/93‑1 (MQ=255)

                                                                                    ccAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGaa     >  1:1219663/1‑93 (MQ=255)

                                                                                    ccAGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGaa     <  1:3515176/93‑1 (MQ=255)

                                                                                      aGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGAAtc   >  1:674891/1‑93 (MQ=255)

                                                                                      aGCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGAAtc   <  2:112268/93‑1 (MQ=255)

                                                                                       gCACGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGAAtct  >  1:2253972/1‑93 (MQ=255)

ATTGACCGATCTGATGTTTAATTCAGTGTTTGCCAGCGTAAATAATCTCTACCATCGACCACAGCAAAAAAAGCAGATCGATCGCCAGCATGCGCGGATCTACAACGCGGTGTTGCAGCGGCTGCCGCACGTCGCCCAGCGCGCAGCACGCGATCATGTGCGGACCGTGAAAAAGAATCT  >  NC_000913/3126802‑3126981

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