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.6%	1191.6 / ‑3.6	262	coding (25/1185 nt)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (136/125); minor base T (0/1); total (136/126)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCACA  >  NC_000913/2487111‑2487293

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

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

   aGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaa                                                                                         <  2:207908/93‑1 (MQ=255)

   aGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGACATTTGCTTATGCCTTGaaaa                                                                                         >  2:6834179/1‑93 (MQ=255)

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

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

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

     gTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaag                                                                                       >  1:5389081/1‑92 (MQ=255)

      tCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagc                                                                                      <  1:2631197/93‑3 (MQ=255)

       ccATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagcc                                                                                     <  2:4923729/93‑4 (MQ=255)

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

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

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

          ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGaa                                                                                 <  1:5571043/93‑1 (MQ=255)

           tAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAg                                                                                <  2:2967765/93‑1 (MQ=255)

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

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

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

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                         >  2:1158433/1‑93 (MQ=255)

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                         >  1:5434679/1‑93 (MQ=255)

                   tgtgTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAg                                                                        >  2:4107031/1‑93 (MQ=255)

                   tgtgTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAg                                                                        >  1:4397376/1‑93 (MQ=255)

                   tgtgTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAg                                                                        >  2:7368039/1‑93 (MQ=255)

                    gtgtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGa                                                                       <  2:7743026/93‑1 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGacc                                                                     >  2:4436921/1‑92 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGaca                                                                     <  1:7180577/93‑1 (MQ=255)

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

                          tatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGAAAAAAGACACTTc                                                                 <  2:5003056/93‑1 (MQ=255)

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

                              gATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCa                                                             >  1:7522890/1‑93 (MQ=255)

                              gATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCa                                                             >  1:865410/1‑93 (MQ=255)

                               aTAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAt                                                            <  2:2496294/93‑1 (MQ=255)

                               aTAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAt                                                            >  1:6138271/1‑93 (MQ=255)

                                 aaCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAttt                                                          <  2:5395599/93‑1 (MQ=255)

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

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

                                   cTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTTCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGTCACTTCGGCATTTAt                                                        >  2:7030897/1‑93 (MQ=255)

                                   cTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTAt                                                        >  1:3354432/1‑93 (MQ=255)

                                   cTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTAt                                                        <  1:6555296/93‑1 (MQ=255)

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

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

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

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTT‑TAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                      <  1:6362540/93‑1 (MQ=255)

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

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                   >  1:4197986/1‑93 (MQ=255)

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

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

                                          gCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAAt                                                 <  1:3718429/93‑1 (MQ=255)

                                          gCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAAt                                                 >  2:2316749/1‑93 (MQ=255)

                                          gCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAAt                                                 <  2:167776/93‑1 (MQ=255)

                                           gcTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATc                                                <  2:6717436/92‑1 (MQ=255)

                                            cTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCa                                               <  1:1502086/93‑1 (MQ=255)

                                             ttCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCaa                                              >  2:7865569/1‑93 (MQ=255)

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

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

                                                aaaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                           <  1:5803717/93‑1 (MQ=255)

                                                aaaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                           <  1:3925486/93‑1 (MQ=255)

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAAc                                          <  1:1663052/93‑1 (MQ=255)

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

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

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

                                                    ggCCCTTTGCTTTCATTATTTGTAATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                       <  2:4337947/93‑1 (MQ=255)

                                                     gCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATgt                                      >  1:6224663/1‑92 (MQ=255)

                                                      cccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGat                                     <  1:3617790/93‑1 (MQ=255)

                                                       ccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatt                                    >  1:1307717/1‑93 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                   >  2:1409963/1‑93 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                   >  2:7336938/1‑93 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                   >  1:4786625/1‑93 (MQ=255)

                                                          ttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaat                                 >  2:2664603/1‑93 (MQ=255)

                                                          ttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaat                                 >  2:4878046/1‑93 (MQ=255)

                                                          ttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaat                                 <  1:7886299/93‑1 (MQ=255)

                                                            gCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                               >  1:3269514/1‑93 (MQ=255)

                                                            gCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                               <  2:7764013/93‑1 (MQ=255)

                                                             cTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTAt                              <  2:7511504/93‑1 (MQ=255)

                                                                 cATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaaa                          >  1:7277829/1‑93 (MQ=255)

                                                                 cATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaaa                          <  1:1844606/93‑1 (MQ=255)

                                                                 cATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaaa                          <  1:3522868/93‑1 (MQ=255)

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

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

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

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

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

                                                                     attTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAg                      >  1:4582220/1‑93 (MQ=255)

                                                                     attTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAg                      >  2:1857536/1‑93 (MQ=255)

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

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

                                                                          tCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGAc                 <  1:5218628/93‑1 (MQ=255)

                                                                            aTTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTg               >  2:3239380/1‑93 (MQ=255)

                                                                              ttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGAc             <  1:221371/93‑1 (MQ=255)

                                                                              ttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGAc             <  2:3065796/93‑1 (MQ=255)

                                                                               tGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACa            >  2:4702819/1‑93 (MQ=255)

                                                                                 cTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATa          <  2:5777045/93‑1 (MQ=255)

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa        <  1:2866473/93‑1 (MQ=255)

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

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa        <  2:7295936/93‑1 (MQ=255)

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

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc      <  1:1208829/93‑1 (MQ=255)

                                                                                      tccTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATcc     <  2:6181613/92‑1 (MQ=255)

                                                                                      gCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATcc     <  1:5510411/93‑1 (MQ=255)

                                                                                      gCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATcc     >  1:5091659/1‑93 (MQ=255)

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCACA  >  NC_000913/2487111‑2487293

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