BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,567,922	A→G	F198L (TTC→CTC)	xasA ←	acid sensitivity protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,567,922	0	A	G	99.2%	776.3 / ‑3.4	240	F198L (TTC→CTC)	xasA	acid sensitivity protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (112/125); total (113/127)
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.26e-01

Only 100 of 249 total aligned reads displayed.

TCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATG  >  NC_000913/1567832‑1568011

tcatttcattGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGt                                                                                         <  1:5770792/93‑1 (MQ=255)

tcatttcattGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGt                                                                                         <  2:4228943/93‑1 (MQ=255)

 catttcattGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCATGGTGCCCACTTTAGAGAGGTc                                                                                        <  1:1464395/93‑1 (MQ=255)

 catttcattGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTc                                                                                        <  2:3500564/93‑1 (MQ=255)

   tttcattGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAg                                                                                      >  1:2595022/1‑93 (MQ=255)

      cattGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGga                                                                                   <  2:1116010/93‑1 (MQ=255)

       attGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGgaa                                                                                  <  1:553732/93‑1 (MQ=255)

          gACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGgaagaa                                                                               <  1:1082337/93‑1 (MQ=255)

          gACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGgaagaa                                                                               <  2:3564835/93‑1 (MQ=255)

            cGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAgg                                                                             <  2:937008/93‑1 (MQ=255)

             gTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGt                                                                            <  1:2052194/93‑1 (MQ=255)

              tGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTc                                                                           <  2:1245044/93‑1 (MQ=255)

              tGGATTGTGGATGTTTCTACGCTCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTTCTTACTTTAGAGAGTTCAGGGAAGAAGGTc                                                                           <  1:4859317/93‑1 (MQ=255)

               gggTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCt                                                                          <  2:4888671/93‑1 (MQ=255)

                ggTTGCGGATGCTTCTACGCGCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCtt                                                                         <  1:344415/93‑1 (MQ=255)

                  ttGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCg                                                                       <  2:1368757/93‑1 (MQ=255)

                  ttGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCg                                                                       >  2:5580050/1‑93 (MQ=255)

                   tGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGa                                                                      <  2:4949715/93‑1 (MQ=255)

                    gCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGaa                                                                     >  2:6539301/1‑93 (MQ=255)

                     cGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGCAt                                                                    >  1:6091091/1‑93 (MQ=255)

                       gATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTCTAGAGAGGTCAGGGAAGAAGGTCTTCGAATcc                                                                  >  2:4439562/1‑93 (MQ=255)

                          gCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCAtt                                                               >  2:3205105/1‑93 (MQ=255)

                          gCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCAtt                                                               <  1:6312579/93‑1 (MQ=255)

                           cTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCAttt                                                              <  2:1366319/93‑1 (MQ=255)

                            ttCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTc                                                             <  1:2089770/93‑1 (MQ=255)

                            ttCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTc                                                             >  2:4329802/1‑93 (MQ=255)

                            ttCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTc                                                             <  2:1547581/93‑1 (MQ=255)

                              cTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGa                                                           <  2:182143/93‑1 (MQ=255)

                               tACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGAt                                                          <  1:6033025/93‑1 (MQ=255)

                                 cGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAg                                                        >  2:186083/1‑93 (MQ=255)

                                  gCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGc                                                       >  1:5141960/1‑93 (MQ=255)

                                   cccATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATATCa                                                      >  1:3630342/1‑93 (MQ=255)

                                    ccATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCaa                                                     >  1:873153/1‑93 (MQ=255)

                                    caatatAACTCAAAATGAAGGCAACAAATACTACCAGGCTGACCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCaa                                                     <  1:1157896/91‑1 (MQ=255)

                                     cATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAAc                                                    >  2:3521339/1‑93 (MQ=255)

                                      atatAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGTGAAGAAGGTCTTCGAATCCATTTCGATAGCATCg                                                   >  1:979619/1‑93 (MQ=255)

                                      atatAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACg                                                   >  1:4678925/1‑93 (MQ=255)

                                       tataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACgg                                                  >  2:5947599/1‑93 (MQ=255)

                                        ataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACggg                                                 <  2:3260463/93‑1 (MQ=255)

                                        ataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACggg                                                 <  2:2092594/93‑1 (MQ=255)

                                        ataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACggg                                                 <  1:2257065/93‑1 (MQ=255)

                                        ataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACggg                                                 >  2:1026786/1‑93 (MQ=255)

                                        ataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACggg                                                 <  2:6637680/93‑1 (MQ=255)

                                        ataACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACggg                                                 <  1:121882/93‑1 (MQ=255)

                                           aCTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGc                                              >  2:540664/1‑93 (MQ=255)

                                           aCTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGc                                              >  1:5466191/1‑93 (MQ=255)

                                               aaaaTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCg                                          >  1:6647566/1‑93 (MQ=255)

                                               aaaaTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCg                                          >  2:3370245/1‑93 (MQ=255)

                                               aaaaTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCg                                          >  2:836347/1‑93 (MQ=255)

                                                 aaTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGa                                        <  1:2762440/93‑1 (MQ=255)

                                                 aaTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGa                                        >  2:2326422/1‑93 (MQ=255)

                                                 aaTGAAGGCAACAAATACTACCAGGGTGCCCAATGTAGAGAGGTTAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGa                                        <  2:5733104/93‑1 (MQ=255)

                                                  aTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAg                                       <  2:5247550/93‑1 (MQ=255)

                                                  aTGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAg                                       >  2:709730/1‑93 (MQ=255)

                                                    gAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGACACCGGAGTg                                     >  2:5450665/1‑93 (MQ=255)

                                                      aGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCa                                   >  2:4029872/1‑93 (MQ=255)

                                                      aGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCa                                   >  1:2345413/1‑93 (MQ=255)

                                                      aGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCa                                   >  1:4159888/1‑93 (MQ=255)

                                                        gCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGa                                 <  2:5649239/93‑1 (MQ=255)

                                                        gCAACAAATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGa                                 <  1:256675/93‑1 (MQ=255)

                                                         caacaaATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGAt                                <  2:4590490/93‑1 (MQ=255)

                                                          aacaaATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATa                               >  1:1797651/1‑93 (MQ=255)

                                                          aacaaATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATa                               >  1:6417846/1‑93 (MQ=255)

                                                           acaaATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATaa                              >  1:3452710/1‑93 (MQ=255)

                                                           acaaATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATaa                              <  2:5324107/93‑1 (MQ=255)

                                                           acaaATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATaa                              >  2:5902742/1‑93 (MQ=255)

                                                            caaATACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATaaa                             <  2:4819692/93‑1 (MQ=255)

                                                              aaTACTACCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATa                           <  1:5208136/93‑1 (MQ=255)

                                                                tactacCAGTGTGCCCACTTTAGAGATGTCAGGGAAGAAGGTCTTCTAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGc                         <  2:220952/93‑1 (MQ=255)

                                                                tactacCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGc                         >  1:4545435/1‑93 (MQ=255)

                                                                tactacCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGc                         <  2:1592025/93‑1 (MQ=255)

                                                                    acCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCt                     >  1:3592027/1‑93 (MQ=255)

                                                                    acCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCt                     >  1:6103552/1‑93 (MQ=255)

                                                                    acCAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCt                     <  1:392415/93‑1 (MQ=255)

                                                                      cAGGGTGCTCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTaa                   <  2:2386127/93‑1 (MQ=255)

                                                                      cAGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTTCAGATAAATAGCCGCTaa                   <  2:2860374/93‑1 (MQ=255)

                                                                       aGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAAt                  <  2:2508703/93‑1 (MQ=255)

                                                                       aGGGTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAAt                  >  1:3871413/1‑93 (MQ=255)

                                                                        gggTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATg                 >  1:4644168/1‑93 (MQ=255)

                                                                        gggTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATg                 >  1:1927695/1‑93 (MQ=255)

                                                                        gggTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATg                 <  2:4702807/93‑1 (MQ=255)

                                                                        gggTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATg                 >  1:4987621/1‑93 (MQ=255)

                                                                        gggTGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATg                 <  2:3171313/93‑1 (MQ=255)

                                                                           tGCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGa              <  2:696209/93‑1 (MQ=255)

                                                                            gCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGAt             >  2:2678838/1‑93 (MQ=255)

                                                                            gCCCACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGAt             >  1:5345402/1‑93 (MQ=255)

                                                                             cccACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATc            <  2:6599959/93‑1 (MQ=255)

                                                                              ccACTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCa           <  1:6542017/93‑1 (MQ=255)

                                                                                aCTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCaaa         <  2:2477659/93‑1 (MQ=255)

                                                                                aCTTTAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCaaa         >  1:5971500/1‑93 (MQ=255)

                                                                                  tttAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaat       >  2:3736237/1‑93 (MQ=255)

                                                                                  tttAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaat       <  1:1602479/93‑1 (MQ=255)

                                                                                   ttAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaata      >  2:5531573/1‑93 (MQ=255)

                                                                                    tAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataa     <  1:566905/93‑1 (MQ=255)

                                                                                    tAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataa     <  2:6184008/93‑1 (MQ=255)

                                                                                    tAGAGAGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataa     >  1:5328867/1‑93 (MQ=255)

                                                                                     agagagGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataaa    >  2:3607689/1‑93 (MQ=255)

                                                                                      gagaGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataaat   <  1:4403239/93‑1 (MQ=255)

                                                                                      gagaGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataaat   <  1:3063791/93‑1 (MQ=255)

                                                                                      gagaGGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAaaataaat   <  1:4624817/93‑1 (MQ=255)

TCATTTCATTGACGTGGGTTGCGGATGCTTCTACGCCCATATAACTCAAAATGAAGGCAACAAATACTACCAGGGTGCCCACTTTAGAGAAGTCAGGGAAGAAGGTCTTCGAATCCATTTCGATAGCAACGGGGGCACCGGAGTGCAGATAAATAGCCGCTAATGCGATCAAAATAAATG  >  NC_000913/1567832‑1568011

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