D3A :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,230,311	C→T	G107E (GGA→GAA)	SAUSA300_RS11385 ←	serine hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,230,311	0	C	T	99.2%	1389.3 / ‑5.0	496	G107E (GGA→GAA)	SAUSA300_RS11385	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand): ref base C (0/4); new base T (6/486); total (6/490)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 498 total aligned reads displayed.

TCAACGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTCCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCTTTAGCACGATCA  >  NC_007793/2230212‑2230406

tCAACGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCGTTc                                                                                                 <  1:5526803/100‑1 (MQ=255)

  aaCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCt                                                                                               <  1:230623/100‑1 (MQ=255)

  aaCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCt                                                                                               <  1:2067556/100‑1 (MQ=255)

  aaCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCt                                                                                               <  1:1977861/100‑1 (MQ=255)

  aaCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCt                                                                                               <  1:2861233/100‑1 (MQ=255)

   aCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:1067943/100‑1 (MQ=255)

   aCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:1231726/100‑1 (MQ=255)

   aCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:2659511/100‑1 (MQ=255)

   aCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:5839054/100‑1 (MQ=255)

   aCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:3745642/100‑1 (MQ=255)

   aCGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:3538765/100‑1 (MQ=255)

   aCGAAATTATAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTa                                                                                              <  1:3079673/100‑1 (MQ=255)

    cGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTag                                                                                             <  1:6312214/100‑2 (MQ=255)

    cGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTaa                                                                                             <  1:1302227/100‑1 (MQ=255)

    cGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTaa                                                                                             <  1:6464993/100‑1 (MQ=255)

      aaaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTaacc                                                                                           <  1:6521142/100‑3 (MQ=255)

      aaaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTaaca                                                                                           <  1:2615225/100‑3 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCCCCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:4701621/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:1003026/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:2975222/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:3223511/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:3964700/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:4155841/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:2445321/100‑1 (MQ=255)

       aaTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGc                                                                                          <  1:1850865/100‑1 (MQ=255)

        aTTGTAGAATTTACTACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGcc                                                                                         <  1:4883029/100‑2 (MQ=255)

        aTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAGTGCa                                                                                         <  1:4404755/100‑1 (MQ=255)

        aTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCa                                                                                         <  1:1448754/100‑1 (MQ=255)

        aTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCa                                                                                         <  1:5672347/100‑1 (MQ=255)

        aTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCa                                                                                         <  1:340089/100‑1 (MQ=255)

        aTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCa                                                                                         <  1:4273210/100‑1 (MQ=255)

        aTTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCa                                                                                         <  1:2733151/100‑1 (MQ=255)

         ttGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCag                                                                                        <  1:3789892/100‑2 (MQ=255)

         ttGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCag                                                                                        <  1:1638708/100‑2 (MQ=255)

         ttGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCag                                                                                        <  1:627853/100‑2 (MQ=255)

          tGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCaaa                                                                                       <  1:5074379/100‑2 (MQ=255)

          tGTAGAATTTACAACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCaag                                                                                       <  1:6074896/100‑2 (MQ=255)

           gTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACt                                                                                      <  1:5155290/100‑1 (MQ=255)

            ttcatcctACCGCTAAAATTAACAGGCGCTCCATGAGTCAAGTGACCACCATGACTCAAATTCATTCCTAAATCTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:1916065/93‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAGCTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:5987078/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGccggta                                                                                     <  1:5412670/100‑6 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACtt                                                                                     <  1:3632098/100‑2 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:404832/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:4779477/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:3697851/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:6141591/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:2109135/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:5898911/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGCGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:5569347/100‑1 (MQ=255)

            tAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACGCAAATTCATACCTAAATCTGTGTCGCCCATTTCTAATGCAACTa                                                                                     <  1:1189608/100‑1 (MQ=255)

               aaTTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGt                                                                                  <  1:3473659/100‑1 (MQ=255)

                aTTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTa                                                                                 <  1:1726389/100‑1 (MQ=255)

                   tACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTaagg                                                                              <  1:4239915/100‑3 (MQ=255)

                       cttccgatctaCAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGcc                                                                          <  1:5452882/90‑1 (MQ=255)

                       aCTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGcc                                                                          <  1:4448770/100‑1 (MQ=255)

                       aCTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGcc                                                                          <  1:1594078/100‑1 (MQ=255)

                       aCTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGca                                                                          <  1:766398/100‑2 (MQ=255)

                       aCTAAAATCAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCAGACCTAAAACTGTGTCGACCATTTCTAATGCAACTAAGTAAACAGca                                                                          <  1:1134655/100‑2 (MQ=255)

                        cTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATGCCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCa                                                                         <  1:6597302/100‑1 (MQ=255)

                        cTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGccg                                                                         <  1:1523360/100‑2 (MQ=255)

                         tAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCac                                                                        <  1:4122100/100‑2 (MQ=255)

                         tAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCac                                                                        <  1:3332079/100‑2 (MQ=255)

                         tAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCAt                                                                        <  1:2358364/100‑1 (MQ=255)

                          aaaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCagt                                                                       <  1:571521/100‑3 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCGCCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATgg                                                                      <  1:6390748/100‑2 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATgg                                                                      <  1:845574/100‑2 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATga                                                                      <  1:622797/100‑2 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGt                                                                      <  1:4614908/100‑1 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGt                                                                      <  1:4538033/100‑1 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGt                                                                      <  1:2570597/100‑1 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGt                                                                      <  1:4220897/100‑1 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAAGGCAACTAAGTAAACAGCCATga                                                                      <  1:6626005/100‑2 (MQ=255)

                           aaaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCGAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATga                                                                      <  1:6569610/100‑2 (MQ=255)

                            aaTTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGGGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATgag                                                                     <  1:6363703/100‑3 (MQ=255)

                               tAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGc                                                                  <  1:3866514/100‑1 (MQ=255)

                                  cAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTtg                                                               <  1:1362823/100‑1 (MQ=255)

                                   aGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTtgg                                                              <  1:2024832/100‑2 (MQ=255)

                                    ggCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTtgtg                                                             <  1:4659878/100‑1 (MQ=255)

                                    ggCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTtgtg                                                             <  1:3132666/100‑1 (MQ=255)

                                    ggCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTtgtg                                                             <  1:3903640/100‑1 (MQ=255)

                                    ggCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTtgtc                                                             <  1:5016071/100‑2 (MQ=255)

                                     gcgcTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGa                                                            <  1:1604910/100‑1 (MQ=255)

                                          ccATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTg                                                       <  1:5270726/100‑1 (MQ=255)

                                                   aaGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGtt                                              <  1:3840243/100‑1 (MQ=255)

                                                    aGTGACCACCATGATTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTg                                             <  1:4991949/100‑1 (MQ=255)

                                                    aGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTg                                             <  1:408415/100‑1 (MQ=255)

                                                    aGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTg                                             <  1:2472498/100‑1 (MQ=255)

                                                    aGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTg                                             <  1:2220721/100‑1 (MQ=255)

                                                                 aCTCAAATTCATACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGt                                <  1:1791098/100‑1 (MQ=255)

                                                                            tACCTAAAACTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTcccat                     <  1:70427/100‑4 (MQ=255)

                                                                                  aaaCTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGc               <  1:4825244/100‑1 (MQ=255)

                                                                                    aCTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCGGCTCCAAACAATGCtt             <  1:5370777/100‑1 (MQ=255)

                                                                                    aCTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCtt             <  1:3186962/100‑1 (MQ=255)

                                                                                     cTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCttt            <  1:5187192/100‑1 (MQ=255)

                                                                                     cTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCttt            <  1:1036642/100‑1 (MQ=255)

                                                                                     cTGTGTCGCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCttg            <  1:2876841/100‑2 (MQ=255)

                                                                                            gCCCATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCTTTAGCACgg     <  1:429479/100‑2 (MQ=255)

                                                                                              ccATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCTTTAGCACGATc   >  1:2155925/1‑100 (MQ=255)

                                                                                               cATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCTTTAGCACGATCa  <  1:5557258/100‑1 (MQ=255)

                                                                                               cATTTCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCTTTAGCACGATCa  <  1:4334046/100‑1 (MQ=255)

TCAACGAAATTGTAGAATTTACCACTAAAATTAACAGGCGCTCCATGTGTCAAGTGACCACCATGACTCAAATTCATACCTAAAACTGTGTCGCCCATTCCTAATGCAACTAAGTAAACAGCCATGTTCGCTTGTGAACCTGAATGTGGTTGAACATTGACATGTTCAGCTCCAAACAATGCTTTAGCACGATCA  >  NC_007793/2230212‑2230406

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