BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,923,969	C→G	15.0%	E49Q (GAG→CAG)	yqcC ←	DUF446 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,969	0	C	G	15.0%	324.4 / 46.1	127	E49Q (GAG→CAG)	yqcC	DUF446 family protein
Reads supporting (aligned to +/- strand): ref base C (60/46); new base G (8/11); total (68/59)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 129 total aligned reads displayed.

CATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACG  >  NC_000913/2923871‑2924067

cATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTcc                                                                                                   <  1:1962605/100‑1 (MQ=255)

 aTTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCa                                                                                                  >  1:2037474/1‑100 (MQ=255)

    tCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTg                                                                                               <  1:1674176/100‑1 (MQ=255)

     cATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTgg                                                                                              <  2:112330/100‑1 (MQ=255)

        aataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTc                                                                                           >  1:2571328/1‑100 (MQ=255)

        aataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTc                                                                                           <  1:200164/100‑1 (MQ=255)

         ataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGttgc                                                                                          <  2:1266592/100‑3 (MQ=255)

         ataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGttgc                                                                                          <  2:751525/100‑3 (MQ=255)

          taGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCa                                                                                         <  2:939142/100‑1 (MQ=255)

          taGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCa                                                                                         <  1:1352833/100‑1 (MQ=255)

          taGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCa                                                                                         >  2:1102913/1‑100 (MQ=255)

           aGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCCGCACCCACTGCAACCACTCCAGTGGTTCCAt                                                                                        >  2:1776470/1‑100 (MQ=255)

           aGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCag                                                                                        <  1:640733/100‑2 (MQ=255)

           aGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCa                                                                                         <  1:179547/99‑1 (MQ=255)

           aGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCAt                                                                                        >  2:2101213/1‑100 (MQ=255)

           aGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGAACCCACTGCAACCACTCCAGTGGTTCCAt                                                                                        >  2:1783329/1‑100 (MQ=255)

            gggTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCagt                                                                                       <  1:307498/100‑3 (MQ=255)

            gggTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATg                                                                                       >  2:1630006/1‑100 (MQ=255)

            gggTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATg                                                                                       >  2:1394940/1‑100 (MQ=255)

             ggTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGttgcagtg                                                                                      <  2:123671/100‑7 (MQ=255)

             ggTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATgg                                                                                      <  1:77279/100‑1 (MQ=255)

              gTGCAACGGCAAAAGCGCCCGGTATCGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATgg                                                                                      <  1:1114889/99‑1 (MQ=255)

               tGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGttgcagtggg                                                                                    <  1:1878485/100‑9 (MQ=255)

               tGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGgtg                                                                                    >  2:911675/1‑100 (MQ=255)

                gCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGgtgt                                                                                   <  2:1670346/100‑1 (MQ=255)

                 cAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGttgcagtgggtg                                                                                  <  2:1178067/100‑11 (MQ=38)

                  aaCGGAAAAAGCGCCCGTCTACGGATATTTGTTATCTAACAAATCGTGCATTCGGTGTGTCAGCACCCAGTGCAACCTCTTCAGTGGTTCCATGGTGTct                                                                                 <  1:1133427/100‑2 (MQ=255)

                   aCGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCa                                                                                <  2:763064/100‑1 (MQ=255)

                     ggCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATa                                                                              <  2:563359/100‑1 (MQ=255)

                       cAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCAGaaa                                                                            >  1:2585978/1‑100 (MQ=255)

                         aaaGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGa                                                                          >  1:101456/1‑100 (MQ=255)

                              gcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGt                                                                     >  2:1595220/1‑100 (MQ=255)

                                ccGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTg                                                                   >  2:145568/1‑100 (MQ=255)

                                  ggTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTggg                                                                 >  2:1723851/1‑100 (MQ=255)

                                  ggTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGAGCAGCACCCACTGCAACCGCTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTgg                                                                  <  1:2173527/99‑1 (MQ=255)

                                   gTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACGCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGt                                                                <  2:1557239/100‑1 (MQ=255)

                                     aaCGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTAc                                                              <  1:1288786/100‑1 (MQ=255)

                                       cGGCTGTTTGTTATCAAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCAATGGTGTCCATAAAGAACGGTTGGGTACTa                                                            >  1:2069537/1‑100 (MQ=255)

                                       cGGCTGATTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCAGGGTGTCCATAAAGAACGGTTGGGGACTa                                                            >  1:2071045/1‑100 (MQ=255)

                                        ggCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTAt                                                           >  1:942083/1‑100 (MQ=255)

                                        cacgatttgttAGATAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTAt                                                           >  1:1744325/6‑100 (MQ=255)

                                         gCTGTTTGGTATCTAACAAATCGTGCATTGGCGGGATCAGCACCCACTGCAACCACTGCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTAtt                                                          <  1:1004561/100‑1 (MQ=255)

                                          cTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTa                                                         <  1:1905625/100‑1 (MQ=255)

                                           tgtttgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaa                                                        >  2:1382917/1‑100 (MQ=255)

                                            gtttgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaaa                                                       <  1:3274/100‑1 (MQ=255)

                                               tgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTg                                                    <  1:1459397/100‑1 (MQ=255)

                                               ngttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTCTCCATAAAGAACGGTTGGGTACTATTAAAtt                                                     >  1:2375205/2‑99 (MQ=255)

                                                gttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGa                                                   <  2:590959/100‑1 (MQ=255)

                                                  tATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTAGTTCCATGCTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATg                                                 <  1:1385771/100‑1 (MQ=255)

                                                   aTCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATACAGAACGCTTGGGTACTATTAAATTGATGc                                                >  1:2193524/1‑100 (MQ=255)

                                                     cTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCgg                                              >  2:1163998/1‑100 (MQ=255)

                                                      tAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGc                                             >  2:700351/1‑100 (MQ=255)

                                                       atcAAATCGTGCAGTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCt                                            >  2:1758084/3‑100 (MQ=255)

                                                            aTCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGt                                       <  2:2131178/100‑1 (MQ=255)

                                                             tCGTGCGTTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGtt                                      >  2:2041152/1‑100 (MQ=255)

                                                              cGTGCATTCGCGGGATCAGCGCCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGTGTACTATTAAATTGATGCGGCTGGGGTTc                                     >  1:1580715/1‑100 (MQ=255)

                                                              cGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTc                                     <  1:1163994/100‑1 (MQ=255)

                                                               gTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCa                                    <  2:2176145/100‑1 (MQ=255)

                                                                tGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCAt                                   <  2:1406029/100‑1 (MQ=255)

                                                                  cgttCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATc                                  <  1:2041140/97‑1 (MQ=255)

                                                                   aTTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGt                                >  2:463583/1‑100 (MQ=255)

                                                                      cgcgGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTc                             >  1:1251709/1‑100 (MQ=255)

                                                                       gcgGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGTGGCTGGGGTTCATCGTTTcn                            <  1:427951/100‑2 (MQ=255)

                                                                       gcgGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCg                            >  1:540615/1‑100 (MQ=255)

                                                                         gggATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGcc                          <  2:731121/100‑1 (MQ=255)

                                                                         gggATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGcc                          <  1:1089388/100‑1 (MQ=255)

                                                                           gATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAg                        >  1:1643317/1‑100 (MQ=255)

                                                                           gATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAg                        >  1:1012544/1‑100 (MQ=255)

                                                                               aGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCt                    >  2:1069819/1‑100 (MQ=255)

                                                                               aGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCt                    <  1:1863390/100‑1 (MQ=255)

                                                                                gaaCCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTg                   <  1:1783320/98‑1 (MQ=255)

                                                                                 ctcccACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGa                  <  2:1643326/98‑1 (MQ=255)

                                                                                 cACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGa                  <  1:535479/100‑1 (MQ=255)

                                                                                 cACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGa                  <  2:2280122/100‑1 (MQ=255)

                                                                                  aCCCACTGCAACCACGCCACTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTGTTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGAt                 >  1:1091527/1‑100 (MQ=255)

                                                                                   cccACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATg                >  1:1754274/1‑100 (MQ=255)

                                                                                     ccatggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGtt              >  2:2361102/7‑100 (MQ=255)

                                                                                     ccatggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGtt              >  1:2141435/7‑100 (MQ=255)

                                                                                     ccatggaaCCACTGGAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGtt              >  1:2234087/7‑100 (MQ=255)

                                                                                     ccatggaaCCACTGGAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTCGGGTTCATCGTTTCGCCAGTGCTGATGtt              >  1:1554022/7‑100 (MQ=255)

                                                                                      aCTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTc             >  2:295376/1‑100 (MQ=255)

                                                                                      aCTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTc             >  1:1371692/1‑100 (MQ=255)

                                                                                       atggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCa            >  2:1227714/5‑100 (MQ=255)

                                                                                       atggaaCCACTGGAGTGGTTCCATGGTGTCCATCAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTTCTGATGTTCa            >  1:518524/5‑100 (MQ=255)

                                                                                       atggaaCCACTGGAGTGGTTCCATGGTGTCCATAAAGAACGATTGGGTACTATTAAATTGAGGCGGCTGGGGTTTATCGTTTCGCCAGTGCTGATGTTCa            >  1:1413497/5‑100 (MQ=255)

                                                                                        tGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCAc           >  1:1012210/1‑100 (MQ=255)

                                                                                        tGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCAc           >  1:1205129/1‑100 (MQ=255)

                                                                                         gCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACg          <  2:1922557/100‑1 (MQ=255)

                                                                                          cAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGc         <  1:1394932/100‑1 (MQ=255)

                                                                                          cAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGc         >  1:715608/1‑100 (MQ=255)

                                                                                            aCCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAg       <  2:1012547/100‑1 (MQ=255)

                                                                                            aCCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAg       >  1:46755/1‑100 (MQ=255)

                                                                                             ccACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCGGTGAGGATGTTCACGCa        >  1:1864565/1‑98 (MQ=255)

                                                                                             ccACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGt      <  1:1858675/100‑1 (MQ=255)

                                                                                             ccACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGt      >  1:2218181/1‑100 (MQ=255)

                                                                                             ccACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGt      >  1:605066/1‑100 (MQ=255)

                                                                                             ccACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATGAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAg       <  1:562140/99‑1 (MQ=255)

                                                                                               aCTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTaa    >  1:2276953/1‑100 (MQ=255)

                                                                                                cTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAAc   >  1:2166428/1‑100 (MQ=255)

                                                                                                 tCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACg  <  1:1069815/100‑1 (MQ=255)

CATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACG  >  NC_000913/2923871‑2924067

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