BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,029,120	+C	coding (73/1191 nt)	b0967 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,029,116	1	.	C	99.2%	357.0 / ‑3.7	117	coding (77/1191 nt)	b0967	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (63/53); total (63/54)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 128 total aligned reads displayed.

CGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑C‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029025‑1029203

cGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c                                                                                         <  1:3367154/93‑1 (MQ=255)

 gAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c‑c                                                                                       >  1:1397713/1‑93 (MQ=255)

  aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c‑cc                                                                                      <  2:3615090/93‑1 (MQ=255)

      aGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c‑ccccgg                                                                                  >  1:4317235/1‑90 (MQ=255)

       gcgcCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c‑ccccgga                                                                                 >  1:800571/1‑89 (MQ=255)

       gcgcCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c‑ccccgga                                                                                 >  2:2508758/1‑89 (MQ=255)

       gcgcCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c‑ccccgga                                                                                 >  2:3414077/1‑89 (MQ=255)

         gccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGGAACGGAC‑ACCGGaaa                                                                                >  1:870142/1‑93 (MQ=255)

         gccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑C‑CCCCGgaaa                                                                               <  2:2897713/93‑1 (MQ=255)

           cgcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAg                                                                              >  1:3359257/1‑93 (MQ=255)

            gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGa                                                                             >  1:1048120/1‑93 (MQ=255)

            gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGa                                                                             <  1:3576250/93‑1 (MQ=255)

            gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGa                                                                             <  1:4663880/93‑1 (MQ=255)

            gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGa                                                                             >  1:2613466/1‑93 (MQ=255)

              gTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGAcc                                                                           >  1:1125837/1‑93 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGAccc                                                                          >  2:2746421/1‑93 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGAccc                                                                          >  1:1416935/1‑93 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGAccc                                                                          >  2:3973042/1‑93 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGAccc                                                                          >  1:168657/1‑93 (MQ=255)

               tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGAccc                                                                          >  2:1564173/1‑93 (MQ=255)

                gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCa                                                                         >  1:1835667/1‑93 (MQ=255)

                gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCa                                                                         <  1:3925491/93‑1 (MQ=255)

                  tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACg                                                                       >  2:1934530/1‑93 (MQ=255)

                  tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACg                                                                       >  2:4528406/1‑93 (MQ=255)

                   aaCCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACgg                                                                      >  1:399852/1‑93 (MQ=255)

                     ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGCCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGAt                                                                    <  2:3058133/93‑1 (MQ=255)

                      cATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATg                                                                   >  1:3147872/1‑93 (MQ=255)

                        ttttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATggg                                                                 >  2:78760/1‑92 (MQ=255)

                        ttttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGc                                                                 <  2:1694764/93‑1 (MQ=255)

                         tttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGcg                                                                >  2:1500116/1‑93 (MQ=255)

                          ttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGcga                                                               <  1:1243058/93‑1 (MQ=255)

                           tCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGcgac                                                              <  1:3960115/93‑1 (MQ=255)

                             ccTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGcgacga                                                            >  2:138820/1‑93 (MQ=255)

                             ccTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGcgacga                                                            >  2:980779/1‑93 (MQ=255)

                              cTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGaa                                                           <  2:960174/93‑1 (MQ=255)

                                gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGt                                                         <  2:2127251/93‑1 (MQ=255)

                                gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGt                                                         >  1:40023/1‑93 (MQ=255)

                                 atgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTa                                                        >  2:1671278/1‑93 (MQ=255)

                                  tgatCAACAATATCGATTGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTaa                                                       <  2:3176284/93‑1 (MQ=255)

                                  tgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTaa                                                       >  2:4041413/1‑93 (MQ=255)

                                     tCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                    <  1:1093449/93‑1 (MQ=255)

                                     tCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                    <  1:2444785/93‑1 (MQ=255)

                                        acaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttt                                                 <  2:4111365/93‑1 (MQ=255)

                                         caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtttt                                                >  1:1967389/1‑93 (MQ=255)

                                         caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtttt                                                <  1:1388370/93‑1 (MQ=255)

                                         caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtttt                                                <  2:1352494/93‑1 (MQ=255)

                                          aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttttt                                               <  2:3618618/93‑1 (MQ=255)

                                          aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttttt                                               >  1:2937587/1‑93 (MQ=255)

                                           atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                              <  1:213731/93‑1 (MQ=255)

                                           atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGAAATGATTTTTc                                              >  1:3229244/1‑93 (MQ=255)

                                            tatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGCAGTAATGATTTTTcg                                             >  1:4608346/1‑93 (MQ=255)

                                            tatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcg                                             >  2:2247687/1‑93 (MQ=255)

                                            tatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcg                                             >  1:4108442/1‑93 (MQ=255)

                                              tCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcg                                           >  1:2912044/1‑93 (MQ=255)

                                              tCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcg                                           >  2:795500/1‑93 (MQ=255)

                                                 aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        >  2:3182301/1‑93 (MQ=255)

                                                     tttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTg                                    <  1:301893/93‑1 (MQ=255)

                                                      ttCACCGAGGCTGGCTTTGCCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTgg                                   <  1:3371056/93‑1 (MQ=255)

                                                       tCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGc                                  >  2:4561561/1‑93 (MQ=255)

                                                         aCCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                <  2:2226834/93‑1 (MQ=255)

                                                         aCCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                >  1:2158085/1‑93 (MQ=255)

                                                          ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               >  1:2133187/1‑93 (MQ=255)

                                                          ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               >  2:2620315/1‑93 (MQ=255)

                                                          ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               >  2:412101/1‑93 (MQ=255)

                                                          ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               <  2:2013185/93‑1 (MQ=255)

                                                          ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               <  2:3449044/93‑1 (MQ=255)

                                                           cGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAac                              <  2:575557/93‑1 (MQ=255)

                                                            gAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                             >  2:1781260/1‑93 (MQ=255)

                                                            gAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                             >  1:2109331/1‑93 (MQ=255)

                                                            gAGGCTGGCTTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                             >  2:3264634/1‑93 (MQ=255)

                                                              ggctggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                           >  2:2989904/1‑93 (MQ=255)

                                                              ggctggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                           >  1:1624493/1‑93 (MQ=255)

                                                              ggctggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAAAACt                           >  1:298427/1‑93 (MQ=255)

                                                               gctggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTa                          >  2:3450604/1‑93 (MQ=255)

                                                                 tggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaaa                        >  1:4322657/1‑93 (MQ=255)

                                                                 tggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaaa                        <  1:2022065/93‑1 (MQ=255)

                                                                 cggctTTGCCTTCCATGCGGACAACGG‑CTCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCTCTCTTTGGCTAACTCTaaa                        <  2:698532/92‑1 (MQ=255)

                                                                  ggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                       <  2:792709/93‑1 (MQ=255)

                                                                  ggctTTACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                       >  1:902399/1‑93 (MQ=255)

                                                                      ttACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                   <  2:179899/93‑1 (MQ=255)

                                                                      ttACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                   >  1:4011996/1‑93 (MQ=255)

                                                                      ttACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                   <  1:323067/93‑1 (MQ=255)

                                                                       tACCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTaca                  <  2:2926208/93‑1 (MQ=255)

                                                                        aCCTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTacac                 <  2:168844/93‑1 (MQ=255)

                                                                         ccTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACt                >  1:3658699/1‑93 (MQ=255)

                                                                          cTTCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTc               <  1:1367423/93‑1 (MQ=255)

                                                                           ttCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCa              >  2:2618809/1‑93 (MQ=255)

                                                                           ttCCATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCa              <  1:4335572/93‑1 (MQ=255)

                                                                              cATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           <  2:4286947/93‑1 (MQ=255)

                                                                              cATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           >  1:4669538/1‑93 (MQ=255)

                                                                              cATGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           <  1:2520614/93‑1 (MQ=255)

                                                                               aTGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt          <  1:3431483/93‑1 (MQ=255)

                                                                               aTGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt          <  1:358280/93‑1 (MQ=255)

                                                                                tGCGGGCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAtt         >  1:985826/1‑93 (MQ=255)

                                                                                   gggCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTAc      <  1:2025386/93‑1 (MQ=255)

                                                                                    ggCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt     >  1:3663627/1‑93 (MQ=255)

                                                                                     gCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtt    <  2:3028413/93‑1 (MQ=255)

                                                                                     gCAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtt    <  1:2442515/93‑1 (MQ=255)

                                                                                      cAACGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACttt   <  2:91239/93‑1 (MQ=255)

                                                                                       aaCGGCC‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  <  1:358303/93‑1 (MQ=255)

CGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑C‑CCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029025‑1029203

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