BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,215	(G)_8→9	noncoding (72/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,207	1	.	G	92.9%	360.6 / 20.7	112	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (2/6); new base G (49/55); total (51/61)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

GTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAAT  >  NC_000913/4606109‑4606300

gTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑g                                                                                               <  2:2911161/100‑1 (MQ=255)

   cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑gggg                                                                                            <  1:1841271/100‑1 (MQ=255)

        ggACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATGTTACTTTTTTTGATATCGCCTCGTGTGGATC‑G‑C‑‑ggtgg‑ggg                                                                                       >  2:528273/1‑99 (MQ=21)

        ggACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑ggggg‑gggg                                                                                      >  1:1584756/1‑99 (MQ=255)

           cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑ggggg‑gggggac                                                                                   <  2:2478934/100‑5 (MQ=255)

           cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑ggggg‑ggggact                                                                                   >  1:318299/1‑96 (MQ=255)

             cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTg                                                                                  >  2:1586366/1‑100 (MQ=255)

             cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTg                                                                                  >  1:3235282/1‑100 (MQ=255)

             cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTg                                                                                  >  1:3050193/1‑100 (MQ=255)

             cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTg                                                                                  <  2:2192009/100‑1 (MQ=255)

             cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑ggggg‑ggggccttg                                                                                 >  1:2080017/1‑94 (MQ=255)

              tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCACCGCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTg                                                                                  <  2:508728/99‑1 (MQ=25)

               aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGaa                                                                                <  1:2821826/100‑1 (MQ=255)

               aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGaa                                                                                >  1:1058215/1‑100 (MQ=255)

                 cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAATCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAAcc                                                                              <  1:3354097/100‑1 (MQ=255)

                  acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGGCATCTTACTTTCTTTGCTCTCGGCTGGTTTGTTGC‑G‑atgt                                                                                               >  2:857567/1‑81 (MQ=18)

                  acCTGAAGCCAGCGAGTCTACCAGATCCGCCACCTTCGCACAGTCATCTTACGTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAAcc                                                                              <  2:1159933/99‑1 (MQ=37)

                   ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AGGGGGGG‑GGGACTTGAAccc                                                                             <  1:290629/100‑1 (MQ=255)

                    cTGAAGCTAGCGCGTCTACCAATTCCGCCACGTTCGCAAGAGCGTCTTACTGGTTTTGAGAGTGCCTCTGTTTGTGT‑G‑A‑‑gggtg‑ggcgccgt                                                                                  >  2:1829357/1‑86 (MQ=2)

                    cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCa                                                                          <  1:3295323/100‑1 (MQ=255)

                     tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCa                                                                          <  1:1939856/100‑1 (MQ=255)

                     tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCa                                                                          <  1:2528267/100‑1 (MQ=255)

                      gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCAc                                                                         <  1:1546566/100‑1 (MQ=255)

                      gAAGCTAGCGCGTCAACCAATTCCGCCACCTTCGCACAGGCAACTTACTTTTCTTGATATCGCCTCGTTTGGTGC‑G‑GG‑GGGGG‑GGGACTTGAACCCCCa                                                                          <  2:1888157/99‑1 (MQ=37)

                       aaccTAGCGCGTCTACCAATTCCGCCACCTCCGCAAACTCATCTTACTTTTTCTGATATCGCCCCCTTTGCTCC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACg                                                                        <  2:3229072/97‑1 (MQ=25)

                       aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTGTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACg                                                                        <  1:969211/100‑1 (MQ=255)

                        aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACg                                                                        <  2:430003/99‑1 (MQ=255)

                        aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGt                                                                       >  1:1098720/1‑100 (MQ=255)

                         gCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTc                                                                      <  1:1586369/100‑1 (MQ=255)

                         gCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCTCACGTcc                                                                     <  1:1787196/100‑1 (MQ=255)

                          cttccgatcttACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTcc                                                                     <  1:2736232/90‑1 (MQ=38)

                          cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGCCTTGAACCCCCACGTcc                                                                     >  1:1946115/1‑100 (MQ=255)

                           tAGCGCGTCTACCAATTCCGCCACCTTCGCACACTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCg                                                                    <  2:577113/100‑1 (MQ=255)

                                cgTCTATCAATTCCGCCACCTTCGCACAGGCACCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTaa                                                                 <  2:674398/98‑1 (MQ=37)

                                cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGa                                                               <  1:1500589/100‑1 (MQ=255)

                                cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGa                                                               >  2:2367599/1‑100 (MQ=255)

                                cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGa                                                               >  2:2366441/1‑100 (MQ=255)

                                 gTCTACCAATTCCGCCACCTACGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGCCCGTAAGa                                                               <  1:809199/99‑1 (MQ=38)

                                  tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAac                                                             <  2:3050181/100‑1 (MQ=255)

                                  tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAac                                                             >  1:787525/1‑100 (MQ=255)

                                   cTACCACTTCCGCAACCTTCTCACAGTCAACTTACGTGTAATGAAATCAACGAATTTGGTGC‑GAA‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAac                                                             <  2:778212/99‑1 (MQ=12)

                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTGCTTTTTTTGACATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAac                                                             <  1:866816/98‑1 (MQ=37)

                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTccagatcggaag                                                           >  2:2736222/1‑90 (MQ=38)

                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAacac                                                           <  2:318298/100‑1 (MQ=255)

                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAacac                                                           >  2:2735428/1‑100 (MQ=255)

                                     aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GCG‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGACCACt                                                          >  2:2311126/1‑100 (MQ=38)

                                     aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACt                                                          >  2:830046/1‑100 (MQ=255)

                                     aCCAATTCCGCCACCTTCGCACAGTCATCTTAC‑TTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                         >  2:3262127/1‑100 (MQ=37)

                                        aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAac                                                       <  1:2898973/100‑1 (MQ=255)

                                        aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAac                                                       >  1:793513/1‑100 (MQ=255)

                                          ttCCGCCACCTTCGCCCAGTCATCTTACTTTCTTTGATATCGCCTCGGTTGGGGG‑G‑GG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                                                      <  1:2262939/99‑1 (MQ=25)

                                           tCCGCCACCTTCGCACAGTCATCTTACTTCTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                    <  1:1992348/100‑1 (MQ=255)

                                            ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                   >  2:1330269/1‑100 (MQ=255)

                                            ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                   >  2:2250762/1‑100 (MQ=255)

                                              gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                 <  1:3089713/100‑1 (MQ=255)

                                              gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                 >  2:2947750/1‑100 (MQ=255)

                                               ccgatctCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                <  2:8909/94‑1 (MQ=25)

                                                 aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                              >  1:2525742/1‑100 (MQ=255)

                                                  ccctCGAACTGTGATCTTAGTTTTTTTGATCAACCCTGGTTTGGTGCGG‑G‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                                 <  2:3235269/93‑1 (MQ=2)

                                                  ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                             <  1:1090363/100‑1 (MQ=255)

                                                     tctcacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                          <  2:107665/97‑1 (MQ=255)

                                                     tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCGCGTCCGTAAGAACACTAGCACCTGAAGCTAgc                                          >  2:2319586/1‑100 (MQ=25)

                                                     tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                           <  2:1098718/99‑1 (MQ=255)

                                                      cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaaagatcgg                                         >  1:8909/1‑93 (MQ=25)

                                                       gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGGCCGGAACAACCCTAACACCTGACGCTAgc                                          >  1:877044/1‑98 (MQ=25)

                                                        cacaGTCGTCTTACTCTGTTTCATATCCCCTCCCTTCGTGA‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                         <  2:196873/98‑1 (MQ=14)

                                                        cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcaga                                       >  1:107664/1‑97 (MQ=255)

                                                        cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                          >  2:1809857/1‑97 (MQ=255)

                                                        cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                       >  1:392590/1‑100 (MQ=255)

                                                        cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                       >  1:2676456/1‑100 (MQ=255)

                                                          caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                       >  2:3047836/1‑98 (MQ=255)

                                                              cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                 <  1:1068731/100‑1 (MQ=255)

                                                              cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                 >  1:2140026/1‑100 (MQ=255)

                                                               aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                <  2:262763/100‑1 (MQ=255)

                                                                tCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                               >  2:2000565/1‑100 (MQ=255)

                                                                 cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                              >  2:3334670/1‑100 (MQ=255)

                                                                   tACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                            >  2:124542/1‑100 (MQ=255)

                                                                     cTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTCAGAACACTAACACAAGAAGCTAGCGCGTCTACCAATTcc                           >  2:89573/1‑99 (MQ=18)

                                                                       ttttttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                        <  1:528274/100‑1 (MQ=255)

                                                                        tttttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGTACCCCCACGTCCGTAAGACCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                       <  1:1829361/100‑1 (MQ=17)

                                                                           ttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                    <  1:2366448/100‑1 (MQ=255)

                                                                           ttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                    <  1:2367606/100‑1 (MQ=255)

                                                                           ttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                    >  1:3235494/1‑100 (MQ=255)

                                                                             gATATCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                  >  1:1915129/1‑100 (MQ=255)

                                                                               tatCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGCCTTAACCCCCCACGTCCGTCAGCACACTAGGACAGAAAGCGAGCCCGTGTACAACTTCCTCCACCTTCg                 >  2:1804317/1‑99 (MQ=9)

                                                                               tatCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                <  1:760004/100‑1 (MQ=255)

                                                                                atCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc               <  1:1889752/100‑1 (MQ=255)

                                                                                atCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc               <  1:1889645/100‑1 (MQ=255)

                                                                                atCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc               >  2:140254/1‑100 (MQ=255)

                                                                                 tCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc              <  1:2795662/100‑1 (MQ=255)

                                                                                 tCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc              <  2:787525/100‑1 (MQ=255)

                                                                                 tCGCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATGCCGCCACCTTGGccc              >  2:653748/1‑100 (MQ=14)

                                                                                   gCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg            >  1:2113165/1‑100 (MQ=255)

                                                                                   gCCTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg            >  1:2911172/1‑100 (MQ=255)

                                                                                    ccTCGTTTGGTGC‑G‑A‑‑GGGGGCGCGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg            <  1:846954/99‑1 (MQ=255)

                                                                                     cTCGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg            <  1:857567/98‑1 (MQ=255)

                                                                                       cGTTTGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg        <  1:2207213/100‑1 (MQ=255)

                                                                                        gTTTGGTGC‑G‑AG‑GGGGG‑GGGCCTTGGACCCCCACGTCCGTGAGAAGGGTAAGAGCTTGGGCTAGCGTGTCTACCGATTGCGCCACCTTCGCCCAGTGCGa       >  2:2720186/1‑100 (MQ=9)

                                                                                         tttGGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg      <  1:2947762/100‑1 (MQ=255)

                                                                                             gTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAGGGAGACGAGCACGTGAAGCTAGCGCGTCTACGAATATCGCCGCCGTCGCCCAGTGCGAGCAat  >  2:1235204/1‑100 (MQ=9)

GTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAAT  >  NC_000913/4606109‑4606300

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