BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,869,012	T→C	I116V (ATT→GTT)	b3691 ←	D‑galactonate transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,869,012	0	T	C	98.4%	607.0 / ‑1.3	191	I116V (ATT→GTT)	b3691	D‑galactonate transport
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (105/83); total (107/84)
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.53e-01

Only 100 of 202 total aligned reads displayed.

TATAAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAA  >  NC_000913/3868922‑3869101

  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCCCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       >  1:732480/1‑93 (MQ=255)

  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       <  1:2101092/93‑1 (MQ=255)

   aaaaaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCgg                                                                                      <  1:4013357/93‑1 (MQ=255)

     aaaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGtt                                                                                    >  1:1186001/1‑93 (MQ=255)

     aaaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGtt                                                                                    >  1:2506425/1‑93 (MQ=255)

      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   <  1:3711362/93‑1 (MQ=255)

       aCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTAt                                                                                  <  1:4012758/93‑1 (MQ=255)

        ccAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATc                                                                                 <  2:4718249/93‑1 (MQ=255)

         cAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcg                                                                                >  1:63071/1‑93 (MQ=255)

         cAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcg                                                                                >  2:43611/1‑93 (MQ=255)

          aaCGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcgc                                                                               >  2:3601743/1‑93 (MQ=255)

              gCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCgcgcgc                                                                           >  1:4450572/1‑93 (MQ=255)

              gCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCgcgcgc                                                                           <  1:1323549/93‑1 (MQ=255)

               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  2:1114713/93‑1 (MQ=255)

                aGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAg                                                                         >  2:1868331/1‑93 (MQ=255)

                aGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAg                                                                         <  2:1412899/93‑1 (MQ=255)

                 gAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGa                                                                        >  1:181019/1‑93 (MQ=255)

                 gAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGa                                                                        >  2:2940972/1‑93 (MQ=255)

                 gAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGa                                                                        >  2:2369686/1‑93 (MQ=255)

                  aaGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGAc                                                                       <  2:4279038/93‑1 (MQ=255)

                    gcgcGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCa                                                                     >  1:2018251/1‑93 (MQ=255)

                    gcgcGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCa                                                                     <  2:3969950/93‑1 (MQ=255)

                     cgcgTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAAGACCGGTTATCGCGCGCAGACCaa                                                                    >  1:374891/1‑93 (MQ=255)

                      gcgTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAAt                                                                   >  2:1557957/1‑93 (MQ=255)

                      gcgTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAAt                                                                   <  1:3351801/93‑1 (MQ=255)

                       cgTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAAtt                                                                  >  2:3998731/1‑93 (MQ=255)

                       cgTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAAtt                                                                  >  2:762767/1‑93 (MQ=255)

                        gTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTa                                                                 >  1:1593703/1‑93 (MQ=255)

                          tCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAAt                                                               >  1:1744216/1‑93 (MQ=255)

                           cATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATg                                                              <  2:1073437/93‑1 (MQ=255)

                            aTGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGa                                                             <  1:3215912/93‑1 (MQ=255)

                             tGTTCCGTGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGGTATCGCGCGCACACCAATTAATGag                                                            >  2:313233/1‑92 (MQ=255)

                              gTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGTAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACa                                                           <  2:3401696/93‑1 (MQ=255)

                              gTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACa                                                           <  1:2868983/93‑1 (MQ=255)

                                tCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACAtt                                                         <  2:1665283/93‑1 (MQ=255)

                                 ccGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTa                                                        <  2:2815377/93‑1 (MQ=255)

                                 ccGGGAACCACCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGCCATTa                                                        <  2:629579/93‑1 (MQ=255)

                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       <  1:2575796/93‑1 (MQ=255)

                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       <  1:1440170/93‑1 (MQ=255)

                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       <  2:245741/93‑1 (MQ=255)

                                   gggAACCAGCTGGTCACGCTCCGGTTATTGGTCGGGAAAGCAGGCGGTTCGAAAACACCGGCTATCGCGAGCACACCAATTAATGATATTaac                                                      >  2:3869201/1‑92 (MQ=255)

                                   gggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAg                                                      >  2:1151603/1‑93 (MQ=255)

                                    ggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGc                                                     <  1:167778/93‑1 (MQ=255)

                                    ggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGc                                                     >  1:4725517/1‑93 (MQ=255)

                                     gAACCAGCTGGTCGCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGcc                                                    <  2:3768172/93‑1 (MQ=255)

                                      aaCCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGccc                                                   >  2:3409371/1‑93 (MQ=255)

                                       aCCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCg                                                  >  2:3343597/1‑93 (MQ=255)

                                        ccAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGt                                                 <  1:919156/93‑1 (MQ=255)

                                         cAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTg                                                <  1:4108603/93‑1 (MQ=255)

                                          aGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTgg                                               >  2:2273818/1‑93 (MQ=255)

                                                  aCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGccc                                       >  2:201132/1‑93 (MQ=255)

                                                  aCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGccc                                       >  1:3840846/1‑93 (MQ=255)

                                                  aCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGcac                                       >  2:139350/1‑91 (MQ=255)

                                                    cATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTg                                     <  1:1236404/93‑1 (MQ=255)

                                                     aTCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTgg                                    >  1:1017572/1‑93 (MQ=255)

                                                      tCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGa                                   >  1:1549527/1‑93 (MQ=255)

                                                       ccGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGaa                                  >  2:931264/1‑93 (MQ=255)

                                                          gTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaa                               >  1:970701/1‑93 (MQ=255)

                                                          gTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaa                               >  1:707920/1‑93 (MQ=255)

                                                          gTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaa                               >  1:4689426/1‑93 (MQ=255)

                                                            tATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAg                             <  1:3385518/93‑1 (MQ=255)

                                                             aTTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGGCCGTGGCAAAGCCCTGGAATAAAGt                            >  2:2063093/1‑93 (MQ=255)

                                                             aTTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGt                            >  2:8534/1‑93 (MQ=255)

                                                              ttGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTc                           >  1:2789624/1‑93 (MQ=255)

                                                              ttGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTc                           >  1:2119040/1‑93 (MQ=255)

                                                                 gTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGcc                        <  1:2831701/93‑1 (MQ=255)

                                                                   cGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCAc                      <  1:2287403/93‑1 (MQ=255)

                                                                   cGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCAc                      <  1:2133614/93‑1 (MQ=255)

                                                                   cGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCAc                      >  1:2891903/1‑93 (MQ=255)

                                                                    gggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACt                     <  2:388533/93‑1 (MQ=255)

                                                                     ggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTg                    >  1:356975/1‑93 (MQ=255)

                                                                     ggAAAGCAGGCGCTTCGAAAACACCGGATATCGCGCGCAGACCAATTAATGACATTAAGCACGTGGCAAAGCCCTGGAATACAGTCGCCACTg                    >  1:4067836/1‑93 (MQ=255)

                                                                      gAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGa                   >  1:568021/1‑93 (MQ=255)

                                                                      gAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGa                   >  2:1194814/1‑93 (MQ=255)

                                                                      gAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGa                   <  2:3147178/93‑1 (MQ=255)

                                                                       aaaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAc                  <  2:4456295/93‑1 (MQ=255)

                                                                        aaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAcc                 <  2:3657999/93‑1 (MQ=255)

                                                                         aGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCa                >  1:4168229/1‑93 (MQ=255)

                                                                         aGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCa                >  2:1861131/1‑93 (MQ=255)

                                                                           cAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGc              >  1:1394219/1‑93 (MQ=255)

                                                                            aGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGcc             >  2:4169991/1‑93 (MQ=255)

                                                                            aGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGcc             >  2:1885797/1‑93 (MQ=255)

                                                                            aGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGcc             <  2:350418/93‑1 (MQ=255)

                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  1:945683/1‑93 (MQ=255)

                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  2:4048650/1‑93 (MQ=255)

                                                                             ggCGCTTCGAAAACAACGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  2:2785817/1‑93 (MQ=255)

                                                                                gcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAat         <  2:421312/93‑1 (MQ=255)

                                                                                gcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAat         >  1:3497079/1‑93 (MQ=255)

                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        <  2:2308440/93‑1 (MQ=255)

                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        >  2:1562474/1‑93 (MQ=255)

                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       >  1:2081564/1‑93 (MQ=255)

                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       >  2:1136623/1‑93 (MQ=255)

                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       <  2:564732/93‑1 (MQ=255)

                                                                                   tCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATc      >  1:2302703/1‑93 (MQ=255)

                                                                                     gAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGc    >  1:1865793/1‑93 (MQ=255)

                                                                                     gAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGc    >  2:2619929/1‑93 (MQ=255)

                                                                                     gAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGc    >  1:2981847/1‑93 (MQ=255)

                                                                                      aaaaCACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCa   >  2:133616/1‑93 (MQ=255)

                                                                                      aaaaCACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCa   >  2:2551148/1‑93 (MQ=255)

                                                                                       aaaCACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCaa  >  1:453564/1‑93 (MQ=255)

TATAAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAA  >  NC_000913/3868922‑3869101

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