Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	322,719	0	T	C	58.6%	19.4 / 97.8	70	G80G (GGA→GGG)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base C (22/19);  total (37/33)
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
Rejected as consensus: Frequency below/above cutoff threshold.

TAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGC  >  NZ_CP009273/322585‑322855                                                                                                                                       |                                                                                                                                         tagtagGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGc                                                                                                                                      <  2:572304/139‑1 (MQ=255)   gtagGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCa                                                                                                                                     <  1:68464/138‑1 (MQ=255)    tagGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACt                                                                                                                                   <  1:265262/139‑1 (MQ=255)    tagGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACt                                                                                                                                   <  1:283604/139‑1 (MQ=255)                ttgAGGTAGCTCCAATTCTCCAGACCGTGTTATTACGCCCAGTTATCGATATCCAGCGCATTGCCCCGGATGTAGCACATGCCGTTGATCAGCGAAGATCCACCCAGGCCTTTACCGCGTCCGCACTCCATGCGGCGGt                                                                                                                       <  2:297380/137‑1 (MQ=255)                tCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGt                                                                                                                       <  1:375299/139‑1 (MQ=255)                        gCTCCAGTTCTCCCGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCAt                                                                                                                >  1:501686/1‑138 (MQ=255)                        gCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATa                                                                                                               <  1:237412/139‑1 (MQ=255)                          tCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATaaa                                                                                                             >  1:340648/1‑139 (MQ=255)                          tCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATaaa                                                                                                             >  1:574333/1‑139 (MQ=255)                            cAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAAc                                                                                                            <  1:569080/138‑1 (MQ=255)                                    cAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGt                                                                                                   >  1:444428/1‑139 (MQ=255)                                       aCCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTcc                                                                                                >  1:121172/1‑139 (MQ=255)                                        ccGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCg                                                                                               <  2:347060/139‑1 (MQ=255)                                           ggTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGttt                                                                                            >  1:124215/1‑139 (MQ=255)                                           ggTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAgg                                                                                                    >  2:262511/1‑131 (MQ=255)                                           ggTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAgg                                                                                                    <  1:262511/131‑1 (MQ=255)                                            gTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGttt                                                                                            >  2:518179/1‑138 (MQ=255)                                               cTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCAt                                                                                         <  2:412493/138‑1 (MQ=255)                                                 tGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAgg                                                                                      >  2:163456/1‑139 (MQ=255)                                                  gcgcCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCCTAAACGGTTCAGGTTCCGTTTCATAGGc                                                                                     >  1:529107/1‑139 (MQ=255)                                                   cgcCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGTcc                                                                                    <  1:27577/139‑1 (MQ=255)                                                   cgcCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGcc                                                                                    <  2:483713/139‑1 (MQ=255)                                                   cgcCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGcc                                                                                    <  1:352074/139‑1 (MQ=255)                                                   cgcCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGcc                                                                                    <  1:144051/139‑1 (MQ=255)                                                          ttATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                                             <  2:258427/139‑1 (MQ=255)                                                               gagaTCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGt                                                                        <  2:121172/139‑1 (MQ=255)                                                               gagaTCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGt                                                                        <  2:501686/139‑1 (MQ=255)                                                                    ccAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTAcc                                                                   >  2:382331/1‑139 (MQ=255)                                                                             tGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGGTCCGTTTCCTACGCCCCGGTGTAGCGTTTACCCTGTAGCgg                                                          >  2:412638/1‑139 (MQ=255)                                                                             tGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCgg                                                          >  1:31789/1‑139 (MQ=255)                                                                               ccACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCggg                                                         <  1:518179/138‑1 (MQ=255)                                                                               ccACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCggg                                                         <  2:31789/138‑1 (MQ=255)                                                                               ccACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCggg                                                         <  2:383614/138‑1 (MQ=255)                                                                                         aGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAggg                                               >  1:163293/1‑138 (MQ=255)                                                                                          gCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCa                                             <  2:356006/139‑1 (MQ=255)                                                                                            acaTGCCGTTGATCAGCGCCGATCCACCGAAGCCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGc                                           <  1:157445/139‑1 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAAGGTTCAGGTTCCGTTTCATAGGCGCAGTTGTAGGGTTTAACCTGTAGCGGGAAAGCCAGGGCAGCggg                                        >  2:124373/1‑139 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCCTAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCa                                             >  2:138086/1‑134 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCa                                             <  1:138086/134‑1 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                                        >  1:189577/1‑139 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                                        >  1:221203/1‑139 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                                        >  1:27888/1‑139 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                                        >  2:242883/1‑139 (MQ=255)                                                                                               tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGGTCAGGTACCGTTTCATAGGCCCAGTTTTAGCGTTTACCCTGTAGCGGGGATGCCAGGGCAGCggg                                        >  1:132560/1‑139 (MQ=255)                                                                                                gCCGTTTATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCa                                                  >  1:526872/1‑128 (MQ=255)                                                                                                gCCGTTTATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCa                                                  <  2:526872/128‑1 (MQ=255)                                                                                                   gTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCCTAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATc                                    >  2:449908/1‑139 (MQ=255)                                                                                                    ttGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCt                                   >  1:462293/1‑139 (MQ=255)                                                                                                            cgacgaTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGt                                                                                              <  2:222848/72‑1 (MQ=255)                                                                                                            cgacgaTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGt                                                                                              >  1:222848/1‑72 (MQ=255)                                                                                                            cgacgaTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCgg                           >  2:192559/1‑139 (MQ=255)                                                                                                            cgacgaTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGGTCCGTTTCATAGGCCCCGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCgg                           >  2:231544/1‑139 (MQ=255)                                                                                                              acgaTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGGTCCGTTTCATAGGCCCAGTTGTAGCGTTTAACCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTgggggggg                           >  2:243204/1‑132 (MQ=255)                                                                                                                 aTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCa                                             >  2:33874/1‑116 (MQ=255)                                                                                                                 aTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCa                                             <  1:33874/116‑1 (MQ=255)                                                                                                                 aTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGGTGTAGCGTTTACCCTGTAGCGGGGATGCCAGGGCAGCGGGCATCTGGGGGCGGAAgg                       >  2:386481/1‑137 (MQ=255)                                                                                                                  tCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCa                     <  2:310865/139‑1 (MQ=255)                                                                                                                      cccAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGGATGCCAGGGCAGCGGGCATCTGGGTGCGGGAGTCAAAGc                 >  2:168992/1‑139 (MQ=255)                                                                                                                      cccAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTAcc                                                                   <  1:64497/89‑1 (MQ=255)                                                                                                                      cccAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTAcc                                                                   >  2:64497/1‑89 (MQ=255)                                                                                                                       ccAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCg                <  2:574333/139‑1 (MQ=255)                                                                                                                         aGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGGTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGGAGTCAAAGCGa               >  2:439301/1‑138 (MQ=255)                                                                                                                           aCCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCCGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGGGCGGAAGTCAAAGCGATAg            >  2:116205/1‑139 (MQ=255)                                                                                                                           aCCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAg            >  1:388360/1‑139 (MQ=255)                                                                                                                           aCCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAg            >  1:46195/1‑139 (MQ=255)                                                                                                                                     cgTCCGCGCTCCATGCGGCGGTTATTCATAAACGATTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGgccgc  <  2:159269/139‑1 (MQ=255)                                                                                                                                     cgTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTcc         <  1:471227/132‑1 (MQ=255)                                                                                                                                     cgTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGgccgc  <  2:340648/139‑1 (MQ=255)                                                                                                                                     cgTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGGTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTcc         >  2:471227/1‑132 (MQ=255)                                                                                                                                       |                                                                                                                                         TAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGC  >  NZ_CP009273/322585‑322855

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