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	.	95.5%	356.3 / 7.8	91	coding (185/681 nt)	ycbF	putative chaperone
Reads supporting (aligned to +/- strand): ref base A (0/2); major base . (62/23); minor base G (1/3); total (63/28)
Fisher's exact test for biased strand distribution p-value = 7.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 110 total aligned reads displayed.

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAAC  >  NC_000913/1003292‑1003474

agaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGa                                                                                             >  1:922855/1‑93 (MQ=255)

   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGagcg                                                                                          >  1:4322616/1‑90 (MQ=255)

   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGagcg                                                                                          >  1:1402341/1‑90 (MQ=255)

   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGagcg                                                                                          >  1:310881/1‑90 (MQ=255)

      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCg                                                                                      >  1:4340909/1‑93 (MQ=255)

      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCg                                                                                      >  2:1431091/1‑93 (MQ=255)

       atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑TGCGTCGa                                                                                     >  2:3858482/1‑93 (MQ=255)

        tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGAc                                                                                    >  1:4519475/1‑93 (MQ=255)

        tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGAc                                                                                    >  1:3166946/1‑93 (MQ=255)

          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAg                                                                                  >  2:2494345/1‑93 (MQ=255)

           cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGt                                                                                 >  1:4349263/1‑93 (MQ=255)

            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTg                                                                                >  2:4662436/1‑93 (MQ=255)

            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTg                                                                                >  1:4538546/1‑93 (MQ=255)

             aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGc                                                                               >  2:3664404/1‑93 (MQ=255)

              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGcc                                                                              >  2:3027765/1‑93 (MQ=255)

               tttACTGACTAATACCAGTCAGAATTCCTGGCTTATTAAAAGAAAAACCACCCGCCCACCGC‑GTGGGGCCGGGGGGG‑AGCGTCGACAGTGCCa                                                                             <  2:2120608/93‑1 (MQ=255)

                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAg                                                                            >  2:236384/1‑93 (MQ=255)

                 tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGc                                                                           >  2:438669/1‑93 (MQ=255)

                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCa                                                                          <  2:780517/93‑1 (MQ=255)

                   cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAc                                                                         >  1:390125/1‑93 (MQ=255)

                   cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGCGCCAGCAc                                                                         >  1:390126/1‑93 (MQ=255)

                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAcc                                                                        >  1:697215/1‑93 (MQ=255)

                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCAcc                                                                        >  1:1540643/1‑93 (MQ=255)

                       cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                      <  1:3143260/93‑1 (MQ=255)

                        tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTa                                                                    >  2:140857/1‑93 (MQ=255)

                         aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑c                                                                  >  1:3938014/1‑93 (MQ=255)

                         aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑c                                                                  >  1:689641/1‑93 (MQ=255)

                          aTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑Ct                                                                 >  1:2916633/1‑93 (MQ=255)

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTg                                                                >  2:4487526/1‑93 (MQ=255)

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTg                                                                >  2:1609327/1‑93 (MQ=255)

                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTgg                                                               >  2:812593/1‑93 (MQ=255)

                             ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGc                                                              <  2:1328127/93‑1 (MQ=255)

                             ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGc                                                              >  2:2714085/1‑93 (MQ=255)

                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGcc                                                             >  1:4395939/1‑93 (MQ=255)

                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGcc                                                             >  2:1858063/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGc                                                           >  1:3191155/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGc                                                           >  2:1203061/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGc                                                           >  1:3862587/1‑93 (MQ=255)

                                gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGc                                                           >  1:1477807/1‑93 (MQ=255)

                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTcc                                                        >  1:260452/1‑93 (MQ=255)

                                    ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCg                                                       >  2:2559991/1‑93 (MQ=255)

                                     gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGc                                                      >  1:1329361/1‑93 (MQ=255)

                                     gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGc                                                      >  1:4097341/1‑93 (MQ=255)

                                        tcgcggCTTAATACTAGAAAACTCAGCCG‑CCACCGCAGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCAc                                                   <  2:2816715/89‑1 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCAc                                                   >  2:1499839/1‑93 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCAc                                                   <  1:3049288/93‑1 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCAc                                                   <  1:678905/93‑1 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCAc                                                   >  1:3733301/1‑93 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCAc                                                   >  2:3017961/1‑93 (MQ=255)

                                         ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACt                                                  <  1:419967/93‑1 (MQ=255)

                                           tGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTa                                                >  1:3776578/1‑93 (MQ=255)

                                           tGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTa                                                >  1:1472881/1‑93 (MQ=255)

                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTAt                                               >  1:3047856/1‑93 (MQ=255)

                                              cTTATTAACAGTAAAACCAACCCCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATtc                                             <  1:3672681/93‑1 (MQ=255)

                                               ttattaACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGGAG‑GTCGAAGAGGGCAACAACATTACCTGGCGCCTCCG‑CACTTAttgt                                            >  1:1386033/1‑91 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATtctc                                           >  2:1458197/1‑93 (MQ=255)

                                                tattaACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCACTA‑CTGGCCGCTCCGCCACTTCTtctc                                           >  1:3226824/1‑93 (MQ=255)

                                                tataaACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATtctc                                           <  2:3832027/93‑1 (MQ=255)

                                                  ttaACAGTAAAATCACCCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCt                                         <  1:2158763/93‑1 (MQ=255)

                                                  ttaACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCt                                         <  1:1539799/93‑1 (MQ=255)

                                                   taACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCTGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTg                                        >  2:3969334/1‑93 (MQ=255)

                                                   taACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTg                                        >  2:438653/1‑93 (MQ=255)

                                                    aaCAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGa                                       >  2:2758651/1‑93 (MQ=255)

                                                     aCAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑CGAGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGaa                                      >  2:2778527/1‑93 (MQ=255)

                                                       aGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                    >  2:2566322/1‑93 (MQ=255)

                                                       aGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                    >  2:3426750/1‑93 (MQ=255)

                                                       aGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                    <  1:1364401/93‑1 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                   >  1:705135/1‑93 (MQ=255)

                                                        gTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                   >  1:4233860/1‑93 (MQ=255)

                                                          aaaaTCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGTCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTTAAGGCCg                                 >  2:45345/1‑93 (MQ=255)

                                                          aaaaTCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                                 <  1:1379888/93‑1 (MQ=255)

                                                           aaaTCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                                >  2:562665/1‑93 (MQ=255)

                                                            aaTCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                               <  1:4437661/93‑1 (MQ=255)

                                                             aTCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                              <  1:1396331/93‑1 (MQ=255)

                                                             aTCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                              >  2:627653/1‑93 (MQ=255)

                                                              tCAACCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                             >  1:4685851/1‑93 (MQ=255)

                                                                aaCCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGGCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCCGTACGa                           >  2:2602896/1‑93 (MQ=255)

                                                                 aCCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGACGCACGGTAAGAc                          >  1:1889461/1‑93 (MQ=255)

                                                                 aCCGCCCAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                          >  1:3635689/1‑93 (MQ=255)

                                                                   cGCCCAACAC‑GTAGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                        <  2:4633600/93‑1 (MQ=255)

                                                                     cccAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                      >  1:4374466/1‑93 (MQ=255)

                                                                     cccAACGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                      >  1:1530873/1‑93 (MQ=255)

                                                                      ccAACGG‑GGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                      <  2:2694615/93‑1 (MQ=255)

                                                                         aCGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                  >  1:3228682/1‑93 (MQ=255)

                                                                         aCGC‑GTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                  >  2:1336058/1‑93 (MQ=255)

                                                                            c‑gTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc               >  2:4224343/1‑93 (MQ=255)

                                                                               tggggCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc             <  1:2616788/91‑1 (MQ=255)

                                                                                 gggCCGGGGGGG‑AGCGTCGACAGTGCCAGCGCCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt           <  2:3735575/93‑1 (MQ=255)

                                                                                 gggCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt           <  1:584023/93‑1 (MQ=255)

                                                                                  ggCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg          >  1:849476/1‑93 (MQ=255)

                                                                                   gCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc         <  2:3556813/93‑1 (MQ=255)

                                                                                   gCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc         <  2:1921642/93‑1 (MQ=255)

                                                                                    gcGGGGGGGGAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc         <  1:236384/92‑1 (MQ=255)

                                                                                     ccggggggGAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGACCGGTACGACAGGCACGTTGCGCTTGCt        >  1:696040/3‑93 (MQ=255)

                                                                                       ggggggGAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga      <  1:143848/93‑1 (MQ=255)

                                                                                        ggggggAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag     >  1:1849297/1‑93 (MQ=255)

                                                                                        ggggggAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag     >  1:287495/1‑93 (MQ=255)

                                                                                        ggggggAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag     >  1:3560471/1‑93 (MQ=255)

                                                                                          ggggAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa   >  1:2646852/1‑93 (MQ=255)

                                                                                           ggAAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAAc  <  1:2574698/93‑1 (MQ=255)

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC‑GTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTA‑CTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAAC  >  NC_000913/1003292‑1003474

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