BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,113,649	T→G	K330N (AAA→AAC)	glcA ←	glycolate permease GlcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,113,649	0	T	G	99.2%	497.7 / NA	131	K330N (AAA→AAC)	glcA	glycolate permease GlcA
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (48/82); minor base C (0/1); total (48/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTA  >  NZ_CP009273/3113564‑3113733

ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                  >  2:209660/1‑90 (MQ=255)

 cGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaaa                                                                                 >  2:647191/1‑90 (MQ=255)

 cGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaaa                                                                                 >  1:1562606/1‑90 (MQ=255)

  gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAAc                                                                                <  2:1748238/90‑1 (MQ=255)

  gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGAAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAAc                                                                                <  2:2481347/90‑1 (MQ=255)

     ttCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGc                                                                             >  1:2433337/1‑90 (MQ=255)

         acacTTGTTGATGCAAATGAGGGATCTGGAACTTAATCCCCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCa                                                                         <  1:2857816/90‑1 (MQ=255)

         acacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCa                                                                         <  2:570350/90‑1 (MQ=255)

         acacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCa                                                                         <  2:2034724/90‑1 (MQ=255)

          cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  1:2368663/90‑1 (MQ=255)

          cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  2:245363/90‑1 (MQ=255)

          cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  1:2402914/90‑1 (MQ=255)

          cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  2:2533351/90‑1 (MQ=255)

          cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  1:1659542/90‑1 (MQ=255)

           acTTGTTGATGCAAATGAGGGTTCTGGAAATTAATCAACAGGGAATAAAACGCGCCGCCCGGAGCAAACAACGCGTTAAGCGGCTTCATg                                                                       <  1:2551909/90‑1 (MQ=255)

           acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                       <  2:2483130/90‑1 (MQ=255)

           acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                       <  2:2953787/90‑1 (MQ=255)

           acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                       <  2:417910/90‑1 (MQ=255)

            cTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATgg                                                                      >  2:2143332/1‑90 (MQ=255)

              tgttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTc                                                                    >  1:1939551/1‑90 (MQ=255)

              tgttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTc                                                                    >  1:2410271/1‑90 (MQ=255)

               gttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTcc                                                                   >  1:1321056/1‑90 (MQ=255)

               gttgATGCAAATGAGGGATCTGGAAATTAATCACCAGAGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTcc                                                                   <  1:1531549/90‑1 (MQ=255)

                ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                  <  1:485706/90‑1 (MQ=255)

                ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                  <  2:1193057/90‑1 (MQ=255)

                   aTGCAAATGAGGGATCTGGAAATTAATCCCCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                               <  2:1666317/90‑1 (MQ=255)

                   aTGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                               <  2:1932017/90‑1 (MQ=255)

                     gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                             <  2:1939551/90‑1 (MQ=255)

                     gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                             <  1:1382511/90‑1 (MQ=255)

                      cAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  2:2075699/78‑1 (MQ=255)

                       aaaTGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTg                                                           >  2:2308132/1‑90 (MQ=255)

                             gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                     <  1:119166/90‑1 (MQ=255)

                             gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                     <  2:633507/90‑1 (MQ=255)

                             gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                     <  2:1773593/90‑1 (MQ=255)

                             gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                     <  1:1632804/90‑1 (MQ=255)

                              ggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCa                                                    <  2:1124903/90‑1 (MQ=255)

                              ggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCa                                                    >  1:564935/1‑90 (MQ=255)

                               gATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAc                                                   <  2:855538/90‑1 (MQ=255)

                                 tCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                             <  1:52809/78‑1 (MQ=255)

                                  cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGt                                                >  2:2946867/1‑90 (MQ=255)

                                    ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                              <  2:1354956/90‑1 (MQ=255)

                                    ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                              <  2:766410/90‑1 (MQ=255)

                                    ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                              <  1:2929817/90‑1 (MQ=255)

                                    ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                              >  2:3011125/1‑90 (MQ=255)

                                     gAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTaa                                             <  1:2068541/90‑1 (MQ=255)

                                     gAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTaa                                             <  1:1121481/90‑1 (MQ=255)

                                     gAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTaa                                             <  2:1677650/90‑1 (MQ=255)

                                     gAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTaa                                             <  2:816997/90‑1 (MQ=255)

                                      aaaTTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAg                                            >  1:1234465/1‑90 (MQ=255)

                                       aaTTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGa                                           <  1:2946867/90‑1 (MQ=255)

                                        aTTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGAt                                          >  2:1020767/1‑90 (MQ=255)

                                         ttAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGAtt                                         >  2:2089901/1‑90 (MQ=255)

                                           aaTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaa                                       >  1:2283823/1‑90 (MQ=255)

                                           aaTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaa                                       >  1:1998912/1‑90 (MQ=255)

                                           aaTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaa                                       >  2:444102/1‑90 (MQ=255)

                                            aTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                       <  1:1817675/57‑1 (MQ=255)

                                            aTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                       >  2:1817675/1‑57 (MQ=255)

                                              cACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaaaaa                                    >  2:2167604/1‑90 (MQ=255)

                                                     gAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                             >  1:1114158/1‑58 (MQ=255)

                                                     gAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaaa                                      >  1:2481741/1‑81 (MQ=255)

                                                        tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                          <  1:2281786/90‑1 (MQ=255)

                                                        tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                          <  1:771999/90‑1 (MQ=255)

                                                        tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                          <  1:2779205/90‑1 (MQ=255)

                                                        tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                          <  2:676721/90‑1 (MQ=255)

                                                        tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                          <  1:209660/90‑1 (MQ=255)

                                                          aaaCGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGc                        >  1:1261145/1‑90 (MQ=255)

                                                          aaaCGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGc                        >  1:2778742/1‑90 (MQ=255)

                                                          aaaCGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGc                        <  2:2988896/90‑1 (MQ=255)

                                                            aCGCGCCGCCCGGAGCAAATAACGCCTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTc                      <  1:1100696/90‑1 (MQ=255)

                                                             cgcgCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCa                                                                         <  1:176375/38‑1 (MQ=38)

                                                              gcgcCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGa                    <  2:1088937/90‑1 (MQ=255)

                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  2:924830/1‑90 (MQ=255)

                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  1:1662636/1‑90 (MQ=255)

                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  1:203422/1‑90 (MQ=255)

                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  2:549469/1‑90 (MQ=255)

                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  1:2629746/1‑90 (MQ=255)

                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  2:2674359/1‑90 (MQ=255)

                                                                 ccgccCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATg                 >  1:959290/1‑90 (MQ=255)

                                                                    cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  2:1841387/90‑1 (MQ=255)

                                                                    cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  2:2410271/90‑1 (MQ=255)

                                                                    cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  1:444102/90‑1 (MQ=255)

                                                                    cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  1:2089901/90‑1 (MQ=255)

                                                                    cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  1:220579/90‑1 (MQ=255)

                                                                     ccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAttt             <  2:1470612/90‑1 (MQ=255)

                                                                       ggAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGc           <  2:1148211/90‑1 (MQ=255)

                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAATGGTGACCACGCTCGAATGATTTGccc         <  1:2452422/90‑1 (MQ=255)

                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc         <  1:2487766/90‑1 (MQ=255)

                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc         <  1:2115896/90‑1 (MQ=255)

                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc         <  1:1020767/90‑1 (MQ=255)

                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc         <  2:1335123/90‑1 (MQ=255)

                                                                           cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa       <  1:2449586/90‑1 (MQ=255)

                                                                           cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa       <  2:2464043/90‑1 (MQ=255)

                                                                           cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa       <  1:181726/90‑1 (MQ=255)

                                                                             aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa     >  1:2272186/1‑90 (MQ=255)

                                                                             aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa     >  1:2983041/1‑90 (MQ=255)

                                                                             aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa     >  1:1180527/1‑90 (MQ=255)

                                                                              aTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc    <  2:412726/90‑1 (MQ=255)

                                                                              aTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc    <  1:1482699/90‑1 (MQ=255)

                                                                                aaCGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACta  <  1:549469/90‑1 (MQ=255)

                                                                                aaCGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACta  <  2:1261145/90‑1 (MQ=255)

CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTA  >  NZ_CP009273/3113564‑3113733

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