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	94.4%	427.2 / 17.4	143	coding (77/1191 nt)	b0967	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (6/2); new base C (74/61); total (80/63)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑C‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCT  >  NC_000913/1029024‑1029205

  gAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCCAC‑G‑G‑c‑c                                                                                          >  1:433700/1‑93 (MQ=255)

   aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGTTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑c‑c‑c                                                                                        <  2:1764535/93‑1 (MQ=255)

   aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCACC‑G‑G‑c‑c‑c                                                                                        >  2:298651/1‑93 (MQ=255)

   aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑c‑c‑c                                                                                        >  2:1570651/1‑93 (MQ=255)

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCCCAGAGGCTGCCTTTACCTTCCATGCGGGCACA‑G‑G‑c‑c‑cc                                                                                       >  1:1642515/1‑93 (MQ=255)

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑c‑c‑cc                                                                                       >  1:2939977/1‑93 (MQ=255)

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑c‑c‑cc                                                                                       >  2:3459417/1‑93 (MQ=255)

     tAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑c‑c‑ccc                                                                                      >  2:3799655/1‑92 (MQ=255)

          gccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGCCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑C‑C‑CCCGgaaa                                                                                 >  1:1635077/1‑93 (MQ=255)

             gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGa                                                                               >  2:3132093/1‑93 (MQ=255)

               gTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGAcc                                                                             >  1:321928/1‑93 (MQ=255)

                tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGAccc                                                                            >  2:671456/1‑93 (MQ=255)

                tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGAcca                                                                            >  1:3725571/1‑92 (MQ=255)

                tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTCCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGAccc                                                                            >  2:2135592/1‑93 (MQ=255)

                tGCTAACCATTTTCCCTGATGATCAACAATATCCATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑GTC‑C‑C‑CCTGCAAAGAccc                                                                            >  1:802138/1‑93 (MQ=255)

                  cTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCAc                                                                          <  1:4089231/93‑1 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACg                                                                         >  2:1173563/1‑93 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACg                                                                         <  1:3976556/93‑1 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACg                                                                         >  1:1479917/1‑93 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACg                                                                         <  1:773061/93‑1 (MQ=255)

                     aCCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGa                                                                       >  2:2924449/1‑93 (MQ=255)

                     aCCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGa                                                                       <  2:1921022/93‑1 (MQ=255)

                      ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGAt                                                                      >  2:4431501/1‑93 (MQ=255)

                       cATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATg                                                                     >  1:331160/1‑93 (MQ=255)

                        aTTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATgg                                                                    <  2:1426656/93‑1 (MQ=255)

                          tttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGcg                                                                  >  1:2301764/1‑93 (MQ=255)

                          tttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGcg                                                                  >  1:238511/1‑93 (MQ=255)

                           ttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGcga                                                                 <  1:3152023/93‑1 (MQ=255)

                            tCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGcgac                                                                <  1:3891597/93‑1 (MQ=255)

                                tgatgaTCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGCAAAGACCCACGGATGCCGACGAAg                                                            >  1:3241367/1‑93 (MQ=255)

                                 gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGt                                                           <  1:2701354/93‑1 (MQ=255)

                                 gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGt                                                           <  2:1270293/93‑1 (MQ=255)

                                    gatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAAt                                                        <  1:4582397/93‑1 (MQ=255)

                                     atCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                       >  2:2411262/1‑93 (MQ=255)

                                      ttaacaaTATCGATGGTTTCACCGAGGCGGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                      <  1:2401914/91‑1 (MQ=255)

                                       caacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAt                                                     <  1:2021315/93‑1 (MQ=255)

                                       caacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAt                                                     >  2:2105275/1‑93 (MQ=255)

                                       caacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAt                                                     <  2:2342114/93‑1 (MQ=255)

                                        aacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtt                                                    >  1:3401221/1‑93 (MQ=255)

                                        aacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtt                                                    >  1:549865/1‑93 (MQ=255)

                                         acaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGGAAC‑G‑GCC‑C‑CCGGAAAAGACCCCCGGATGGCGACGAATTAATGAttt                                                   >  2:21860/1‑93 (MQ=255)

                                         acaaTATCGATGGTTTCACCGAGGCTGGCATTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGAACCACGGATGGCGACGAAGTAATGAttt                                                   >  2:2244422/1‑93 (MQ=255)

                                           aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAATAGACCCACGGATGGCGACGAAGTAATGAttttt                                                 >  1:4571373/1‑93 (MQ=255)

                                           aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttttt                                                 <  1:1318201/93‑1 (MQ=255)

                                           aaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑CGC‑CCCGAAAAGAGCCACGGATGGCGACGAAGTAATGAgtttt                                                 >  1:4247257/1‑93 (MQ=255)

                                           aaTATCGATGGTTTCACCGAGGCTGGCTTTACCGTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAttttt                                                 <  1:1557975/93‑1 (MQ=255)

                                            atatTGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                                <  2:419134/93‑1 (MQ=255)

                                            atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                                >  2:827367/1‑93 (MQ=255)

                                            atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                                <  2:4131530/93‑1 (MQ=255)

                                               tCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcg                                             >  1:1707380/1‑93 (MQ=255)

                                                cGATGGTTGCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCgcgc                                            <  1:1745769/93‑1 (MQ=255)

                                                 gATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGcc                                           >  2:3301120/1‑93 (MQ=255)

                                                 gATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGcc                                           <  1:4379043/93‑1 (MQ=255)

                                                 gATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGcc                                           >  1:298213/1‑93 (MQ=255)

                                                 gATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑C‑CACCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGcc                                           >  1:1274520/1‑93 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                          <  2:4264300/93‑1 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                          <  1:2587071/93‑1 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACACACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                          >  2:2582377/1‑93 (MQ=255)

                                                   tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                         <  1:4305915/93‑1 (MQ=255)

                                                   tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                         >  2:4059146/1‑93 (MQ=255)

                                                   tGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCt                                         >  2:1593724/1‑93 (MQ=255)

                                                    ggTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACGCACGGATGGCGACGAAGGAATGATTTTTCGCGCCCtt                                        >  2:2764814/1‑93 (MQ=255)

                                                     gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                       <  2:1840405/93‑1 (MQ=255)

                                                      tttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTg                                      >  1:4674453/1‑93 (MQ=255)

                                                         cACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCt                                   <  2:1840666/93‑1 (MQ=255)

                                                          aCCGAGGCTGGCTTTACCTTCCATGCGGGCAACTG‑C‑C‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                  <  2:3968364/93‑1 (MQ=255)

                                                          aCCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                  >  1:2420685/1‑93 (MQ=255)

                                                          aCCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                  <  2:309343/93‑1 (MQ=255)

                                                          aCCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                  >  1:609050/1‑93 (MQ=255)

                                                           ccGAGGCCGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                                 <  1:1669329/93‑1 (MQ=255)

                                                            cGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAac                                <  2:3429841/93‑1 (MQ=255)

                                                             gAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAacc                               >  1:82162/1‑92 (MQ=255)

                                                             gAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAaca                               >  2:4326183/1‑93 (MQ=255)

                                                              aGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAacac                              >  1:1208599/1‑93 (MQ=255)

                                                              aGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAacac                              <  1:3641985/93‑1 (MQ=255)

                                                               ggctggctTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                             >  2:1105034/1‑93 (MQ=255)

                                                                gctggctTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTa                            >  2:1025448/1‑93 (MQ=255)

                                                                   ggctTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                         >  2:590693/1‑93 (MQ=255)

                                                                   ggctTTACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                         <  2:531364/93‑1 (MQ=255)

                                                                      tttACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTACACGTa                      >  1:3040578/1‑93 (MQ=255)

                                                                      tttACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTa                      >  2:3574252/1‑93 (MQ=255)

                                                                      tttACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTa                      <  1:1807236/93‑1 (MQ=255)

                                                                       ttACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                     >  1:3952987/1‑93 (MQ=255)

                                                                       ttACCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                     >  1:3141994/1‑93 (MQ=255)

                                                                         tccTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTacac                   <  1:1138609/92‑1 (MQ=255)

                                                                         aCCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTacac                   <  1:513077/93‑1 (MQ=255)

                                                                         aCCTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAAACCCACGGATGTCGGCGAAGGAATGATTTTTCGCGCCCTGTTGCTAACACTACACGTacac                   >  2:2718937/1‑93 (MQ=255)

                                                                          ccTTCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACt                  <  2:880711/93‑1 (MQ=255)

                                                                             tCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCAt               >  1:2893536/1‑93 (MQ=255)

                                                                             tCCATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCAt               >  2:4531487/1‑93 (MQ=255)

                                                                               cATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa             <  1:347943/93‑1 (MQ=255)

                                                                               cATGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa             >  2:3738269/1‑93 (MQ=255)

                                                                                aTGCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt            >  1:1746922/1‑93 (MQ=255)

                                                                                  gCGGGCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAttt          >  2:466718/1‑93 (MQ=255)

                                                                                    gggCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTAc        >  2:373907/1‑93 (MQ=255)

                                                                                    gggCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTAc        >  1:2306652/1‑93 (MQ=255)

                                                                                     ggCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt       <  1:890979/93‑1 (MQ=255)

                                                                                     ggCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt       <  2:2488830/93‑1 (MQ=255)

                                                                                      gCAAC‑G‑GCC‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtt      <  2:685715/93‑1 (MQ=255)

                                                                                         aC‑G‑G‑C‑C‑CCGGGAAAGGCCCACGGAAGGCGACCAAGTAATGGTTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCt  <  2:2720649/93‑1 (MQ=255)

GCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAAC‑G‑G‑C‑C‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCT  >  NC_000913/1029024‑1029205

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