BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,486,295	T→C	D446G (GAC→GGC)	yjgR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,486,295	0	T	C	90.6%	956.3 / 68.8	372	D446G (GAC→GGC)	yjgR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (18/17); new base C (156/181); total (174/198)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 385 total aligned reads displayed.

TTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCG  >  NC_000913/4486203‑4486384

   tGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCAt                                                                                        >  2:3935263/1‑93 (MQ=255)

    ggCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATc                                                                                       <  2:4812590/93‑1 (MQ=255)

       cATTGTTTGCACCACGCCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCAc                                                                                    <  1:8201957/93‑1 (MQ=255)

       cATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCAc                                                                                    <  2:1313525/93‑1 (MQ=255)

        aTTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCAcc                                                                                   >  1:8899797/1‑93 (MQ=255)

        aTTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCAcc                                                                                   <  2:9357931/93‑1 (MQ=255)

        aTTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCAcc                                                                                   <  2:4503453/93‑1 (MQ=255)

         ttGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCg                                                                                  >  1:2650950/1‑93 (MQ=255)

         ttGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCAccc                                                                                  >  1:8171837/1‑92 (MQ=255)

         ttGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCg                                                                                  >  1:3116538/1‑93 (MQ=255)

         gtGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCg                                                                                  >  2:6161415/2‑93 (MQ=255)

             ttGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTAc                                                                              <  2:8198729/93‑1 (MQ=255)

               gCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCACTCCACCAAGAATGCCGTCATCCACCGCTACct                                                                            >  1:5285969/1‑93 (MQ=255)

               gCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACct                                                                            >  1:1484921/1‑93 (MQ=255)

                caccacaccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACctc                                                                           >  1:9383804/1‑93 (MQ=255)

                 accacaccaTCTTTCTTCCTGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCCATCCACCAAGAATGCCGCCATCCACCGCTACctct                                                                          >  2:3577987/1‑93 (MQ=255)

                 accacaccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACctct                                                                          >  2:8709029/1‑93 (MQ=255)

                   caccccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCttt                                                                        <  2:3776436/93‑1 (MQ=255)

                   cacaccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCttt                                                                        >  2:5171938/1‑93 (MQ=255)

                     caccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTcc                                                                      <  2:5826338/93‑1 (MQ=255)

                      cccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTccc                                                                     <  1:4475328/92‑1 (MQ=255)

                      cccaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTccc                                                                     <  1:6712838/92‑1 (MQ=255)

                      accaTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTccc                                                                     >  2:9710660/1‑93 (MQ=255)

                         aTCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCttt                                                                  <  2:5558050/93‑1 (MQ=255)

                          tctttcttCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTc                                                                 >  2:7485461/1‑93 (MQ=255)

                          tctttcttCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTc                                                                 >  1:7669127/1‑93 (MQ=255)

                            tttcttCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTcgc                                                               >  1:6588650/1‑93 (MQ=255)

                            tttcttCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTcgc                                                               >  1:1900916/1‑93 (MQ=255)

                            tttcttCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTcgc                                                               <  1:7226333/93‑1 (MQ=255)

                             ttcttcCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTcgcg                                                              >  2:4287181/1‑93 (MQ=255)

                              tcttcCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCgg                                                             <  1:1656592/93‑1 (MQ=255)

                                ttcCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGgga                                                           <  1:9080345/93‑1 (MQ=255)

                                 tcCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGggag                                                          <  1:4979968/93‑1 (MQ=255)

                                   ccgccgCGTGGTCCGGTAGTGCCGAACAAAATATCTTTCAATCCACCAACAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGggagga                                                        <  1:9986652/93‑1 (MQ=255)

                                   ccgccgCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGggagga                                                        <  1:7837824/93‑1 (MQ=255)

                                     gccgcGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGatt                                                      <  1:989287/93‑1 (MQ=255)

                                     gccgcGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGatt                                                      <  2:1758419/93‑1 (MQ=255)

                                     gccgcGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGatt                                                      <  2:6896098/93‑1 (MQ=255)

                                      ccgcGTGGTCCGGTAGTGTCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGatta                                                     <  2:2244614/93‑1 (MQ=255)

                                       cgcgTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGattat                                                    >  1:8388590/1‑93 (MQ=255)

                                        gcgTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGattatt                                                   >  2:341638/1‑93 (MQ=255)

                                        gcgTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGattatt                                                   <  1:8550001/93‑1 (MQ=255)

                                        gcgTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGattatt                                                   <  1:9127494/93‑1 (MQ=255)

                                         cgTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTAttt                                                  <  2:4704477/93‑1 (MQ=255)

                                          gTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTAtttt                                                 >  2:8513720/1‑93 (MQ=255)

                                               ccGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTtgctgc                                            >  1:1356012/1‑93 (MQ=255)

                                               ccGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTtgctgc                                            <  1:6971658/93‑1 (MQ=255)

                                               ccGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTtgctgc                                            <  1:5025680/93‑1 (MQ=255)

                                                 ggTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTc                                          >  2:2459996/1‑93 (MQ=255)

                                                  gTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCg                                         >  1:10027134/1‑93 (MQ=255)

                                                  gTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCg                                         <  1:5139757/93‑1 (MQ=255)

                                                       gCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACt                                    <  1:4148808/93‑1 (MQ=255)

                                                       gCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCATTTCGCGGGAGGATTATTTTGCTGCTCGGTACt                                    >  2:2967424/1‑93 (MQ=255)

                                                         cGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTgg                                  >  1:3521336/1‑93 (MQ=255)

                                                         cGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTgg                                  <  1:6276427/93‑1 (MQ=255)

                                                         cGAACAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTgg                                  <  1:8277379/93‑1 (MQ=255)

                                                           aaCAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGcc                                <  1:1182429/93‑1 (MQ=255)

                                                           aaCAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGcc                                <  2:2600506/93‑1 (MQ=255)

                                                            aCAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCt                               >  2:910498/1‑93 (MQ=255)

                                                            aCAAAATATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGATACTGGCCt                               >  1:5603740/1‑93 (MQ=255)

                                                              aaaaTATCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGa                             >  1:7454784/1‑93 (MQ=255)

                                                                   atCCTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGcc                        >  1:2808160/1‑93 (MQ=255)

                                                                      cTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCttt                     <  1:6801563/93‑1 (MQ=255)

                                                                      cTTCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCttt                     <  2:6546936/93‑1 (MQ=255)

                                                                       ttCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCtttt                    >  1:7978004/1‑93 (MQ=255)

                                                                       ttCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCtttt                    <  2:602650/93‑1 (MQ=255)

                                                                       ttCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCtttt                    <  2:3193448/93‑1 (MQ=255)

                                                                       ttCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCtttt                    <  2:8272383/93‑1 (MQ=255)

                                                                        tCAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCttttt                   >  2:3048546/1‑93 (MQ=255)

                                                                         cAATCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTATTg                  <  2:7239423/93‑1 (MQ=255)

                                                                          aaTCCTCCAAGACTGCACCCATCCACCGCTACCTCTTTCCCTGTCGCGGGAGTATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGc                 <  2:1618581/93‑1 (MQ=255)

                                                                          aaTCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGc                 >  2:1488361/1‑93 (MQ=255)

                                                                          aaTCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGc                 >  1:3825342/1‑93 (MQ=255)

                                                                          aaTCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGACTTTTCGc                 >  2:2385958/1‑93 (MQ=255)

                                                                           aTCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCTCGGGAGGATGATTTTGGTGCTCGGTACTGGCCTGAAAGCCTTTTTGCa                >  2:7751605/1‑93 (MQ=255)

                                                                            tCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCTGTACTGGCCTGAAAAACTTTTTTCaa               >  1:4845821/1‑93 (MQ=255)

                                                                            tCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCaa               >  2:7988947/1‑93 (MQ=255)

                                                                            tCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCaa               >  1:8680512/1‑93 (MQ=255)

                                                                            tCCACCAAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCaa               <  2:2889332/93‑1 (MQ=255)

                                                                              caccaAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACa             <  2:3118902/93‑1 (MQ=255)

                                                                               accaAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAt            >  2:7079138/1‑93 (MQ=255)

                                                                                ccaAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAtc           <  1:6224945/93‑1 (MQ=255)

                                                                                ccaAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAtc           <  2:7878377/93‑1 (MQ=255)

                                                                                 caAGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAtct          <  1:2514045/93‑1 (MQ=255)

                                                                                 caAGAATGCCGCCATCCACCGCCACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAtct          >  1:5810655/1‑93 (MQ=255)

                                                                                  aaGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAtctc         <  1:106850/93‑1 (MQ=255)

                                                                                  aaGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACAtctc         >  1:2231701/1‑93 (MQ=255)

                                                                                   aGAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCa        <  1:2324275/93‑1 (MQ=255)

                                                                                    gAATGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCAt       <  2:2026603/93‑1 (MQ=255)

                                                                                     aaggccgccATCCACCGCTACCTCTTTCCCTTTGGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATa      <  2:88791/90‑1 (MQ=255)

                                                                                     aaTGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATa      <  2:5866682/93‑1 (MQ=255)

                                                                                     aaTGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATa      <  2:4953742/93‑1 (MQ=255)

                                                                                     aaTGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATa      >  1:9686753/1‑93 (MQ=255)

                                                                                      aTGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAg     <  2:2226615/93‑1 (MQ=255)

                                                                                      aTGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAg     <  1:6401393/93‑1 (MQ=255)

                                                                                      aTGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAg     >  1:2023677/1‑93 (MQ=255)

                                                                                       tGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAgg    >  1:3027293/1‑93 (MQ=255)

                                                                                       tGCCGCCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAgg    <  2:6662438/93‑1 (MQ=255)

                                                                                        gccgccATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGc   >  1:5480123/1‑93 (MQ=255)

                                                                                        gCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTGCTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGc   <  1:3740032/93‑1 (MQ=255)

TTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAATATCCTTCAATCCACCAAGAATGCCGTCATCCACCGCTACCTCTTTCCCTTTCGCGGGAGGATTATTTTGCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCG  >  NC_000913/4486203‑4486384

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