BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,181,290	A→T	18.4%	E290D (GAA→GAT)	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,181,290	0	A	T	18.4%	477.0 / 89.1	199	E290D (GAA→GAT)	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base A (99/63); new base T (18/19); total (117/82)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

ACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGC  >  NC_000913/3181191‑3181388

aCCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGaa                                                                                                    >  1:3287210/1‑100 (MQ=255)

 ccGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAg                                                                                                   >  2:3539062/1‑100 (MQ=255)

  cGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAga                                                                                                  <  1:2893008/100‑1 (MQ=255)

   ngCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgat                                                                                                 >  1:3343750/2‑100 (MQ=255)

   ggCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgat                                                                                                 <  1:3937111/100‑1 (MQ=255)

    gCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCATCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatg                                                                                                <  1:2030238/100‑1 (MQ=255)

    gCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTAATTTTGCGGAAgatg                                                                                                >  1:4125963/1‑100 (MQ=255)

      gTGGAATAGCATTATCTACAACAAGGCACTTCTACCGCTACTGTGGGATATGTTCCTGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                              <  2:1823475/100‑1 (MQ=255)

      gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                              <  1:1023703/100‑1 (MQ=255)

       tGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatc                                                                                             >  1:1112873/1‑100 (MQ=255)

         gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcat                                                                                           >  1:669347/1‑100 (MQ=255)

            gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcatc                                                                                          <  2:4196423/98‑1 (MQ=255)

             agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGCTGATCATCCGc                                                                                       >  2:878605/1‑100 (MQ=255)

             agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                       <  2:3284563/100‑1 (MQ=255)

                aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                    <  2:2654344/100‑1 (MQ=255)

                  tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                  <  2:244288/100‑3 (MQ=39)

                  tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                  <  1:1737571/100‑3 (MQ=38)

                  tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                  <  1:1208102/100‑1 (MQ=255)

                  tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                  <  1:722605/100‑1 (MQ=255)

                   aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                 <  1:4164305/100‑4 (MQ=38)

                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                              <  1:494879/100‑1 (MQ=255)

                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGAAAATGGaaa                                                                              >  1:1302063/1‑100 (MQ=255)

                       ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                             <  1:514245/100‑8 (MQ=37)

                       ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                              <  2:826069/99‑7 (MQ=37)

                        caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaaa                                                                            >  1:4021266/1‑100 (MQ=255)

                        caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaaa                                                                            >  1:2613309/1‑100 (MQ=255)

                         aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcggg                                                                           <  2:3589354/100‑10 (MQ=37)

                          acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                          <  2:4032433/100‑1 (MQ=255)

                          acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTCATTTTGCGGATGATCATCTTCCGCaaaataagcgggc                                                                          <  2:2241059/100‑11 (MQ=25)

                           caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                         >  2:2325466/1‑100 (MQ=255)

                           caaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                         >  2:682547/1‑100 (MQ=255)

                            aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                        <  2:3420407/100‑1 (MQ=255)

                            aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                        <  2:1951185/100‑1 (MQ=255)

                             aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                       >  1:3258075/1‑100 (MQ=255)

                              ggCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgtg                                                                      >  1:183436/1‑100 (MQ=255)

                                 aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                   <  2:1729790/100‑1 (MQ=255)

                                  cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCACCCGCAAATGGAAAAATATGTGGTTa                                                                  >  2:1310602/1‑100 (MQ=255)

                                     cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                               <  1:217600/100‑1 (MQ=255)

                                      tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                              >  2:3931102/1‑100 (MQ=255)

                                      tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                              >  2:3090814/1‑100 (MQ=255)

                                      tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATAATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                              >  2:3030341/1‑100 (MQ=255)

                                       aCCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCg                                                             <  2:2688152/100‑1 (MQ=255)

                                        ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTACCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                            >  2:4166284/1‑100 (MQ=255)

                                         cGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAt                                                           >  1:395695/1‑100 (MQ=255)

                                          gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                          <  1:1511846/100‑1 (MQ=255)

                                          gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                          <  1:2974350/100‑1 (MQ=255)

                                          gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCGGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                          <  2:2489563/100‑1 (MQ=255)

                                               tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCACATGGAAAAATATGTGGTTAAACCGATCTtctc                                                     >  2:1842711/1‑100 (MQ=255)

                                                gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                    <  2:1207723/100‑1 (MQ=255)

                                                gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                    >  1:3594296/1‑100 (MQ=255)

                                                    gagaTGTTCCCGAATCACCCGAACCTGCTCCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGTAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                <  2:2523523/100‑1 (MQ=255)

                                                        tgtcgccGAAGCACGCGAACCTGCTGCCTGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAAGATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            <  2:3635023/95‑1 (MQ=255)

                                                        tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            <  2:4125934/100‑1 (MQ=255)

                                                        tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            <  1:3617558/100‑1 (MQ=255)

                                                        tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            >  1:542183/1‑100 (MQ=255)

                                                        tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            <  1:3523652/100‑1 (MQ=255)

                                                         gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                           >  1:2560452/1‑100 (MQ=255)

                                                         gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                           <  1:2126285/100‑1 (MQ=255)

                                                           tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGGTCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                         >  1:2419079/1‑100 (MQ=255)

                                                           tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                         <  1:2644817/100‑1 (MQ=255)

                                                           tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                         >  1:63244/1‑100 (MQ=255)

                                                            cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                        >  1:2522162/1‑100 (MQ=255)

                                                             ccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaa                                       >  2:3816130/1‑100 (MQ=255)

                                                                aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                    <  1:3665504/100‑1 (MQ=255)

                                                                    aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                >  2:2905652/1‑100 (MQ=255)

                                                                     cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                               >  1:3392314/1‑100 (MQ=255)

                                                                          aCCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAt                          <  2:2173295/100‑1 (MQ=255)

                                                                            ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                        <  2:1698777/100‑1 (MQ=255)

                                                                            ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                        >  2:620356/1‑100 (MQ=255)

                                                                             tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                       >  1:1218215/1‑100 (MQ=255)

                                                                             tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                       <  1:917212/100‑1 (MQ=255)

                                                                             tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                       >  2:1649250/1‑100 (MQ=255)

                                                                                tgaccGCTTATTTTGCGGATGATGTTCATCCACAAATGGAAAAATATGAGGTTAAACCGATCTTCTCCCGTGATGGCGGAAACGTGTCGATCATTg                        <  2:3602711/93‑1 (MQ=255)

                                                                                 gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                   >  1:1288356/1‑100 (MQ=255)

                                                                                  cccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                  >  1:3733937/1‑100 (MQ=255)

                                                                                  cccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                  <  1:249080/100‑1 (MQ=255)

                                                                                    ttttccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                >  1:1141443/8‑100 (MQ=37)

                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGACAAATATGTGGCTAAACCGCTCTTCTCCCGTGAAGGCGCAACCCTGTCGATCATTGAGAACg                  >  2:1099159/1‑98 (MQ=255)

                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                <  1:2729958/100‑1 (MQ=255)

                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                >  2:3266322/1‑100 (MQ=255)

                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                >  1:3141350/1‑100 (MQ=255)

                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                >  1:3467010/1‑100 (MQ=255)

                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                <  2:2974419/100‑1 (MQ=255)

                                                                                     tttccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCa               >  2:159339/7‑100 (MQ=37)

                                                                                     gCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCa               >  2:113195/1‑100 (MQ=255)

                                                                                      ttccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCACAAGCGTGTCGATCATTGAGAACGCCa               >  2:3151366/6‑99 (MQ=37)

                                                                                       ttATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa             >  1:1472772/1‑100 (MQ=255)

                                                                                        tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa            <  1:2801996/100‑1 (MQ=255)

                                                                                        ccatttGCGGATGATCATCTTCCGCGAATGGAAAAATATGTGGTTAAACCGATCGTCTCCCGTGAAGGCGCAAAGGTGACAATCATTGAGAACGGCaga             >  2:1017302/4‑97 (MQ=37)

                                                                                         catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc           >  1:3844561/3‑100 (MQ=38)

                                                                                         catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAAACGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc           >  2:3796302/3‑100 (MQ=38)

                                                                                         catttGCGGATGAACATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc           >  1:2846758/3‑100 (MQ=255)

                                                                                            ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAt        <  2:1506697/100‑1 (MQ=255)

                                                                                             tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt       <  1:1909930/100‑1 (MQ=255)

                                                                                              gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg      <  1:841851/100‑1 (MQ=255)

                                                                                              gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg      <  2:669346/100‑1 (MQ=255)

                                                                                                ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa    >  2:2742946/1‑100 (MQ=255)

                                                                                                 gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag   <  1:3074243/100‑1 (MQ=255)

                                                                                                 gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag   >  1:2933947/1‑100 (MQ=255)

                                                                                                  aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc  <  2:1190429/100‑1 (MQ=255)

ACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGC  >  NC_000913/3181191‑3181388

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