BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,800,753	A→C	V606V (GTT→GTG)	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,800,753	0	A	C	96.2%	701.7 / 7.9	210	V606V (GTT→GTG)	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (4/4); new base C (94/108); total (98/112)
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.95e-01

Only 100 of 210 total aligned reads displayed.

GTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATT  >  NC_000913/1800656‑1800848

gTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCAc                                                                                               >  2:5368891/1‑100 (MQ=255)

gTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCAc                                                                                               <  2:198047/100‑1 (MQ=255)

   tAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACttt                                                                                            <  2:2312065/100‑1 (MQ=255)

     aGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGc                                                                                          <  1:4701265/100‑1 (MQ=255)

     aGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACATCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGc                                                                                          >  2:3810534/1‑100 (MQ=255)

      gACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTAGCACGGCGGGTGCGAACGGCCACTTTGcc                                                                                         >  1:5057825/1‑100 (MQ=255)

      gACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGcc                                                                                         <  2:3981832/100‑1 (MQ=255)

      gACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGcc                                                                                         <  1:3378266/100‑1 (MQ=255)

      gACTGCGGCTGCGAATCTCTTGTTGCAGCGTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGcc                                                                                         <  1:5250663/100‑1 (MQ=255)

       aCTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCt                                                                                        <  1:4060202/100‑1 (MQ=255)

       aCTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCt                                                                                        <  2:1648190/100‑1 (MQ=255)

       aCTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCt                                                                                        <  1:122690/100‑1 (MQ=255)

        cTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTg                                                                                       >  2:4852459/1‑100 (MQ=255)

        cTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTg                                                                                       <  1:3882610/100‑1 (MQ=255)

                gCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCAcc                                                                               <  2:4109328/100‑1 (MQ=255)

                    aTCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCtt                                                                           >  2:4257398/1‑100 (MQ=255)

                    aTCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCtt                                                                           >  2:3217011/1‑100 (MQ=255)

                     tctcTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCttt                                                                          <  1:640953/100‑1 (MQ=255)

                      ctctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTa                                                                         <  1:897569/100‑1 (MQ=255)

                        ctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATc                                                                       <  2:857237/100‑1 (MQ=255)

                        ctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATc                                                                       <  2:2658153/100‑1 (MQ=255)

                        ctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATc                                                                       >  2:374688/1‑100 (MQ=255)

                         ttgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATca                                                                      <  2:1326530/100‑1 (MQ=255)

                          tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcac                                                                     >  2:3158981/1‑100 (MQ=255)

                          tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcac                                                                     >  2:781413/1‑100 (MQ=255)

                          tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcac                                                                     <  1:2662617/100‑1 (MQ=255)

                          tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcac                                                                     >  1:925852/1‑100 (MQ=255)

                          tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcac                                                                     <  2:493291/100‑1 (MQ=255)

                          tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTGCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCGCCTCTTTATcac                                                                     <  2:5285689/100‑1 (MQ=255)

                           gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcacc                                                                    >  2:1764137/1‑100 (MQ=255)

                           gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcacc                                                                    <  2:3527793/100‑1 (MQ=255)

                             tgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcaccac                                                                  >  1:3177419/1‑100 (MQ=255)

                             tgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcaccac                                                                  >  1:4988308/1‑100 (MQ=255)

                             tgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcaccac                                                                  <  1:4854884/100‑1 (MQ=255)

                             tgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATcaccac                                                                  >  1:4127284/1‑100 (MQ=255)

                              gCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACcaca                                                                 >  1:1085481/1‑100 (MQ=255)

                               cAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAg                                                                <  2:2676992/100‑1 (MQ=255)

                               cAGCTTCTCGATCACTTCATCTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAg                                                                <  2:3337263/100‑1 (MQ=255)

                                 ncTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGAc                                                              >  1:2767300/2‑100 (MQ=255)

                                 gCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGAc                                                              >  2:2932346/1‑100 (MQ=255)

                                   ttCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCa                                                            >  2:1108676/1‑100 (MQ=255)

                                      tcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCa                                                         <  1:1531774/100‑1 (MQ=255)

                                       cGATTACTGCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                        >  2:4655590/1‑100 (MQ=255)

                                       cGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                        >  2:2432663/1‑100 (MQ=255)

                                       cGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                        <  2:2532212/100‑1 (MQ=255)

                                        naTCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCata                                                       >  1:4415927/2‑100 (MQ=255)

                                        gATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCata                                                       >  2:1173815/1‑100 (MQ=255)

                                        gATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCCTTATCACCACAGACCAGCata                                                       >  1:2579360/1‑100 (MQ=255)

                                          tCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCatata                                                     <  2:4771972/100‑1 (MQ=255)

                                           cACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCatatat                                                    <  1:2242743/100‑1 (MQ=255)

                                           cACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCatatat                                                    <  1:1513948/100‑1 (MQ=255)

                                            tcTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATg                                                   <  1:320870/99‑1 (MQ=255)

                                            aCTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTTCGAACGGCCACTTTTCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATg                                                   >  2:622600/1‑100 (MQ=255)

                                            aCTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATg                                                   >  1:3279759/1‑100 (MQ=255)

                                             ctccATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATgg                                                  <  1:4965237/97‑1 (MQ=255)

                                              ttCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATggg                                                 <  2:1921089/100‑1 (MQ=255)

                                              ttCATTTACGTCCATGCTTCCCAGGGCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATggg                                                 <  2:1478415/100‑1 (MQ=255)

                                               tCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGga                                                >  1:4225455/1‑100 (MQ=255)

                                                cATTTACTTCCATGCTTCCCAGGTCTTTACCACGGCGGTTGCGAACGGCCACTTTCCCTCTTTCCACCTCTTTATCACCACAGACCAGCATATATGGgac                                               <  1:1863056/100‑1 (MQ=255)

                                                cATTTACGTCCATGCTTCCCCGGTCTTTTCCACGGCGGGTGAGAACGGCCACTTTGCCTGATTCCACCTTTTTATCACCACAGACCAGCATTTATGGgac                                               <  2:2171182/100‑1 (MQ=255)

                                                cATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgac                                               <  2:4873187/100‑1 (MQ=255)

                                                   ttACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcg                                            >  1:1176852/1‑100 (MQ=255)

                                                    tACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcga                                           >  2:3686687/1‑100 (MQ=255)

                                                     aCGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgac                                          >  1:4098318/1‑100 (MQ=255)

                                                      cGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgacg                                         >  1:1933695/1‑100 (MQ=255)

                                                        tCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCa                                       >  1:545971/1‑100 (MQ=255)

                                                          cATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCaaa                                     >  1:2297973/1‑100 (MQ=255)

                                                             gCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAgtg                                  <  1:3293231/100‑1 (MQ=255)

                                                             gCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAgtg                                  <  2:3440647/100‑1 (MQ=255)

                                                              cTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAgtgt                                 >  2:2918427/1‑100 (MQ=255)

                                                                tCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGc                               <  1:4014199/100‑1 (MQ=255)

                                                                tCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGc                               <  1:1740155/100‑1 (MQ=255)

                                                                tCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGc                               <  2:1185729/100‑1 (MQ=255)

                                                                tCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGc                               <  1:1199644/100‑1 (MQ=255)

                                                                 cccAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCt                              <  2:5478530/100‑1 (MQ=255)

                                                                 cccAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGCGTGCt                              >  2:2092930/1‑100 (MQ=255)

                                                                  ccAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTc                             >  1:2216952/1‑100 (MQ=255)

                                                                   cAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcg                            >  1:247195/1‑100 (MQ=255)

                                                                   cAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcg                            >  1:3525227/1‑100 (MQ=255)

                                                                   cAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcg                            >  2:1852870/1‑100 (MQ=255)

                                                                     ggTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                          <  2:4158883/100‑1 (MQ=255)

                                                                          ttACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGAttt                     <  2:2003038/100‑1 (MQ=255)

                                                                            aCCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATAAGGGACGCGACGCAAAGTGTGCTCGCGGATTTTa                   >  2:3628816/1‑100 (MQ=255)

                                                                             ccACGGCGGGTGCGAACGGCCACTTTGTCTGATTCCACCTCTTTATCACCACACACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTaa                  >  1:5178238/1‑100 (MQ=255)

                                                                             ccACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGATCAGCATATATGGGACGCGACGCAAAGGGTGCTCGCGAATTTTaa                  >  1:5125199/1‑100 (MQ=255)

                                                                             ccACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTaa                  <  1:1419334/100‑1 (MQ=255)

                                                                             ccACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTaa                  >  1:575829/1‑100 (MQ=255)

                                                                             ccACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTATGCTCGCGGATTTTaa                  >  1:719194/1‑100 (MQ=255)

                                                                              cACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACGGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTaaa                 <  2:1716738/100‑1 (MQ=255)

                                                                               aCGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAg                <  1:3814439/100‑1 (MQ=255)

                                                                               aCGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAg                <  1:4329316/100‑1 (MQ=255)

                                                                                     ggTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtc          <  2:5124518/100‑1 (MQ=255)

                                                                                        gCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttc       <  1:2960159/100‑1 (MQ=255)

                                                                                         cGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttct      <  1:2946141/100‑1 (MQ=255)

                                                                                         cGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttct      <  2:4587646/100‑1 (MQ=255)

                                                                                         cGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttct      <  1:3679382/100‑1 (MQ=255)

                                                                                         cGAAAGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttct      >  2:4677246/1‑100 (MQ=255)

                                                                                          gAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTtctc     <  1:3075375/100‑1 (MQ=255)

                                                                                           aaCGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCa    >  2:3111558/1‑100 (MQ=255)

                                                                                             cGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCAtt  <  2:4749673/100‑1 (MQ=255)

GTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATT  >  NC_000913/1800656‑1800848

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