BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,093	A→G	16.8%	I39V (ATC→GTC)	smrA →	DNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,093	0	A	G	16.8%	398.2 / 69.8	178	I39V (ATC→GTC)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base A (85/63); new base G (19/11); total (104/74)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 179 total aligned reads displayed.

AATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATG  >  NC_000913/1405995‑1406192

attCGCTGTTTCTTAACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTAt                                                                                                    >  1:2836376/3‑100 (MQ=255)

 aTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACACGTAACCAACGTGCGCCGCAGCGTATc                                                                                                   >  1:35550/1‑100 (MQ=255)

  tCGCTCTTTCTTGACGCCATGGAAGATCTCCAGCCGCTGAAACGTGCTACCGAGGTCCACTGGCATCAAACGCGTAACCAACGTGCGCCGCAGCGTATCg                                                                                                  <  2:3112382/100‑1 (MQ=255)

   cGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCATCGTGCGCCGCAGCGTATcca                                                                                                 <  2:35550/100‑3 (MQ=255)

   cGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGa                                                                                                 >  1:2592120/1‑100 (MQ=255)

    gCTGTTTCTTGACGCCGTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                >  1:1889269/1‑100 (MQ=255)

    gCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                <  2:1283526/100‑1 (MQ=255)

    gCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                <  1:886113/100‑1 (MQ=255)

           cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                         <  1:2689834/100‑1 (MQ=255)

           cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATATCCACTGGCATTCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                         >  2:2937403/1‑100 (MQ=255)

            ttGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCa                                                                                        >  1:1566008/1‑100 (MQ=255)

            tctaCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGcn                                                                                        <  2:1935725/97‑2 (MQ=255)

             tGACGCCTGGGAAGATGTCCAGCCGCTGAAGCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  2:2419880/100‑3 (MQ=255)

             tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  1:1818799/100‑3 (MQ=255)

             tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  2:3377206/100‑3 (MQ=255)

             tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  1:2857292/100‑3 (MQ=255)

             tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAg                                                                                       <  2:4175156/100‑1 (MQ=255)

              gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                      <  1:3350468/100‑1 (MQ=255)

               ncGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCaaga                                                                                     >  1:1935728/2‑97 (MQ=255)

               aCGCCATGGAAGCTGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                     >  1:2762289/1‑100 (MQ=255)

               aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGcg                                                                                     <  2:3956118/100‑2 (MQ=255)

               aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                     >  1:1941805/1‑100 (MQ=255)

               aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                     >  1:1136938/1‑100 (MQ=255)

                cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCtt                                                                                    <  1:671675/100‑1 (MQ=255)

                 gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGcgca                                                                                   <  2:1644707/100‑4 (MQ=39)

                 gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGcgca                                                                                   <  2:2881661/100‑4 (MQ=39)

                 gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTg                                                                                   >  2:638605/1‑100 (MQ=255)

                  ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                  >  1:1038059/1‑100 (MQ=255)

                     tGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATaa                                                                               >  1:3432506/1‑100 (MQ=255)

                        aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAttt                                                                            <  1:918833/100‑1 (MQ=255)

                            tGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTc                                                                        >  2:2291362/1‑100 (MQ=255)

                             gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                       <  2:3289715/100‑1 (MQ=255)

                             gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                       >  1:2187126/1‑100 (MQ=255)

                             gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                       <  1:1224573/100‑1 (MQ=255)

                                cAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcacca                                                                    >  2:14685/1‑100 (MQ=255)

                                   ccGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGCAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATCTCCTCACCACCg                                                                 >  1:2154524/1‑100 (MQ=255)

                                    cGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCgg                                                                >  1:2505658/1‑100 (MQ=255)

                                    cGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCgg                                                                <  1:563941/100‑1 (MQ=255)

                                    cGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCgg                                                                >  1:4101787/1‑100 (MQ=255)

                                     gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                               >  2:3046863/1‑100 (MQ=255)

                                            cGTGCTACCGGTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGGTTGATAATTTCCTCACCACCGGATTTCTCg                                                        <  2:2762282/100‑1 (MQ=255)

                                            cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                        <  1:4059447/100‑1 (MQ=255)

                                            cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                        <  2:1136938/100‑1 (MQ=255)

                                            cGTGCTAACCATGTCCACTGGCATCCAACGCGTAAACAACGTGCGCCGCATCGTATCGACACGCTGCAGCTTGCTAATTTCCTCACCACCGGATTTCTCg                                                        <  2:1032726/100‑1 (MQ=255)

                                             gTGCTACTGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGAAACGCTGAACCTTGATAATTTCCTCACCACCGGATATCTCGa                                                       >  2:1822233/1‑100 (MQ=255)

                                             gTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGa                                                       >  1:3646611/1‑100 (MQ=255)

                                              tGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAc                                                      <  2:4101758/100‑1 (MQ=255)

                                               gCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAca                                                     >  1:1469610/1‑100 (MQ=255)

                                                 tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatc                                                   >  1:3296842/1‑100 (MQ=255)

                                                 tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGAGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatc                                                   >  2:1426340/1‑100 (MQ=255)

                                                     gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcc                                               >  2:2170179/1‑100 (MQ=255)

                                                     gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcc                                               <  2:2454982/100‑1 (MQ=255)

                                                      aTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                              <  1:1042902/100‑1 (MQ=255)

                                                      aTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                              <  1:638606/100‑1 (MQ=255)

                                                       tGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCa                                             >  2:1025993/1‑100 (MQ=255)

                                                        gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                            <  2:3874431/100‑1 (MQ=255)

                                                        gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                            >  1:3187861/1‑100 (MQ=255)

                                                        gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                            >  1:1095328/1‑100 (MQ=255)

                                                        gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                            <  1:124987/100‑1 (MQ=255)

                                                         tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                           <  1:3955052/100‑1 (MQ=255)

                                                         tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                           >  2:82846/1‑100 (MQ=255)

                                                          ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                           >  2:542207/1‑99 (MQ=255)

                                                           cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTaa                                         >  2:3907787/1‑100 (MQ=255)

                                                                 cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACAGGCTGCAGCTTGATAATTTCCTCACCACCGGATTTGTCGACATCATCCCACTAAGTCAGc                                   >  1:382133/1‑100 (MQ=255)

                                                                 cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAg                                    >  2:1685930/1‑99 (MQ=255)

                                                                    ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                >  2:1760649/1‑100 (MQ=255)

                                                                    ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                >  2:2096006/1‑100 (MQ=255)

                                                                    ccAACGCGTAACCAACGTGCGCCGCAGCGTAGCGACACGCTGCAGCTTGATAATTTCCTCACCCCCGGATTTCTCGACATCATCCCACTAAGTGAGCCg                                 >  2:3775726/1‑99 (MQ=255)

                                                                     cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCt                               >  1:982968/1‑100 (MQ=255)

                                                                       ncgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTgg                             >  2:3613153/2‑100 (MQ=255)

                                                                         gcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                           <  2:2684494/100‑1 (MQ=255)

                                                                         gcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGATTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                           <  2:1466789/100‑1 (MQ=255)

                                                                          cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGt                          >  1:373343/1‑100 (MQ=255)

                                                                           gTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCTCCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                           <  2:2517630/98‑1 (MQ=255)

                                                                           gTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGtt                         <  2:2979692/100‑1 (MQ=255)

                                                                            tAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGttt                        >  1:2627208/1‑100 (MQ=255)

                                                                              aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcg                      <  1:2319142/100‑1 (MQ=255)

                                                                               ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                     >  2:87743/1‑100 (MQ=255)

                                                                               ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                     <  1:2360145/100‑1 (MQ=255)

                                                                               ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                     >  1:3718568/1‑100 (MQ=255)

                                                                               ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                     >  1:3158763/1‑100 (MQ=255)

                                                                                cAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggc                    <  2:2966379/100‑1 (MQ=255)

                                                                                  aCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggcgg                  <  1:2309658/100‑1 (MQ=255)

                                                                                    gTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGa                <  1:3036166/100‑1 (MQ=255)

                                                                                     caagctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa               >  1:3080666/7‑100 (MQ=38)

                                                                                     caagctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTCTCGGCGGGaa               >  1:3216797/7‑100 (MQ=37)

                                                                                       atctgcAGCGTGTAGATACGCTGCAGCTTGAAAATTTCCTCACCACCGGATTTCTCGACATCATCACACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg             >  1:2822802/5‑100 (MQ=255)

                                                                                       agctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg             >  1:2530891/5‑100 (MQ=38)

                                                                                       agctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg             >  2:3475478/5‑100 (MQ=255)

                                                                                       agctgcAGCGTGTCGATACGCGGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGCAAgg             >  2:2063890/5‑100 (MQ=255)

                                                                                        gctgcAGCGTGTCGATACGCGGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTACGTCAGCCGCTGGAGTTTcggcgg                  >  2:538433/4‑94 (MQ=255)

                                                                                        gccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg            >  2:3594793/1‑100 (MQ=255)

                                                                                         ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt           >  1:3581594/3‑100 (MQ=255)

                                                                                         ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt           >  1:3737414/3‑100 (MQ=255)

                                                                                          tgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt          >  2:2845725/2‑100 (MQ=255)

                                                                                          tgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt          >  1:3603562/2‑100 (MQ=255)

                                                                                          cgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt          >  2:3583060/1‑100 (MQ=255)

                                                                                           gcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTg         >  2:3070687/1‑100 (MQ=255)

                                                                                              gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCaa      >  1:768701/1‑100 (MQ=255)

                                                                                                  aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg  >  2:3142392/1‑100 (MQ=255)

AATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATG  >  NC_000913/1405995‑1406192

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