BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	235,988	T→C	V5V (GTA→GTG)	rnhA ←	RNase HI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	235,988	0	T	C	99.7%	1146.7 / NA	327	V5V (GTA→GTG)	rnhA	RNase HI
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (169/156); minor base G (0/1); total (170/157)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

AAAGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGC  >  NC_000913/235898‑236077

 aaGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCAcc                                                                                        <  2:5389393/93‑1 (MQ=255)

  aGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCt                                                                                       >  2:8343335/1‑93 (MQ=255)

  aGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCt                                                                                       <  1:8368509/93‑1 (MQ=255)

  aGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCt                                                                                       <  1:228978/93‑1 (MQ=255)

    gtgttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGt                                                                                     <  2:2506688/90‑1 (MQ=255)

     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAGAATTTCCACCTGtt                                                                                    >  1:3549262/1‑93 (MQ=255)

     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                    >  2:5784592/1‑93 (MQ=255)

     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                    >  1:4869134/1‑93 (MQ=255)

     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                    <  2:1105560/93‑1 (MQ=255)

     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                    >  1:105743/1‑93 (MQ=255)

      tttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGttt                                                                                   >  2:7103950/1‑93 (MQ=255)

      tttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGttt                                                                                   <  1:4159405/93‑1 (MQ=255)

        tctcGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTaa                                                                                 <  1:5101373/93‑1 (MQ=255)

          tcGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGc                                                                               >  1:7605001/1‑93 (MQ=255)

           cgcgTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCa                                                                              >  2:484964/1‑93 (MQ=255)

           cgcgTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCa                                                                              <  1:454663/93‑1 (MQ=255)

              gTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCAtct                                                                           <  1:4641776/93‑1 (MQ=255)

               tCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCAtctc                                                                          >  1:2835856/1‑93 (MQ=255)

               gccGTGATAGCGGAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCGACCTGTTTAAGCAtctc                                                                          <  1:3595490/92‑1 (MQ=255)

                ccGCGATAGCGTAAAATAGCGCCTTAACCCCCAGGTCCTTGATGGCCCAGACACGAACCATCGGGTAAAATTTCCACCTTTTTAAGCATctct                                                                         <  2:5991389/93‑1 (MQ=255)

                ccGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATctct                                                                         <  1:451195/93‑1 (MQ=255)

                 cgcgATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTg                                                                        >  2:6545999/1‑93 (MQ=255)

                  gcgATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTgg                                                                       <  1:6193434/93‑1 (MQ=255)

                  gcgATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTgg                                                                       <  2:1396821/93‑1 (MQ=255)

                   cgATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGt                                                                      <  2:2459094/93‑1 (MQ=255)

                    gATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCATCTGTTTAAGCATCTCTGGTa                                                                     >  1:2353151/1‑93 (MQ=255)

                    gATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTa                                                                     <  2:1789961/93‑1 (MQ=255)

                     aTAGCGTAAAATAGCGCCGTAACCCGCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAg                                                                    <  1:615012/93‑1 (MQ=255)

                     aTAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAg                                                                    <  2:5290558/93‑1 (MQ=255)

                     aTAGCGTAAAATAGCGCCGTAACCCCCAGGCCGTGGATTGCCCAGCGACGATCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAg                                                                    <  2:4129517/93‑1 (MQ=255)

                      tAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGa                                                                   >  2:4280188/1‑93 (MQ=255)

                       aGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGAc                                                                  >  1:6699360/1‑93 (MQ=255)

                       aGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGAc                                                                  >  2:7369306/1‑93 (MQ=255)

                       aGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGAc                                                                  >  1:1189837/1‑93 (MQ=255)

                        gCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACt                                                                 >  1:3253958/1‑93 (MQ=255)

                           tAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTcc                                                              >  1:1735853/1‑93 (MQ=255)

                             aaaTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTg                                                            <  1:6396183/93‑1 (MQ=255)

                              aaTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGt                                                           <  2:7884904/93‑1 (MQ=255)

                               aTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTGCTGTa                                                          >  1:8696585/1‑93 (MQ=255)

                                tAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTaa                                                         >  2:3546090/1‑93 (MQ=255)

                                 aGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAAt                                                        >  1:2024319/1‑93 (MQ=255)

                                  gcgcCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAAtt                                                       <  2:3091713/93‑1 (MQ=255)

                                   cgcCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTg                                                      <  1:939292/93‑1 (MQ=255)

                                     ccGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGaa                                                    >  2:7497414/1‑93 (MQ=255)

                                      cGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAAt                                                   <  2:4905808/93‑1 (MQ=255)

                                         aaCCCCCAGGTCCAGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGa                                                >  1:5804899/1‑93 (MQ=255)

                                           cccccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAAc                                              >  2:4498399/1‑93 (MQ=255)

                                           cccccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAAc                                              <  2:5422910/93‑1 (MQ=255)

                                           cccccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAAc                                              <  2:7658440/93‑1 (MQ=255)

                                           cccccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAAc                                              >  1:5486351/1‑93 (MQ=255)

                                             cccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTg                                            <  2:88808/93‑1 (MQ=255)

                                              ccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGt                                           >  1:6039127/1‑93 (MQ=255)

                                              ccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGt                                           >  1:3465565/1‑93 (MQ=255)

                                                aGGTCCTGGATTGCCCAGACACGAACCATCGGCGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATAGAACTGTaa                                         >  2:1469870/1‑93 (MQ=255)

                                                 ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaa                                        >  2:7419282/1‑93 (MQ=255)

                                                 ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaa                                        >  2:283329/1‑93 (MQ=255)

                                                 ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaa                                        >  1:5278605/1‑93 (MQ=255)

                                                  gTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaaa                                       >  1:3981620/1‑93 (MQ=255)

                                                  gTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaaa                                       <  2:3762023/93‑1 (MQ=255)

                                                   tCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAAc                                      <  1:7990274/93‑1 (MQ=255)

                                                   tCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAAc                                      <  1:3852651/93‑1 (MQ=255)

                                                    ccTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACg                                     >  2:2745401/1‑93 (MQ=255)

                                                     cTGGACTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGa                                    <  2:7752156/93‑1 (MQ=255)

                                                      tGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGAc                                   <  2:8142614/93‑1 (MQ=255)

                                                       ggATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACa                                  <  2:764048/93‑1 (MQ=255)

                                                       ggATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACa                                  <  1:4799566/93‑1 (MQ=255)

                                                       ggATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACa                                  <  2:5555538/93‑1 (MQ=255)

                                                        gATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACaa                                 >  1:6648044/1‑93 (MQ=255)

                                                        gATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACaa                                 >  2:6250428/1‑93 (MQ=255)

                                                        gATGGCCCAGCCACGACCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACaa                                 <  2:1975897/93‑1 (MQ=255)

                                                          ttGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGt                               >  1:650224/1‑93 (MQ=255)

                                                               cAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGAc                          <  1:691059/93‑1 (MQ=255)

                                                                aGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACa                         >  1:7908456/1‑93 (MQ=255)

                                                                   cacGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATaa                      >  2:5025038/1‑93 (MQ=255)

                                                                    acGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATaaa                     >  1:7454134/1‑93 (MQ=255)

                                                                     cGAACCATCGGTGAAAATTTCCGCCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAAt                    <  2:857522/93‑1 (MQ=255)

                                                                      gAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATg                   <  2:4017433/93‑1 (MQ=255)

                                                                      gAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATg                   >  1:6953783/1‑93 (MQ=255)

                                                                        aCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAc                 <  2:4551900/93‑1 (MQ=255)

                                                                        aCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAc                 <  2:3357271/93‑1 (MQ=255)

                                                                        aCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAc                 <  2:6192005/93‑1 (MQ=255)

                                                                        aCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAc                 <  2:6828833/93‑1 (MQ=255)

                                                                         ccATCGGTGAAAATTTCCCCCTGTTTAAGCATTTCTGGGAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAcc                >  1:5035629/1‑93 (MQ=255)

                                                                         ccATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAcc                >  1:7031933/1‑93 (MQ=255)

                                                                          cATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGATCg               >  1:849031/1‑93 (MQ=255)

                                                                           aTCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGc              >  1:1368909/1‑93 (MQ=255)

                                                                            tCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCt             <  1:1464474/93‑1 (MQ=255)

                                                                             cGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTa            <  1:2531000/93‑1 (MQ=255)

                                                                               gTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATg          <  1:7715157/93‑1 (MQ=255)

                                                                                tGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAGCGACAAGTCTGACATAAATGACCGCTATGa         <  2:4700072/93‑1 (MQ=255)

                                                                                 gAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAg        <  2:7285491/93‑1 (MQ=255)

                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc       >  2:1885671/1‑93 (MQ=255)

                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc       >  2:6971931/1‑93 (MQ=255)

                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc       <  2:6389703/93‑1 (MQ=255)

                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc       <  1:4686394/93‑1 (MQ=255)

                                                                                   aaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCa      <  1:1917765/93‑1 (MQ=255)

                                                                                      tttCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTg   >  1:2210311/1‑93 (MQ=255)

                                                                                      tttCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTg   <  1:808519/93‑1 (MQ=255)

                                                                                      tttCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTg   >  1:5223956/1‑93 (MQ=255)

                                                                                       ttCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGc  <  2:2315034/93‑1 (MQ=255)

AAAGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGC  >  NC_000913/235898‑236077

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