BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,158,202	+C	coding (961/975 nt)	ygiRQ ←	Fe‑S oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,158,201	1	.	C	99.5%	692.5 / ‑4.5	218	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (98/119); total (98/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 241 total aligned reads displayed.

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑C‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

 cGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑c                                                                                         >  1:472824/1‑93 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cc                                                                                        <  1:6730820/93‑2 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cc                                                                                        <  2:6271095/93‑2 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cc                                                                                        >  1:1981742/1‑92 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cc                                                                                        >  2:6336197/1‑92 (MQ=255)

   cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cct                                                                                       <  1:3779198/93‑3 (MQ=255)

    tGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cctg                                                                                      <  2:5370519/93‑4 (MQ=255)

    tGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cctg                                                                                      >  2:6701479/1‑90 (MQ=255)

    tGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cctg                                                                                      <  1:4108994/93‑4 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCTTCAATGTGCGGGCACAGCC‑TGggtt                                                                                  >  2:6717675/1‑93 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggtt                                                                                  >  1:3341983/1‑93 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttg                                                                                 >  2:1579012/1‑93 (MQ=255)

          tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttgg                                                                                <  1:818724/93‑1 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttggtt                                                                              <  1:2869592/93‑1 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGgggtggtt                                                                              >  1:1417262/1‑93 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGttt                                                                             >  2:2544233/1‑93 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTc                                                                            >  1:408240/1‑93 (MQ=255)

                cGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgt                                                                          <  2:5962398/93‑1 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtg                                                                         >  1:284781/1‑93 (MQ=255)

                   tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtggt                                                                       >  2:4442638/1‑93 (MQ=255)

                    aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtggtg                                                                      >  1:5073257/1‑93 (MQ=255)

                    aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtggtg                                                                      <  1:5879639/93‑1 (MQ=255)

                     gAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGc                                                                     >  2:1472816/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCa                                                                    <  1:4000125/93‑1 (MQ=255)

                       aaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAg                                                                   >  2:173148/1‑93 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGaa                                                                 >  2:4061300/1‑93 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAAt                                                                >  2:3270352/1‑93 (MQ=255)

                           ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATa                                                               >  2:2522743/1‑93 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAc                                                              >  2:2362807/1‑93 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAc                                                              >  1:4908546/1‑93 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAc                                                              >  1:2123820/1‑93 (MQ=255)

                             aaTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAcg                                                             <  2:2647145/93‑1 (MQ=255)

                              aTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAcgc                                                            <  1:4761615/93‑1 (MQ=255)

                               tAGCAGGCGGGTTGATCCACACATAGCGGTGGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑CTTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           <  1:1260539/93‑1 (MQ=255)

                                  cAGGCGGGTTGATCCACACATAGCGGTCTCCGTGTTTTTGCATCAAGGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                                        <  2:2032627/93‑1 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCTCTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      >  2:4689491/1‑93 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      >  1:3838949/1‑93 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGCCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGCTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      >  2:4185247/1‑93 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:2374653/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     >  1:3433309/1‑93 (MQ=255)

                                      cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    <  1:3982303/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   <  1:1089536/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   <  2:3156919/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   >  1:1061348/1‑93 (MQ=255)

                                         gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                 <  1:6339834/93‑1 (MQ=255)

                                          ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                <  2:4972896/93‑1 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                              <  2:5153451/93‑1 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          >  1:150090/1‑93 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          >  1:2673632/1‑93 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          >  2:3109023/1‑93 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          >  2:102166/1‑93 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCa                                         >  1:4278214/1‑93 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCa                                         <  2:6513505/93‑1 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCa                                         >  1:2451459/1‑93 (MQ=255)

                                                   acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                       >  2:855094/1‑93 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      >  2:375772/1‑93 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  1:5136153/93‑1 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  1:901106/93‑1 (MQ=255)

                                                    aaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  1:5317433/92‑1 (MQ=255)

                                                        gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATc                                  <  1:4528360/93‑1 (MQ=255)

                                                          ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATAGGCGTACAGCACTTTATCGc                                >  1:36601/1‑93 (MQ=255)

                                                            tCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGccc                              >  1:5858168/1‑93 (MQ=255)

                                                            tCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGccc                              <  1:6643543/93‑1 (MQ=255)

                                                             cccggtgtTTTTGCACAAATGCGCGGGCACAGAC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCt                             <  1:6237451/89‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  2:6005862/93‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  2:2829702/93‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  2:2928163/93‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           >  2:6606142/1‑93 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  1:3552811/93‑1 (MQ=255)

                                                               ccGTGCTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  2:2730005/93‑1 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                         >  2:2450859/1‑93 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                         >  1:4168802/1‑93 (MQ=255)

                                                                  tgtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                        >  2:4675393/1‑93 (MQ=255)

                                                                  tgtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                        >  2:3024962/1‑93 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                       <  1:6433088/93‑1 (MQ=255)

                                                                    tttttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACttt                      <  2:1449591/93‑1 (MQ=255)

                                                                     ttttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtttt                     >  2:590072/1‑92 (MQ=255)

                                                                     ttttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                     <  1:923785/93‑1 (MQ=255)

                                                                      tttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtct                    <  1:2090939/93‑1 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                   <  1:4831998/93‑1 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                   <  2:1302822/93‑1 (MQ=255)

                                                                        tGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCg                  >  2:4352796/1‑93 (MQ=255)

                                                                          cATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGac                >  1:1315949/1‑92 (MQ=255)

                                                                          cATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaa                >  1:6700365/1‑93 (MQ=255)

                                                                           aTCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa               <  1:6174081/93‑1 (MQ=255)

                                                                           aTCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa               >  2:183944/1‑93 (MQ=255)

                                                                           aTCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCCCGaaa               <  2:4481561/93‑1 (MQ=255)

                                                                            tCAATGCGCGGGCACAGCC‑TGGGATGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaag              >  1:1935911/1‑93 (MQ=255)

                                                                              aaTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            <  2:789116/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            <  2:5575326/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            >  2:5736100/1‑93 (MQ=255)

                                                                              aaTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGCACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            >  2:1542467/1‑93 (MQ=255)

                                                                               aTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaag           <  2:608349/93‑1 (MQ=255)

                                                                               aTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaag           <  1:5553902/93‑1 (MQ=255)

                                                                                 gcgcGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc         >  1:164937/1‑93 (MQ=255)

                                                                                     gggCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagca     <  2:4384394/93‑1 (MQ=255)

                                                                                     gggCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACATTATCGCCCTTCACTTTCTCGAAAGAAGgccgca     >  1:2783119/1‑93 (MQ=255)

                                                                                      ggCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCaa    >  1:1831595/1‑93 (MQ=255)

                                                                                        cacaGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAaca  <  2:1151967/93‑1 (MQ=255)

                                                                                        cacaGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAaca  <  1:3210867/93‑1 (MQ=255)

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑C‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

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