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	.	90.6%	612.1 / 49.2	171	coding (185/681 nt)	ycbF	putative chaperone
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (105/50); minor base G (10/6); total (115/56)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

Only 100 of 192 total aligned reads displayed.

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1003292‑1003473

atattCGTTATCTTTTTTCCTGACTAATACCAGCCAGGACTCCTGGCTTATTAACACAAAAAACACCCGCCCAGCGCGTGGGGCCGGGGGGGa                                                                                           <  2:3600475/89‑1 (MQ=255)

     caatatTTATTTTACTGACTAAACCCAGTCAGGAATCCTGGCTTATAAACAGAAAAATCAACCGCCCAACGGGTTGGGCCGGGGGGG‑AGCGTc                                                                                     <  2:7737267/91‑1 (MQ=255)

     cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTc                                                                                     >  2:257961/1‑93 (MQ=255)

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

       atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGa                                                                                   >  2:6244756/1‑93 (MQ=255)

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

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

          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTTGGCCGGGGGGG‑AGCGTCGACAg                                                                                >  1:7421800/1‑93 (MQ=255)

          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑TGCGTCGACAg                                                                                >  1:6011624/1‑93 (MQ=255)

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

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

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

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

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

               tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCa                                                                           >  1:5030502/1‑93 (MQ=255)

               tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCa                                                                           >  1:4712361/1‑93 (MQ=255)

                  tcTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCa                                                                        <  2:2765781/92‑1 (MQ=255)

                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCa                                                                        >  2:4134614/1‑93 (MQ=255)

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

                     gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCa                                                                     >  1:7600195/1‑93 (MQ=255)

                      aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCAt                                                                    <  2:1301883/93‑1 (MQ=255)

                       cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCAtt                                                                   >  2:1917906/1‑93 (MQ=255)

                       cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGCTGGGCCGGGGGGG‑AGAGTCGACAGTGCCAGCACCAtt                                                                   >  1:883190/1‑93 (MQ=255)

                        tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:2725877/1‑93 (MQ=255)

                        tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:1728663/1‑93 (MQ=255)

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

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

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTg                                                               <  1:4765791/93‑1 (MQ=255)

                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTAGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCCCCCTTACTg                                                               >  2:3303935/1‑93 (MQ=255)

                             ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGc                                                             <  1:1962022/93‑1 (MQ=255)

                              cAGTCAGGAATCCTGGTTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGca                                                            <  2:3901021/93‑2 (MQ=255)

                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            <  2:2659472/93‑1 (MQ=255)

                                  cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCCACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTTGTCGCTc                                                        >  2:1043738/1‑93 (MQ=255)

                                  cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTc                                                        >  1:4122158/1‑93 (MQ=255)

                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       >  1:2549845/1‑93 (MQ=255)

                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       >  2:5826626/1‑93 (MQ=255)

                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       >  2:1242885/1‑93 (MQ=255)

                                      aaTCGTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                    <  1:4888038/93‑1 (MQ=255)

                                        tCCTGGCTTATTAACAGTAAAAGCAACCGCCCGCGGCGTT‑GGTGGGGGGGGGAGCCTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                  <  2:1336866/93‑1 (MQ=255)

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

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

                                           tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                               >  1:3071351/1‑93 (MQ=255)

                                           tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                               >  1:2046520/1‑93 (MQ=255)

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

                                             gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCCTTACTGGCCGCTCCGCCACTTAtt                                             >  1:5221090/1‑93 (MQ=255)

                                             gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt                                             >  2:5961133/1‑93 (MQ=255)

                                             gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt                                             >  1:5138603/1‑93 (MQ=255)

                                              cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                            >  2:7576596/1‑93 (MQ=255)

                                               tttttaACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                            <  2:6667507/90‑1 (MQ=255)

                                               ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                           >  1:2339218/1‑93 (MQ=255)

                                                tattaACAGTAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGTCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                          <  2:6316620/93‑1 (MQ=255)

                                                 attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑TGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                         >  1:4564105/1‑93 (MQ=255)

                                                  ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                        <  1:6766771/93‑1 (MQ=255)

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

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

                                                    aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                      >  1:848830/1‑93 (MQ=255)

                                                     aCAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                     <  2:5097155/93‑1 (MQ=255)

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

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

                                                        gaaaaaCCACCCCCCAAAGGCGTTGGG‑CGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                  <  1:2482857/91‑1 (MQ=255)

                                                           aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                               >  1:5330849/1‑93 (MQ=255)

                                                             ataaaCCGCCCAGCGCGTTGGGCCGGGGAGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  2:1175168/90‑1 (MQ=255)

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

                                                             aTCAACCGCCCAACGCGGTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  1:5526337/93‑1 (MQ=255)

                                                              tCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCGCTTATTCTCCTGAAGCCCGGTAc                            >  2:5634882/1‑93 (MQ=255)

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

                                                              tCAACCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                            <  1:6573887/93‑1 (MQ=255)

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

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

                                                                 aCCGCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                         <  2:3651727/93‑1 (MQ=255)

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

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

                                                                   cGCCCAACGCGTTGGGCCGGGGGGG‑AGCGCCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGCCAg                       >  2:432494/1‑93 (MQ=255)

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

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      <  1:6531070/93‑1 (MQ=255)

                                                                    gCCCAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      <  2:5512832/93‑1 (MQ=255)

                                                                    ccccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      <  2:4799676/92‑1 (MQ=255)

                                                                     cccAACGCGTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTCAAGCCCGGTACGACAGGc                     >  1:5745982/1‑93 (MQ=255)

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

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

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

                                                                          ggcgTGGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGtt                <  1:4332957/92‑1 (MQ=255)

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

                                                                           gcgTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               >  2:3859199/1‑93 (MQ=255)

                                                                            cgTTGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc              <  1:7297467/93‑1 (MQ=255)

                                                                             gTTGGGCCGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCTACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc              >  2:2908483/1‑93 (MQ=255)

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

                                                                              ttGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc            >  2:2180780/1‑93 (MQ=255)

                                                                              ttGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc            >  2:4046879/1‑93 (MQ=255)

                                                                              ttGGG‑CGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc            <  1:3586901/93‑1 (MQ=255)

                                                                               tGGGCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           <  2:7923642/93‑1 (MQ=255)

                                                                                  gCCGGGGGGG‑AGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        <  1:3265442/93‑1 (MQ=255)

                                                                                    ccggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt       <  2:2064435/91‑1 (MQ=255)

                                                                                     cggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg      >  1:4412796/2‑93 (MQ=255)

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

                                                                                        ggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgaga   >  2:3071945/1‑93 (MQ=255)

                                                                                        ggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTACTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgaga   >  2:7377625/1‑93 (MQ=255)

                                                                                         ggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  >  1:2171082/1‑93 (MQ=255)

                                                                                         ggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  >  2:728418/1‑93 (MQ=255)

                                                                                         ggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  <  1:3859200/93‑1 (MQ=255)

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1003292‑1003473

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