BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,991,845	T→C	N62N (AAT→AAC)	b2853 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,991,845	0	T	C	98.6%	676.5 / ‑1.3	211	N62N (AAT→AAC)	b2853	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (108/99); total (111/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 224 total aligned reads displayed.

CACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTT  >  NC_000913/2991755‑2991934

cacCAGTGCTGCTGATACAGATCATATACAAGCATTTATAGTATCAACATGGATGGCGCCTTTTCAAAGTGATATGAATTCAAAAGTTAAccc                                                                                         <  2:4227126/93‑1 (MQ=255)

cacCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAAccc                                                                                         <  2:2312441/93‑1 (MQ=255)

cacCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAAccc                                                                                         <  1:3495743/93‑1 (MQ=255)

 acCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCt                                                                                        >  1:4541129/1‑93 (MQ=255)

  ccAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTa                                                                                       <  2:3590039/93‑1 (MQ=255)

  ccAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTa                                                                                       <  2:4438763/93‑1 (MQ=255)

    aGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtc                                                                                     <  2:1923288/93‑1 (MQ=255)

      tGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtctc                                                                                   >  1:1428513/1‑93 (MQ=255)

      tGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtctc                                                                                   >  1:1842260/1‑93 (MQ=255)

         tACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATACCCCTATCTCAcc                                                                                >  1:4666248/1‑93 (MQ=255)

          aCTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCt                                                                               <  2:1704057/93‑1 (MQ=255)

             gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTAc                                                                            <  2:1988986/93‑1 (MQ=255)

             gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTAc                                                                            <  1:2883667/93‑1 (MQ=255)

              aTACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACt                                                                           >  2:4558640/1‑93 (MQ=255)

               tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACta                                                                          >  2:3983377/1‑93 (MQ=255)

               tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACta                                                                          >  2:3790950/1‑93 (MQ=255)

                  aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaa                                                                       <  2:191590/93‑1 (MQ=255)

                  aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaa                                                                       <  1:2148132/93‑1 (MQ=255)

                  aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaa                                                                       >  1:1364476/1‑93 (MQ=255)

                   gATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaaa                                                                      <  1:897921/93‑1 (MQ=255)

                     tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAt                                                                    >  2:2697620/1‑93 (MQ=255)

                     tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAt                                                                    <  2:3299808/93‑1 (MQ=255)

                       aaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAtta                                                                  >  1:689490/1‑92 (MQ=255)

                       aaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTg                                                                  >  1:1258947/1‑93 (MQ=255)

                        aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGa                                                                 >  2:4030584/1‑93 (MQ=255)

                        aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGa                                                                 <  1:1523479/93‑1 (MQ=255)

                         aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAg                                                                >  1:798231/1‑93 (MQ=255)

                          tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGt                                                               >  2:252190/1‑93 (MQ=255)

                           aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTg                                                              <  1:3218064/93‑1 (MQ=255)

                           aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTg                                                              >  1:3604277/1‑93 (MQ=255)

                            cAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTgg                                                             >  2:1660478/1‑93 (MQ=255)

                              aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTa                                                           >  1:1681579/1‑93 (MQ=255)

                               gaaTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAg                                                          >  2:2218714/3‑93 (MQ=255)

                                 aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGag                                                        >  1:200645/1‑92 (MQ=255)

                                 aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGAc                                                        <  1:2352589/93‑1 (MQ=255)

                                  tttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACg                                                       >  2:2999875/1‑93 (MQ=255)

                                   ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGa                                                      <  1:3081713/93‑1 (MQ=255)

                                   ttATAGTTTCAACATGGATGGCCCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTAT‑AAATTGAGTGGTAGACGAt                                                     <  2:1330466/93‑1 (MQ=255)

                                    tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAt                                                     >  1:395528/1‑93 (MQ=255)

                                    tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAt                                                     <  2:3418272/93‑1 (MQ=255)

                                     ataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAtt                                                    >  2:3451427/1‑93 (MQ=255)

                                      taGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTg                                                   >  1:4272119/1‑93 (MQ=255)

                                         tttCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGAcc                                                <  1:1182509/93‑1 (MQ=255)

                                            cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaa                                             >  2:1808059/1‑93 (MQ=255)

                                            cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaa                                             <  2:4474311/93‑1 (MQ=255)

                                            cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaa                                             <  1:3196173/93‑1 (MQ=255)

                                             aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaat                                            <  2:540501/93‑1 (MQ=255)

                                              aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaata                                           <  1:4635664/93‑1 (MQ=255)

                                              aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaata                                           >  1:3995224/1‑93 (MQ=255)

                                                atggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaa                                         <  2:4321034/93‑1 (MQ=255)

                                                atggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaa                                         <  2:350234/93‑1 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                        >  1:2676421/1‑93 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                        <  2:2605422/93‑1 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                        >  2:3327602/1‑93 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                        >  2:4365195/1‑93 (MQ=255)

                                                     tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATa                                    <  2:1445326/93‑1 (MQ=255)

                                                     tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATa                                    <  1:4097991/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   <  2:2053083/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   >  2:3133496/1‑93 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   <  1:2789903/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   <  1:4413678/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTACCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                   <  1:1948323/93‑1 (MQ=255)

                                                       gcgcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAAc                                  >  2:2265385/1‑93 (MQ=255)

                                                       gcgcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAGCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAAc                                  <  1:1162693/93‑1 (MQ=255)

                                                         gcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATTATAACga                                <  1:3677478/93‑1 (MQ=255)

                                                          ccTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACgag                               <  1:337801/93‑1 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             >  2:604514/1‑93 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             <  1:4553365/93‑1 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                             >  1:2054831/1‑93 (MQ=255)

                                                             tttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATa                            >  2:4270258/1‑93 (MQ=255)

                                                             tttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATa                            <  1:343274/93‑1 (MQ=255)

                                                             tttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATa                            <  2:565022/93‑1 (MQ=255)

                                                                 aaaaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                        >  1:4358583/1‑93 (MQ=255)

                                                                  aaaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtttt                       >  1:2726634/1‑93 (MQ=255)

                                                                   aaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTta                      <  2:2777311/93‑1 (MQ=255)

                                                                   aaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTta                      <  1:3345366/93‑1 (MQ=255)

                                                                   aaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTta                      <  1:1338425/93‑1 (MQ=255)

                                                                    aTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                     >  2:1086931/1‑93 (MQ=255)

                                                                     tGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                    <  1:1058024/93‑1 (MQ=255)

                                                                     tGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                    >  2:3072206/1‑93 (MQ=255)

                                                                     tGATAGGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAATGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                    >  2:4571605/1‑93 (MQ=255)

                                                                      gATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                   >  2:1502418/1‑93 (MQ=255)

                                                                        tatGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataata                 >  2:119494/1‑93 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                >  1:2954431/1‑93 (MQ=255)

                                                                          tGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt               >  2:4507600/1‑93 (MQ=255)

                                                                          tGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt               >  1:804186/1‑93 (MQ=255)

                                                                            tATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt             <  1:2430502/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            >  2:3170503/1‑93 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            <  1:2575706/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc            <  1:543047/93‑1 (MQ=255)

                                                                                cAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCg         <  2:4115158/93‑1 (MQ=255)

                                                                                 agaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt        >  2:1754757/1‑93 (MQ=255)

                                                                                 agaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt        >  1:1001773/1‑93 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt       >  1:3781409/1‑93 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt       >  2:3431865/1‑93 (MQ=255)

                                                                                    agaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc     >  2:687577/1‑93 (MQ=255)

                                                                                     gaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCt    >  2:2861177/1‑93 (MQ=255)

                                                                                      aTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt   >  1:4136558/1‑93 (MQ=255)

                                                                                       tAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt  >  1:2694417/1‑93 (MQ=255)

                                                                                       tAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt  <  2:3814689/93‑1 (MQ=255)

CACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTT  >  NC_000913/2991755‑2991934

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