BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,784,344	(T)_8→7	100%	pseudogene (668/1374 nt)	ypjC ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,784,337	0	T	.	67.2%	150.7 / 185.3	119	pseudogene (675/1374 nt)	ypjC	pseudogene
Reads supporting (aligned to +/- strand): ref base T (17/22); new base . (38/42); total (55/64)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 139 total aligned reads displayed.

TAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGA  >  NC_000913/2784238‑2784434

tAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCt                                                                                                   >  1:725209/1‑100 (MQ=255)

 aaaaaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtt                                                                                                  >  2:186209/1‑100 (MQ=255)

 aaaaaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtt                                                                                                  >  1:2407692/1‑100 (MQ=255)

      cTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCttttttt                                                                                             <  1:1705613/100‑1 (MQ=255)

       ttATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttt                                                                                            >  1:1011461/1‑100 (MQ=255)

       ttATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttt                                                                                            >  1:2887963/1‑100 (MQ=255)

       ttATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttc                                                                                            <  2:218179/100‑2 (MQ=255)

        tataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTc                                                                                           <  2:2544679/100‑1 (MQ=255)

         ataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttctt                                                                                          <  2:87030/100‑4 (MQ=255)

         ataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttctt                                                                                          >  1:2265632/1‑97 (MQ=255)

          taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttcttg                                                                                         <  1:3003786/100‑5 (MQ=255)

          taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttcttg                                                                                         <  1:1012008/100‑5 (MQ=255)

          taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttcttg                                                                                         <  1:630653/100‑5 (MQ=255)

          taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttat                                                                                             >  2:398896/1‑94 (MQ=255)

          taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCtt                                                                                         <  2:2742538/100‑1 (MQ=255)

          taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCtt                                                                                         <  2:2601843/100‑1 (MQ=255)

          taGCAACACTCATAAAAATTAAACACAAGCGCTGAAGACCTTGTATTTTTATCAAATTATATCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCtt                                                                                         <  1:1498654/100‑1 (MQ=255)

           aGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTg                                                                                        <  1:1148462/100‑1 (MQ=255)

           aGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTg                                                                                        <  2:2446983/100‑1 (MQ=255)

           aGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGa                                                                                       >  2:1815684/1‑100 (MQ=255)

            gCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttttttctt                                                                                          >  1:2456291/1‑94 (MQ=255)

                   tCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTAt                                                                               <  1:1422299/100‑1 (MQ=255)

                    cATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGCATATCTa                                                                                >  2:1724566/1‑98 (MQ=255)

                     aTAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATgg                                                                             <  2:159452/100‑1 (MQ=255)

                     aTAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATgg                                                                             >  1:845812/1‑100 (MQ=255)

                      tAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATgg                                                                             >  2:148654/1‑100 (MQ=255)

                      tAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATggg                                                                            <  2:2319771/100‑1 (MQ=255)

                      tAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATggg                                                                            >  1:2929429/1‑100 (MQ=255)

                        aaaaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTTTTTTTTTTTACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATggg                                                                            >  1:1004328/1‑99 (MQ=255)

                        aaaaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATgggg                                                                           <  2:142598/100‑1 (MQ=255)

                        aaaaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGgt                                                                          <  1:1319399/100‑1 (MQ=255)

                          aaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGgtg                                                                         >  2:2236394/1‑100 (MQ=255)

                              aaaCAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCt                                                                    <  2:2168002/100‑1 (MQ=255)

                              aaaCAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCt                                                                    >  2:1153164/1‑100 (MQ=255)

                               aaCAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCtt                                                                   <  1:1537301/100‑1 (MQ=255)

                               aaCAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCtt                                                                   >  2:325083/1‑100 (MQ=255)

                                 cAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCttt                                                                  >  2:2569523/1‑100 (MQ=255)

                                 cAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCttt                                                                  <  2:2038806/100‑1 (MQ=255)

                                 cAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTg                                                                 <  2:657163/100‑1 (MQ=255)

                                 cAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTg                                                                 <  2:184614/100‑1 (MQ=255)

                                 cAGAAGCTCTCAAGACCTTGTATTTTTACCCAATACTACCCCCTTGAACTCATTCACAAACGCATC‑TCTTTTTCTTGAATATCTATGGGGTGTTCCTTTg                                                                 <  1:963635/100‑1 (MQ=255)

                                  aGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTg                                                                 <  1:2019742/100‑1 (MQ=255)

                                  aGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGa                                                                >  2:858235/1‑100 (MQ=255)

                                   gAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGAc                                                               <  2:516480/100‑1 (MQ=255)

                                     aGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGActc                                                             <  2:1951895/100‑1 (MQ=255)

                                       ctctGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTg                                                           <  2:1021095/100‑1 (MQ=255)

                                        tctGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGa                                                          <  1:1418289/100‑1 (MQ=255)

                                         ctGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGAc                                                         <  2:294522/100‑1 (MQ=255)

                                          tGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACa                                                        >  2:3060967/1‑100 (MQ=255)

                                          tGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACa                                                        >  1:3061577/1‑100 (MQ=255)

                                           gAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTTAATATCTATGGGGTGTTCCTTTGACTCTGAc                                                         >  1:220546/1‑98 (MQ=255)

                                             aGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAg                                                     <  1:344014/100‑1 (MQ=255)

                                             aGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAg                                                     <  2:146123/100‑1 (MQ=255)

                                              gACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGa                                                    <  2:1675332/100‑1 (MQ=255)

                                               aCCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTATGACTCTTACAAAg                                                     >  2:670545/1‑99 (MQ=255)

                                                ccTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGAt                                                   >  2:233307/1‑100 (MQ=255)

                                                   tGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGaa                                                >  1:450376/1‑99 (MQ=255)

                                                    gTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGaaa                                               <  2:2076737/100‑1 (MQ=255)

                                                      aTTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGAGTGTTCCTTTTACTCTGACAAAGATGAAAt                                              >  1:1729644/1‑99 (MQ=255)

                                                      aTTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGATGTTCCTTTGACTCTGACAAAGATGAAAtc                                             >  1:2451706/1‑99 (MQ=255)

                                                       tttttATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtct                                            <  2:1506981/100‑1 (MQ=255)

                                                       tttttATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtct                                            >  2:1340992/1‑100 (MQ=255)

                                                       tttttATCAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtctc                                           >  2:2289919/1‑100 (MQ=255)

                                                            aTCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTTTTTTTTATCTATGGGGAGTTCCTTTGACTCTGACAAAGATGAAATCTCATCt                                       >  1:1302824/1‑100 (MQ=255)

                                                              cAAATAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAGATTTCATCtttg                                    <  2:1198713/100‑4 (MQ=255)

                                                               aaaTAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTAttt                                   >  2:294191/1‑100 (MQ=255)

                                                               aaaTAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTAttt                                   >  2:669234/1‑100 (MQ=255)

                                                               aaaTAATACCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTAttt                                   >  2:1302668/1‑100 (MQ=255)

                                                                    ataCCCCCTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGttt                              <  1:1368323/100‑1 (MQ=255)

                                                                       ctcacTTGAACTCCTGCACCAACGCCTC‑TTTTTTTCTTGAATATGTATGGGGTGTTCCTTTGGGTCTGACAAAGATGAAATCTCATCTATTTTGttttta                           <  1:1146638/96‑3 (MQ=255)

                                                                       cccccTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGtttt                             <  1:2633065/99‑1 (MQ=255)

                                                                       cccccTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCa                           >  2:2459375/1‑100 (MQ=255)

                                                                         cccTTGAACTCATTCACCAACGCATGTTTTTTTTCTTTAATATCTATGGGGTGTTCTTTTGATTCTGATAAAGATGAAATCTCATCTATTTTTTTTTc                            >  1:883933/1‑98 (MQ=255)

                                                                         cccTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCat                          >  1:2669977/1‑100 (MQ=255)

                                                                         cccTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCat                          <  2:502092/99‑1 (MQ=255)

                                                                          ccTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCa                           <  1:663935/97‑1 (MQ=255)

                                                                           cTTGAACTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATg                       <  1:2764420/100‑1 (MQ=255)

                                                                             tGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGa                      <  2:2877755/100‑1 (MQ=255)

                                                                                aCTCATTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTa                  >  1:684238/1‑100 (MQ=255)

                                                                                 cTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTa                  >  2:1370371/1‑100 (MQ=255)

                                                                                   cATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTAc                 <  2:23994/99‑1 (MQ=255)

                                                                                    aTTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGt               >  2:24058/1‑100 (MQ=255)

                                                                                    aTTCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTg              <  2:378479/100‑1 (MQ=255)

                                                                                      tCACAAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGAt            <  1:1152476/100‑1 (MQ=255)

                                                                                       cacaAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGAtt           >  2:1066429/1‑100 (MQ=255)

                                                                                       cacaAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGAtt           >  2:2526767/1‑100 (MQ=255)

                                                                                        acaAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTc          <  2:437533/100‑1 (MQ=255)

                                                                                         caAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTc          <  1:2236381/100‑1 (MQ=255)

                                                                                         caAACGCATC‑TTTTTTTCTTGAATATCTATGTGCTGTTGCTTTGGCTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTcc         <  1:2577658/100‑1 (MQ=255)

                                                                                         caAACGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTcc         <  2:2407705/100‑1 (MQ=255)

                                                                                          aaaCGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCa        <  2:2265645/100‑1 (MQ=255)

                                                                                           aaCGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCAt       >  1:2588569/1‑100 (MQ=255)

                                                                                            aCGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATa      <  2:2718450/100‑1 (MQ=255)

                                                                                             cGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATa      <  2:3026361/100‑1 (MQ=255)

                                                                                             cGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATa      <  1:2071205/99‑1 (MQ=255)

                                                                                             cGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATAc     <  2:1256044/100‑1 (MQ=255)

                                                                                             cGCATC‑TTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATAc     <  1:1815673/100‑1 (MQ=255)

                                                                                                 tCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGa  <  2:2887980/100‑1 (MQ=255)

                                                                                                 tCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGa  <  2:1011466/100‑1 (MQ=255)

                                                                                                  cTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAg   >  1:308596/1‑98 (MQ=255)

TAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGA  >  NC_000913/2784238‑2784434

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