BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,226,457	A→G	N339S (AAT→AGT)	ygjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,226,457	0	A	G	98.4%	561.0 / ‑0.8	182	N339S (AAT→AGT)	ygjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (64/115); total (66/116)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 202 total aligned reads displayed.

TGAAACCACCTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAATGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTA  >  NC_000913/3226365‑3226549

tGAAACCACCTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGaaaaa                                                                                              <  2:930673/93‑1 (MQ=255)

    accaccTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGcc                                                                                          <  2:2268034/93‑1 (MQ=255)

    accaccTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGcc                                                                                          >  2:2263075/1‑93 (MQ=255)

    accaccTCTTTCAGGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGcc                                                                                          <  2:566550/93‑1 (MQ=255)

      caccTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCctc                                                                                        <  1:3217865/93‑1 (MQ=255)

      caccTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAATGCctc                                                                                        <  2:3863090/93‑1 (MQ=255)

        ccTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCtctc                                                                                      >  2:2356155/1‑93 (MQ=255)

         ctctTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCa                                                                                     <  2:129403/93‑1 (MQ=255)

         ctctTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCa                                                                                     <  1:3906288/93‑1 (MQ=255)

          tctTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCAc                                                                                    <  1:105981/93‑1 (MQ=255)

          tctTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCAc                                                                                    <  2:1170658/93‑1 (MQ=255)

          tctTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCAc                                                                                    >  2:1134123/1‑93 (MQ=255)

            tttCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGa                                                                                  <  2:4191838/93‑1 (MQ=255)

              tCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATc                                                                                >  1:4622421/1‑93 (MQ=255)

              tCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATc                                                                                <  1:310537/93‑1 (MQ=255)

                aTGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATccc                                                                              >  2:1971594/1‑93 (MQ=255)

                 tGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCa                                                                             <  1:4549071/93‑1 (MQ=255)

                 tGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCa                                                                             >  1:908414/1‑93 (MQ=255)

                 tGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCa                                                                             <  1:1777498/93‑1 (MQ=255)

                   cAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGc                                                                           <  1:1168692/93‑1 (MQ=255)

                     gTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGa                                                                         <  2:3260696/93‑1 (MQ=255)

                     gTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAATGCCTCTCACGATCCCAGCGa                                                                         >  2:1677540/1‑93 (MQ=255)

                      tttGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGTCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGAc                                                                        <  2:4337872/93‑1 (MQ=255)

                        tGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCa                                                                      >  1:4329097/1‑93 (MQ=255)

                         ggTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCaa                                                                     <  2:2539457/93‑1 (MQ=255)

                         ggTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCaa                                                                     >  2:315079/1‑93 (MQ=255)

                          gTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCaaa                                                                    <  1:2340829/93‑1 (MQ=255)

                           tATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAAt                                                                   <  2:1462111/93‑1 (MQ=255)

                            aTCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATTAAAAGTGCCTCTCACGATACCCGCGACCAAATc                                                                  >  2:485806/1‑93 (MQ=255)

                            aTCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATc                                                                  <  1:1325434/93‑1 (MQ=255)

                            aTCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATc                                                                  <  1:2957615/93‑1 (MQ=255)

                             tCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCt                                                                 <  1:3620391/93‑1 (MQ=255)

                              cTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCtt                                                                >  2:1922793/1‑93 (MQ=255)

                               tGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTc                                                               >  1:2085991/1‑93 (MQ=255)

                                  tACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTg                                                            <  2:4364136/93‑1 (MQ=255)

                                   aCGAATACGGCGGCGTACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGa                                                           <  1:1664740/93‑1 (MQ=255)

                                    cGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGAt                                                          <  2:3994723/93‑1 (MQ=255)

                                    cGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGAt                                                          >  2:3484828/1‑93 (MQ=255)

                                    cGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGAt                                                          <  2:4456467/93‑1 (MQ=255)

                                        tACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCaa                                                      <  1:4404141/93‑1 (MQ=255)

                                        tACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCaa                                                      <  1:1379495/93‑1 (MQ=255)

                                         aCGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAat                                                     <  1:2859530/93‑1 (MQ=255)

                                            gcggcggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATAtt                                                  <  2:609628/93‑1 (MQ=255)

                                             cggcggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATAttt                                                 <  2:3499190/93‑1 (MQ=255)

                                              ggcggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTc                                                <  1:2692568/93‑1 (MQ=255)

                                              ggcggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTc                                                >  1:3000720/1‑93 (MQ=255)

                                                gggACGTTTCCCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGAGGCAATATTTCTa                                              <  2:724846/92‑1 (MQ=255)

                                                cggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTa                                              >  1:4115101/1‑93 (MQ=255)

                                                 ggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCtcc                                             <  2:3733287/93‑3 (MQ=255)

                                                 ggACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTAc                                             <  2:3324957/93‑1 (MQ=255)

                                                  gACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACt                                            <  1:4502332/93‑1 (MQ=255)

                                                   aCGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTg                                           <  1:891453/93‑1 (MQ=255)

                                                     gTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGt                                         >  1:1849812/1‑93 (MQ=255)

                                                     gTTTCGCCACGCCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGt                                         <  2:3510376/93‑1 (MQ=255)

                                                      tttCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGta                                        >  1:3578382/1‑93 (MQ=255)

                                                      tttCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGta                                        <  2:1102137/93‑1 (MQ=255)

                                                      tttCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGta                                        >  2:4051884/1‑93 (MQ=255)

                                                          gCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGtaataa                                    >  1:4670935/1‑93 (MQ=255)

                                                          gCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGtaataa                                    <  1:1528161/93‑1 (MQ=255)

                                                          gCCACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGtaataa                                    <  1:2279653/93‑1 (MQ=255)

                                                            cACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAgg                                  <  2:864723/93‑1 (MQ=255)

                                                            cACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAgg                                  >  2:256173/1‑93 (MQ=255)

                                                            cACGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAgg                                  <  2:684631/93‑1 (MQ=255)

                                                             aCGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGa                                 <  2:2225552/93‑1 (MQ=255)

                                                             aCGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGa                                 <  2:2165236/93‑1 (MQ=255)

                                                             aCGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGa                                 >  1:3143365/1‑93 (MQ=255)

                                                              cGGCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGAc                                <  2:2531153/93‑1 (MQ=255)

                                                               ggCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACt                               <  1:484591/93‑1 (MQ=255)

                                                               ggCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACt                               >  2:2836842/1‑93 (MQ=255)

                                                               ggCTTTTGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACt                               <  2:4475748/93‑1 (MQ=255)

                                                                gCTTTTGATAGCCGCTTCTACATGAAAAATGCCTCTCACGATGCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTg                              >  1:3303397/1‑93 (MQ=255)

                                                                 cTTTTGATAGCCGCTTCTACATGAAATGTGCCTCGCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGt                             <  1:992261/93‑1 (MQ=255)

                                                                  ttttGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTa                            >  2:4631214/1‑93 (MQ=255)

                                                                    ttGATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAAt                          <  1:1286946/93‑1 (MQ=255)

                                                                      gATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAAttt                        >  2:4641532/1‑93 (MQ=255)

                                                                      gATAGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAAttt                        <  1:1919391/93‑1 (MQ=255)

                                                                         ggCCGCTGCTACATGAAAAGTCCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCAt                     <  2:2590265/92‑1 (MQ=255)

                                                                         aGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCAt                     >  2:3092117/1‑93 (MQ=255)

                                                                         aGCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCAt                     <  1:2473643/93‑1 (MQ=255)

                                                                          gCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATg                    <  2:1572276/93‑1 (MQ=255)

                                                                          gCCGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCATGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATg                    <  1:3293337/93‑1 (MQ=255)

                                                                           ccGCTTCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGa                   <  1:3690088/93‑1 (MQ=255)

                                                                               ttCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAAt               >  1:423148/1‑93 (MQ=255)

                                                                               ttCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAAt               >  1:3793819/1‑93 (MQ=255)

                                                                                tCTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATa              <  1:2569371/93‑1 (MQ=255)

                                                                                 cTACATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATaa             <  1:317699/93‑1 (MQ=255)

                                                                                    cATGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATaaaaa          <  1:3560365/93‑1 (MQ=255)

                                                                                     aTGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAAc         <  1:2221263/93‑1 (MQ=255)

                                                                                     aTGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAAc         >  2:2121643/1‑93 (MQ=255)

                                                                                     aTGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAAc         >  1:3033107/1‑93 (MQ=255)

                                                                                      tGAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACt        <  2:4443346/93‑1 (MQ=255)

                                                                                       taaaagTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTa       <  2:2651192/92‑1 (MQ=255)

                                                                                       gAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTa       <  2:390997/93‑1 (MQ=255)

                                                                                       gAAAAGTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTa       <  1:2138703/93‑1 (MQ=255)

                                                                                        aaaagTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTAt      >  1:1302347/1‑93 (MQ=255)

                                                                                        aaaagTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTAt      >  2:771062/1‑93 (MQ=255)

                                                                                         aaagTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTAtt     <  1:1423537/93‑1 (MQ=255)

                                                                                         aaagTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTAtt     >  1:1315032/1‑93 (MQ=255)

                                                                                         aaagTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTAtt     <  1:1009191/93‑1 (MQ=255)

                                                                                            aTGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTa  >  2:4142835/1‑93 (MQ=255)

TGAAACCACCTCTTTCATGCAGTTTGGTATCTGGTACGAATACGGCGGCGGACGTTTCGCCACGGCTTTTGATAGCCGCTTCTACATGAAAAATGCCTCTCACGATCCCAGCGACCAAATCTTCCTGATGCAATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTA  >  NC_000913/3226365‑3226549

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