BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,003,385	Δ1 bp	coding (186/681 nt)	ycbF →	putative chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,003,384	0	A	.	93.6%	674.3 / 30.7	172	coding (185/681 nt)	ycbF	putative chaperone
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (123/38); minor base G (2/9); total (125/47)
Fisher's exact test for biased strand distribution p-value = 1.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 193 total aligned reads displayed.

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGA  >  NC_000913/1003292‑1003472

agaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGa                                                                                          <  1:6147058/93‑1 (MQ=255)

    tCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagcgt                                                                                      >  1:6304795/1‑89 (MQ=255)

    tCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagcgc                                                                                      >  1:3446006/1‑89 (MQ=255)

     tgATATTTTTTTTACTGACTAATCCCAGTCAGGAATCCTGTTTTATTAACAGAAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGG‑AGCGTc                                                                                    <  2:1620939/92‑1 (MQ=255)

      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCg                                                                                   >  2:2740968/1‑93 (MQ=255)

      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCg                                                                                   >  1:2930698/1‑93 (MQ=255)

      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCg                                                                                   >  1:2246444/1‑93 (MQ=255)

        tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGAc                                                                                 >  1:6095747/1‑93 (MQ=255)

         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACa                                                                                >  1:4263525/1‑93 (MQ=255)

         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACa                                                                                >  2:2127465/1‑93 (MQ=255)

          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAg                                                                               >  2:3393088/1‑93 (MQ=255)

          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAg                                                                               >  1:5530334/1‑93 (MQ=255)

           cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAgg                                                                              >  2:1766248/1‑92 (MQ=255)

           cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGt                                                                              >  1:176111/1‑93 (MQ=255)

            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTg                                                                             >  1:6244955/1‑93 (MQ=255)

              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGcc                                                                           >  1:2458108/1‑93 (MQ=255)

                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑GGCGTCGACAGTGCCAg                                                                         >  2:2843356/1‑93 (MQ=255)

                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAcc                                                                     >  1:4576644/1‑93 (MQ=255)

                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAcc                                                                     >  2:1041411/1‑93 (MQ=255)

                      aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCAt                                                                   >  1:4895586/1‑93 (MQ=255)

                         aaaaCCAGTCAGAAATCCTGGCTTATAAACAGAAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTAc                                                                <  1:4724657/90‑1 (MQ=255)

                         aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTAc                                                                >  1:272171/1‑93 (MQ=255)

                         aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTAc                                                                <  1:2079029/93‑1 (MQ=255)

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTg                                                              >  1:1080774/1‑93 (MQ=255)

                            aCCAGTCAGAAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTgg                                                             <  1:4013098/93‑1 (MQ=255)

                             ccAGTCAGGAATCCTGTCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGc                                                            <  2:767842/93‑1 (MQ=255)

                             ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGc                                                            >  2:5353135/1‑93 (MQ=255)

                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGcc                                                           >  1:792584/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                         >  1:5022397/1‑93 (MQ=255)

                                     gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc                                                    >  2:1159432/1‑93 (MQ=255)

                                     gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc                                                    >  2:5476332/1‑93 (MQ=255)

                                      aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                   >  1:327269/1‑93 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                 >  1:5230508/1‑93 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                 >  2:4298493/1‑93 (MQ=255)

                                         ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                >  1:3665415/1‑93 (MQ=255)

                                         ccTGGCTTATTAACAGGAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                >  2:3185681/1‑93 (MQ=255)

                                          cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                               >  1:5701055/1‑93 (MQ=255)

                                          cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                               <  1:2122585/93‑1 (MQ=255)

                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                             >  2:6207380/1‑93 (MQ=255)

                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                             >  1:590060/1‑93 (MQ=255)

                                              tttattaACAGAAAAATCAACCGCCCATCGCGTTGGCCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                           <  2:2711383/92‑1 (MQ=255)

                                              cgtattaACATAAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                           <  2:3956193/91‑1 (MQ=255)

                                              cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                           <  1:4030254/93‑1 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                         >  2:1606522/1‑93 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                         >  1:3902918/1‑93 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                         >  1:2975446/1‑93 (MQ=255)

                                                  ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                       >  2:2627309/1‑93 (MQ=255)

                                                   taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                      >  1:378865/1‑93 (MQ=255)

                                                    aaCAGTAAAATCAACCGCCCAACGCGATGGGCCGGGGGGG‑AGTGTCGAAAGGGCCAGCACCAGTACTGGCCGCACCGCCCCTTATTCTCCTGa                                     >  2:2228108/1‑93 (MQ=255)

                                                     aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                    <  1:5267187/93‑1 (MQ=255)

                                                      cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAg                                   >  2:877104/1‑93 (MQ=255)

                                                       aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑TGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTGTTCTCCTGAAGc                                  >  2:3105011/1‑93 (MQ=255)

                                                       aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                  >  1:5742083/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                 >  2:1929862/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                 >  2:1152307/1‑93 (MQ=255)

                                                         tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                <  1:3071003/93‑1 (MQ=255)

                                                          aaagTCAACCGGCCAGCGCTTTTGGCC‑GGGGGGAAGCGTCGACAGTTCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                               <  2:6131425/93‑1 (MQ=255)

                                                          aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                               >  2:4082825/1‑93 (MQ=255)

                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                            >  1:233486/1‑93 (MQ=255)

                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                            >  1:1124739/1‑93 (MQ=255)

                                                              tCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                           >  2:717663/1‑93 (MQ=255)

                                                              tCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                           >  2:949317/1‑93 (MQ=255)

                                                               cAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                          >  1:4515003/1‑93 (MQ=255)

                                                                acccGCCCACCGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                         <  1:2620960/91‑1 (MQ=255)

                                                                aaCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACgc                         >  1:1617449/1‑92 (MQ=255)

                                                                 aCCGCCCAACGCGTTGGGCCGGGGGGG‑TGCGTCGACAGTGTCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                        >  2:3721477/1‑93 (MQ=255)

                                                                 aCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                        >  1:1031377/1‑93 (MQ=255)

                                                                 aCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                        >  1:4920417/1‑93 (MQ=255)

                                                                  ccccccAACGGGTGGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                       <  1:4639931/90‑1 (MQ=255)

                                                                  ccGCCCACGGTGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                       <  1:6000030/93‑1 (MQ=255)

                                                                  ccGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                       >  2:5189371/1‑93 (MQ=255)

                                                                  ccGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                       <  1:4196160/93‑1 (MQ=255)

                                                                   cGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                      >  1:47099/1‑93 (MQ=255)

                                                                   cGCCAAAGGCGTGGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                      <  2:626445/93‑1 (MQ=255)

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                     >  2:4447395/1‑93 (MQ=255)

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTAATGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAAAgg                     >  2:300255/1‑93 (MQ=255)

                                                                     cccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                    >  1:1478895/1‑93 (MQ=255)

                                                                     cccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                    >  1:407529/1‑93 (MQ=255)

                                                                     cccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCACCGCCACTTATTCTCCTGAAGCCCGGTGCGACAGGc                    >  2:5870135/1‑93 (MQ=255)

                                                                      ccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                   <  1:1055661/93‑1 (MQ=255)

                                                                      ccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                   >  2:5053515/1‑93 (MQ=255)

                                                                      ccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                   <  2:2001541/93‑1 (MQ=255)

                                                                      ccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                   >  2:226445/1‑93 (MQ=255)

                                                                       cAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                  >  1:942494/1‑93 (MQ=255)

                                                                       cAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                  >  1:71076/1‑93 (MQ=255)

                                                                       cAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                  >  1:1178654/1‑93 (MQ=255)

                                                                       cAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                  >  1:4047010/1‑93 (MQ=255)

                                                                       cAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                  <  1:1934725/93‑1 (MQ=255)

                                                                        aaCGCGTTGGGCCGGGGGGG‑GGCGTCGACAGTGCCAGCACCATTACTGGTCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                 >  2:5564612/1‑93 (MQ=255)

                                                                        aaCGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                 >  2:2696411/1‑93 (MQ=255)

                                                                          cgggTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGtt               <  2:4863546/90‑1 (MQ=255)

                                                                          cgcgTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGtt               <  2:4253458/93‑1 (MQ=255)

                                                                             gTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg            >  1:2109734/1‑93 (MQ=255)

                                                                              ttGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc           <  2:6078585/93‑1 (MQ=255)

                                                                                 gggcGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCTACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt         <  1:3375802/90‑1 (MQ=255)

                                                                                 ggCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg        >  1:4402245/1‑93 (MQ=255)

                                                                                   gcGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc       <  1:2416841/92‑1 (MQ=255)

                                                                                      ggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga    >  2:2491885/1‑93 (MQ=255)

                                                                                      ggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga    <  2:4645506/93‑1 (MQ=255)

                                                                                       ggggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgtg   >  2:2602257/1‑91 (MQ=255)

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGA  >  NC_000913/1003292‑1003472

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: `‑`