BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,487,208	(A)_6→5	coding (20/1185 nt)	yfdE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,487,203	0	A	.	99.2%	1752.9 / 3.9	388	coding (25/1185 nt)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/3); new base . (204/180); total (205/183)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 437 total aligned reads displayed.

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCACA  >  NC_000913/2487111‑2487293

gAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGa                                                                                            <  2:7281782/93‑1 (MQ=255)

 aaaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaa                                                                                           >  1:5869813/1‑93 (MQ=255)

 aaaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGACCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaa                                                                                           >  2:6840832/1‑93 (MQ=255)

  aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                          >  2:3391391/1‑93 (MQ=255)

    ggTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                        <  2:3227761/93‑1 (MQ=255)

    ggTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                        >  2:1984793/1‑93 (MQ=255)

      tCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagc                                                                                      >  1:7082370/1‑91 (MQ=255)

       ccATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagcc                                                                                     >  2:9889137/1‑90 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                    >  2:9000998/1‑89 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                    >  2:7909568/1‑89 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                    <  2:6087447/93‑5 (MQ=255)

         aTTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGa                                                                                  >  1:7994752/1‑93 (MQ=255)

          ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGaa                                                                                 >  2:4491138/1‑93 (MQ=255)

          ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGaa                                                                                 >  2:4902386/1‑93 (MQ=255)

           tAAGGACATGTGTCATATCGATAACTAATACGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAg                                                                                <  1:6028174/93‑1 (MQ=255)

           tAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAg                                                                                >  1:4935439/1‑93 (MQ=255)

            aaGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAgg                                                                               >  2:3070262/1‑93 (MQ=255)

            aaGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAgg                                                                               <  2:2127942/93‑1 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                              <  2:5066825/93‑1 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                              <  1:5173704/93‑1 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                              >  1:9498034/1‑93 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                              <  1:8490410/93‑1 (MQ=255)

              ggACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTg                                                                             <  1:266258/93‑1 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                            >  2:2681929/1‑93 (MQ=255)

                     tgtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGac                                                                      >  2:6025514/1‑93 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGaca                                                                     >  1:2798233/1‑93 (MQ=255)

                       tCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGacac                                                                    >  1:1506687/1‑93 (MQ=255)

                        cATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACt                                                                   <  1:1148627/93‑1 (MQ=255)

                        cATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACt                                                                   >  1:10023538/1‑93 (MQ=255)

                        cATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACt                                                                   <  2:3482003/93‑1 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                  >  2:5700811/1‑93 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                  <  2:6429610/93‑1 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                  <  1:368543/93‑1 (MQ=255)

                          tatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTc                                                                 >  2:8602303/1‑93 (MQ=255)

                          tatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTc                                                                 <  2:8289976/93‑1 (MQ=255)

                           atCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCg                                                                <  2:411459/93‑1 (MQ=255)

                            tCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCgg                                                               <  1:230254/93‑1 (MQ=255)

                             cGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGc                                                              >  2:8788645/1‑93 (MQ=255)

                                tAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAtt                                                           <  2:7400590/93‑1 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                         <  2:5814950/93‑1 (MQ=255)

                                    taataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATg                                                       >  1:734940/1‑93 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATg                                                       <  2:1057146/93‑1 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                      <  2:7759017/93‑1 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                      >  1:429871/1‑93 (MQ=255)

                                       taaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaa                                                    <  1:9399814/93‑1 (MQ=255)

                                       taaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaa                                                    <  1:7909567/93‑1 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                   <  2:4781604/93‑1 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                   <  2:2536785/93‑1 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGGGACAAGACACTTCGGCATTTATGTaaa                                                   >  1:1721095/1‑93 (MQ=255)

                                         aGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaaa                                                  >  1:7227803/1‑93 (MQ=255)

                                            cTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCa                                               >  2:7588649/1‑93 (MQ=255)

                                             ttCAAACGGCCCTTTGCTTTCATTATTTGTGATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCaa                                              >  1:5505329/1‑93 (MQ=255)

                                             ttCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTCACAAGACACTTCGGCATTTATGTAAAATCaa                                              >  2:9819735/1‑93 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                             <  2:6718898/93‑1 (MQ=255)

                                               cAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATa                                            <  1:4552151/93‑1 (MQ=255)

                                               cAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAGTTATGTAAAATCAATa                                            <  1:5430479/93‑1 (MQ=255)

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTTTGCCTTG‑AAAAAGCGGAAGGGGAAAAGAACTTTCGCCTTTTTTGAAAAATCAATaca                                          >  1:4092131/1‑91 (MQ=255)

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAAc                                          >  2:1535710/1‑93 (MQ=255)

                                                  aCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACa                                         <  2:7808729/93‑1 (MQ=255)

                                                  aCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACa                                         >  1:7644812/1‑93 (MQ=255)

                                                    ggCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                       >  1:7713489/1‑93 (MQ=255)

                                                     gCCTTTTGTTTTCATTTTTTGTCATTGGTTTATG‑CTTAAAAAAAGCCGAAGGTAAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGa                                      <  1:1913858/93‑1 (MQ=255)

                                                     gCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGa                                      >  1:8690781/1‑93 (MQ=255)

                                                      cccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGat                                     >  2:5461420/1‑93 (MQ=255)

                                                       ccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatt                                    <  2:7193711/93‑1 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                  >  2:3526640/1‑93 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                  >  2:4217387/1‑93 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                  <  1:216247/93‑1 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                  >  1:5640677/1‑93 (MQ=255)

                                                              tttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATg                             <  1:3171732/93‑1 (MQ=255)

                                                              tttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATg                             >  2:2435521/1‑93 (MQ=255)

                                                              tttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATg                             >  1:6525426/1‑93 (MQ=255)

                                                              tttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCAAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATg                             <  2:4353186/93‑1 (MQ=255)

                                                                 cATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACCCTTCGGAATTTATGTAAAATCAATAACATGATTAATTATGgaa                          >  1:626537/1‑93 (MQ=255)

                                                                  attattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAc                         <  2:3971829/93‑1 (MQ=255)

                                                                   ttattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAcc                        >  2:3893498/1‑93 (MQ=255)

                                                                   ttattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAcc                        <  2:2299296/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                       >  1:1586153/1‑93 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                       <  2:4632149/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                       <  2:9349716/93‑1 (MQ=255)

                                                                      tttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGt                     >  1:4698829/1‑93 (MQ=255)

                                                                      tttGTCAGTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGt                     <  1:8379891/93‑1 (MQ=255)

                                                                       ttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTc                    <  1:9100552/93‑1 (MQ=255)

                                                                        tGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCg                   <  2:1747711/93‑1 (MQ=255)

                                                                         gTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGa                  <  2:8167783/93‑1 (MQ=255)

                                                                          tCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGAc                 >  2:943585/1‑93 (MQ=255)

                                                                           cATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACt                <  2:5263032/93‑1 (MQ=255)

                                                                            aTTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTg               <  2:1843187/93‑1 (MQ=255)

                                                                             tttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa              >  1:6705510/1‑93 (MQ=255)

                                                                             tttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa              <  2:3689793/93‑1 (MQ=255)

                                                                               tGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACa            <  1:490508/93‑1 (MQ=255)

                                                                                 cTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATa          >  1:3065853/1‑93 (MQ=255)

                                                                                  ttATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaa         <  2:1633385/93‑1 (MQ=255)

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa        >  1:8076882/1‑93 (MQ=255)

                                                                                    aTGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAAt       >  1:3887800/1‑93 (MQ=255)

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc      >  1:3422743/1‑93 (MQ=255)

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc      >  1:987129/1‑93 (MQ=255)

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc      >  2:9520585/1‑93 (MQ=255)

                                                                                       ccTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCca    >  2:2724402/1‑93 (MQ=255)

                                                                                         ttGAAAAAA‑CCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCcaca  <  2:8035194/93‑1 (MQ=255)

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCACA  >  NC_000913/2487111‑2487293

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