BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	800,958	A→G	Q326R (CAG→CGG)	ybhH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	800,958	0	A	G	99.5%	635.5 / NA	196	Q326R (CAG→CGG)	ybhH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (86/109); total (87/109)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

CCCGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCAGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTGT  >  NC_000913/800868‑801047

  cGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGAt                                                                                       <  2:3438488/93‑1 (MQ=255)

  cGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGAt                                                                                       <  1:4533152/93‑1 (MQ=255)

    aCAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGc                                                                                     >  2:6116427/1‑93 (MQ=255)

      aaaTCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCa                                                                                   <  2:4873663/93‑1 (MQ=255)

      aaaTCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCa                                                                                   >  2:5206879/1‑93 (MQ=255)

       aaTCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCac                                                                                  >  2:4907228/1‑93 (MQ=255)

           gTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCacgacg                                                                              >  2:1536180/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             <  1:5554847/93‑1 (MQ=255)

             cccTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGtt                                                                            >  1:4283166/1‑93 (MQ=255)

              ccTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTTGTGCGCTCGACGTTCATTTAAGTAATGAAGGTTGGGATGTCACGACGTTa                                                                           >  1:5003450/1‑93 (MQ=255)

              ccTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTa                                                                           <  1:208073/93‑1 (MQ=255)

               cTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAc                                                                          <  2:5607074/93‑1 (MQ=255)

               cTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAc                                                                          >  1:502636/1‑93 (MQ=255)

               cTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAc                                                                          >  1:3121388/1‑93 (MQ=255)

                ttCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAcg                                                                         <  2:4431150/93‑1 (MQ=255)

                  cTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAcgcg                                                                       <  2:3808747/93‑1 (MQ=255)

                    gTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTACGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCa                                                                     >  2:5352379/1‑93 (MQ=255)

                      aGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATc                                                                   <  1:4198542/93‑1 (MQ=255)

                       ggATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCt                                                                  >  1:4992318/1‑93 (MQ=255)

                       ggATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCt                                                                  >  2:160228/1‑93 (MQ=255)

                        gATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTg                                                                 >  1:5685624/1‑93 (MQ=255)

                        gATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGNTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTg                                                                 >  2:4815717/1‑93 (MQ=255)

                         aTACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGt                                                                <  1:5455068/93‑1 (MQ=255)

                          tACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGtt                                                               <  1:1679180/93‑1 (MQ=255)

                           aCGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTa                                                              <  2:5561408/93‑1 (MQ=255)

                             ggCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTAtt                                                            <  1:1925263/93‑1 (MQ=255)

                              gCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTc                                                           <  2:5884530/93‑1 (MQ=255)

                                aatatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCgg                                                         >  2:4624798/1‑93 (MQ=255)

                                aatatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCgg                                                         <  2:4700457/93‑1 (MQ=255)

                                 atatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGga                                                        >  2:4197643/1‑93 (MQ=255)

                                  tatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgac                                                       >  2:4947650/1‑93 (MQ=255)

                                   atcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgacg                                                      <  1:3212727/93‑1 (MQ=255)

                                    tcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgacga                                                     <  1:1277083/93‑1 (MQ=255)

                                     caatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgacgac                                                    <  2:2801102/93‑1 (MQ=255)

                                     caatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgacgac                                                    <  1:418918/93‑1 (MQ=255)

                                      aatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGAcc                                                   >  1:4690704/1‑93 (MQ=255)

                                      aatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGAcc                                                   <  2:1342988/93‑1 (MQ=255)

                                       atatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCa                                                  >  2:2355075/1‑93 (MQ=255)

                                       atatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCa                                                  <  1:3884515/93‑1 (MQ=255)

                                       atatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCa                                                  >  1:3761173/1‑93 (MQ=255)

                                       atatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCa                                                  <  2:1801679/93‑1 (MQ=255)

                                         atTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGa                                                <  1:973871/93‑1 (MQ=255)

                                             aaCACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaaaaa                                            <  1:5070236/93‑1 (MQ=255)

                                             aaCACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaaaaa                                            <  2:1402380/93‑1 (MQ=255)

                                             aaCACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGAGCATCTGTTATTCGGACGACCAGaaaaa                                            >  2:6640425/1‑93 (MQ=255)

                                               ccccccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAata                                          <  1:4779126/91‑1 (MQ=255)

                                                acCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAatat                                         <  1:867534/93‑1 (MQ=255)

                                                acCCCAGTGGTGCGCGCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGGTATGCGGACGACCAGAAAAatat                                         <  1:780702/93‑1 (MQ=255)

                                                  cccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAttt                                       <  1:3502839/93‑1 (MQ=255)

                                                  cccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAttt                                       <  1:4372451/93‑1 (MQ=255)

                                                  cccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACAAGAAAAATAttt                                       >  2:6009642/1‑93 (MQ=255)

                                                   ccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAtttt                                      <  2:722617/93‑1 (MQ=255)

                                                   ccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAtttt                                      <  2:5150852/93‑1 (MQ=255)

                                                     aGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTcc                                    <  1:4908220/93‑1 (MQ=255)

                                                      gtggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCg                                   <  2:376795/93‑1 (MQ=255)

                                                      gtggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCg                                   >  2:6603368/1‑93 (MQ=255)

                                                       tggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCgg                                  >  1:5360497/1‑93 (MQ=255)

                                                        ggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGt                                 >  1:5914811/1‑93 (MQ=255)

                                                        ggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGt                                 <  2:4062397/93‑1 (MQ=255)

                                                           gcgcTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGaa                              >  1:3355778/1‑93 (MQ=255)

                                                            cgcTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAg                             >  2:5515525/1‑93 (MQ=255)

                                                            cgcTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAg                             <  1:4434491/93‑1 (MQ=255)

                                                            cgcTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAg                             >  1:3673839/1‑93 (MQ=255)

                                                             gcTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGt                            >  1:5188346/1‑93 (MQ=255)

                                                             gcTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGt                            >  1:3007975/1‑93 (MQ=255)

                                                              cTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGtt                           <  1:221270/93‑1 (MQ=255)

                                                              cTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGtt                           <  1:4698978/93‑1 (MQ=255)

                                                              cTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGtt                           <  2:3697990/93‑1 (MQ=255)

                                                               tCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGttt                          >  1:2865771/1‑93 (MQ=255)

                                                               tCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGttt                          >  2:4028324/1‑93 (MQ=255)

                                                               tCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGttt                          <  2:2199797/93‑1 (MQ=255)

                                                                cGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTa                         <  2:1011094/93‑1 (MQ=255)

                                                                cGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTa                         <  2:4755825/93‑1 (MQ=255)

                                                                cGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTCCCGGTGAAGTTTa                         <  1:6437191/93‑1 (MQ=255)

                                                                  aCGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATc                       >  2:3374170/1‑93 (MQ=255)

                                                                  aCGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATc                       <  2:5845984/93‑1 (MQ=255)

                                                                  aCGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATc                       >  2:4147489/1‑93 (MQ=255)

                                                                  aCGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATc                       <  1:5976155/93‑1 (MQ=255)

                                                                  aCGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATc                       <  1:6484312/93‑1 (MQ=255)

                                                                     ttCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTc                    <  2:772868/93‑1 (MQ=255)

                                                                     ttCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTc                    <  1:1443852/93‑1 (MQ=255)

                                                                     ttCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTc                    >  2:3130200/1‑93 (MQ=255)

                                                                     ttCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTc                    >  2:2978993/1‑93 (MQ=255)

                                                                     ttCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTc                    >  1:1077533/1‑93 (MQ=255)

                                                                      tCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTcc                   <  1:3169104/93‑1 (MQ=255)

                                                                       cATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTcca                  <  2:1646877/93‑2 (MQ=255)

                                                                         tttAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTg                >  1:2492029/1‑93 (MQ=255)

                                                                            aaGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaa             >  2:470765/1‑93 (MQ=255)

                                                                               tAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaaa          <  1:6638998/93‑1 (MQ=255)

                                                                               tAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaaa          <  1:6721619/93‑1 (MQ=255)

                                                                               tAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaaa          <  1:1627852/93‑1 (MQ=255)

                                                                                aaTGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAAt         <  2:2868737/93‑1 (MQ=255)

                                                                                aaTGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAAt         <  1:1807441/93‑1 (MQ=255)

                                                                                 aTGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAAtt        <  1:6430716/93‑1 (MQ=255)

                                                                                    aaGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGt     >  1:4046345/1‑93 (MQ=255)

                                                                                     aGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGtt    >  2:647418/1‑93 (MQ=255)

                                                                                     aGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGtt    <  1:2707310/93‑1 (MQ=255)

                                                                                      ggTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTg   <  2:2545976/93‑1 (MQ=255)

                                                                                      ggTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTg   >  2:3231241/1‑93 (MQ=255)

                                                                                       gTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTGt  >  1:6199264/1‑93 (MQ=255)

CCCGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCAGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTGT  >  NC_000913/800868‑801047

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