BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,332,695	G→A	19.7%	A870T (GCG→ACG)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,332,695	0	G	A	19.7%	531.7 / 119.2	225	A870T (GCG→ACG)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand): ref base G (70/110); new base A (19/25); total (90/135)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 231 total aligned reads displayed.

TAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAT  >  CP009273/1332546‑1332843

   cgcgTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGa                                                                                                                                                   <  2:3503344/150‑1 (MQ=255)

     cgTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCacgacg                                                                                                                                                 <  2:3885195/150‑1 (MQ=255)

     cgTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACg                                                                                                                                                 >  2:3382607/1‑150 (MQ=255)

     cgTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACg                                                                                                                                                 >  2:3973164/1‑150 (MQ=255)

       tAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCGACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGa                                                                                                                                               <  2:3417866/150‑1 (MQ=255)

        aaaaCGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATAGACACCACGACGGAg                                                                                                                                              >  1:3502098/1‑150 (MQ=255)

          aaCGTTAGGGCTAACCGGGGAAGAGAAGATTGATAGTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGtt                                                                                                                                            <  1:990186/150‑1 (MQ=255)

              ttAGGGCTAACCGGGGAAGAGAAGATTGATATTGGTGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGAcc                                                                                                                                        >  2:3207342/1‑150 (MQ=255)

               tAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACct                                                                                                                                       >  2:704184/1‑150 (MQ=255)

               tAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACct                                                                                                                                       >  2:3133768/1‑150 (MQ=255)

               tAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACct                                                                                                                                       >  2:4773111/1‑150 (MQ=255)

                aGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACcta                                                                                                                                      <  2:1565346/150‑1 (MQ=255)

                aGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACcta                                                                                                                                      <  2:3453767/150‑1 (MQ=255)

                  ggCTAACCGGGGAAGAGAAGATTGATCTTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACctact                                                                                                                                    <  1:1981213/150‑1 (MQ=255)

                    cTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCtactac                                                                                                                                  <  1:2652577/150‑1 (MQ=255)

                    cTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGCGTTGACCtacaac                                                                                                                                  >  1:3341815/1‑150 (MQ=255)

                     tAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTAcc                                                                                                                                 <  2:4749032/150‑1 (MQ=255)

                      aaCCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCa                                                                                                                                <  2:4102871/150‑1 (MQ=255)

                      aaCCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCa                                                                                                                                >  2:2389834/1‑150 (MQ=255)

                       aCCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGGTGACCTACTACCAg                                                                                                                               >  1:1402489/1‑150 (MQ=255)

                        ccGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGa                                                                                                                              >  1:802584/1‑150 (MQ=255)

                         cGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGaa                                                                                                                             >  1:250195/1‑150 (MQ=255)

                         cGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGaa                                                                                                                             <  1:772145/150‑1 (MQ=255)

                          ggggAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAac                                                                                                                            <  1:2156645/150‑1 (MQ=255)

                          ggggAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAac                                                                                                                            <  1:4610332/150‑1 (MQ=255)

                           gggAAGAGAAGATTGATAGTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGGAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAacg                                                                                                                           <  1:2300775/150‑1 (MQ=255)

                                  gaataTTGATATTGGCTATCTGCAAAACCTACAACCCGGCGCTACTGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACGGCGACGGAGTTTACCTACTACCAGAACGACGGCAt                                                                                                                    <  2:1734809/150‑1 (MQ=255)

                                  gaagaTTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCAt                                                                                                                    <  2:3100011/150‑1 (MQ=255)

                                    agaTTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCAttt                                                                                                                  >  1:2561466/1‑150 (MQ=255)

                                    agaTTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCAttt                                                                                                                  >  1:128470/1‑150 (MQ=255)

                                    agaTTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCAttt                                                                                                                  >  1:1938268/1‑150 (MQ=255)

                                      aTTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCTTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTg                                                                                                                <  1:967502/150‑1 (MQ=255)

                                         gATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCAt                                                                                                             >  2:2507109/1‑150 (MQ=255)

                                          atatTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCAtt                                                                                                            <  2:1673914/150‑1 (MQ=255)

                                          atatTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCAtt                                                                                                            <  2:2896753/150‑1 (MQ=255)

                                               ggCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTc                                                                                                       <  1:3133768/150‑1 (MQ=255)

                                               ggCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTc                                                                                                       <  1:1337704/150‑1 (MQ=255)

                                               ggCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTc                                                                                                       <  1:704184/150‑1 (MQ=255)

                                                gCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGCGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCa                                                                                                      >  1:3357533/1‑150 (MQ=255)

                                                 cGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACCCCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCAt                                                                                                     >  1:576426/1‑150 (MQ=255)

                                                   aTCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTc                                                                                                   <  1:4573024/150‑1 (MQ=255)

                                                     cTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGt                                                                                                 >  1:4653862/1‑150 (MQ=255)

                                                        cAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAat                                                                                              <  2:1689245/150‑1 (MQ=255)

                                                        cAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAat                                                                                              >  1:1079845/1‑150 (MQ=255)

                                                        cAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAat                                                                                              <  1:1543777/150‑1 (MQ=255)

                                                        cAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCAATATGTCATTCGTAat                                                                                              >  1:3169844/1‑150 (MQ=255)

                                                          aaaCCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAatat                                                                                            >  1:200849/1‑150 (MQ=255)

                                                          aaaCCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAatat                                                                                            >  1:3477005/1‑150 (MQ=255)

                                                           aaCCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATg                                                                                           >  1:1662965/1‑150 (MQ=255)

                                                           aaCCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATg                                                                                           <  2:261174/150‑1 (MQ=255)

                                                            aCCTACAACCCGGCTCGACGGTTCCGTTGACGCTTACGCGCGCGGATGGTATCCAGGAAGTCGTACACTGCCGTTTTCGTATCGACACCGCGACGGAGATGACCTACTACCAGAACGAAGGCATTTTGCATTATGTCATTCTTAATATGt                                                                                          >  2:1386608/1‑150 (MQ=255)

                                                             ccTACCACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGtt                                                                                         <  2:1971686/150‑1 (MQ=255)

                                                              cTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTg                                                                                        >  1:1471515/1‑150 (MQ=255)

                                                                aCAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGaa                                                                                      <  2:2702061/150‑1 (MQ=255)

                                                                 cAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAg                                                                                     >  2:1200030/1‑150 (MQ=255)

                                                                 cAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACCGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAg                                                                                     >  2:4237009/1‑150 (MQ=255)

                                                                  aaCCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGt                                                                                    >  2:680507/1‑150 (MQ=255)

                                                                  aaCCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGt                                                                                    >  2:1926387/1‑150 (MQ=255)

                                                                         cgcgACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACt                                                                             >  2:3012073/1‑150 (MQ=255)

                                                                          gcgACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTa                                                                            <  1:2140871/150‑1 (MQ=255)

                                                                           cgACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTAt                                                                           <  2:1592206/150‑1 (MQ=255)

                                                                             aCGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTAttt                                                                         <  2:1897242/150‑1 (MQ=255)

                                                                                gTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATttgct                                                                      <  1:3612987/150‑1 (MQ=255)

                                                                                    cGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGcc                                                                  <  2:2750955/150‑1 (MQ=255)

                                                                                           gCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTTCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTAtt                                                           >  1:2874562/1‑150 (MQ=255)

                                                                                           gCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTAtt                                                           >  2:1524228/1‑150 (MQ=255)

                                                                                              tACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCCCGACGGAGTTGTCCTACAAACAGAACGACGGCATTTTGCATTATTTCATTAGTAATATGTTTTATTAACAACAATTTTCTTTCCGGTTAttttt                                                        >  1:877241/1‑149 (MQ=255)

                                                                                                 gcgcgcgGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCtt                                                     <  1:1894847/150‑1 (MQ=255)

                                                                                                 gcgcgcgGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCtt                                                     <  2:2541057/150‑1 (MQ=255)

                                                                                                  cgcgcgGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTc                                                    >  2:420349/1‑150 (MQ=255)

                                                                                                   gcgcgGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTcc                                                   <  2:4461100/150‑1 (MQ=255)

                                                                                                        gATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGgcaa                                              >  1:2558214/1‑150 (MQ=255)

                                                                                                        gATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTCGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGgcaa                                              <  2:4063060/150‑1 (MQ=255)

                                                                                                           ggTAGCCAGGAAGTCGGACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCAACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGgcaagca                                           <  1:4210590/150‑1 (MQ=255)

                                                                                                                ccAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGa                                      <  2:2460849/150‑1 (MQ=255)

                                                                                                                ccAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGa                                      >  1:3879631/1‑150 (MQ=255)

                                                                                                                  aGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAAt                                    >  2:4842626/1‑150 (MQ=255)

                                                                                                                    gAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTTCTTCCGGCAAGCAAATGAATTa                                  >  2:443147/1‑150 (MQ=255)

                                                                                                                     aaGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTac                                 <  2:3943431/150‑1 (MQ=255)

                                                                                                                      aGTCGTACCCTGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTaca                                >  2:762885/1‑150 (MQ=255)

                                                                                                                        tCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTActca                              <  1:4050920/150‑1 (MQ=255)

                                                                                                                          gTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAAt                            <  1:3990571/150‑1 (MQ=255)

                                                                                                                          gTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAAt                            <  1:3988713/150‑1 (MQ=255)

                                                                                                                             cccTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCa                         >  1:1380592/1‑150 (MQ=255)

                                                                                                                                tGCCGTTGTCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAaga                      <  2:2282549/150‑1 (MQ=255)

                                                                                                                                    gTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGt                  <  2:3079567/150‑1 (MQ=255)

                                                                                                                                    gTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGt                  >  1:2753718/1‑150 (MQ=255)

                                                                                                                                        tCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTAtt              <  2:11567/150‑1 (MQ=255)

                                                                                                                                        tCGTATCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTAtt              <  1:1090543/150‑1 (MQ=255)

                                                                                                                                         cGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTAttt             <  1:1040494/150‑1 (MQ=255)

                                                                                                                                            aTCGACACCACGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGtt          >  1:3574430/1‑150 (MQ=255)

                                                                                                                                             tCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTc         >  1:1005666/1‑150 (MQ=255)

                                                                                                                                               gACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAg       <  1:4773111/150‑1 (MQ=255)

                                                                                                                                               gACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAg       <  1:3453027/150‑1 (MQ=255)

                                                                                                                                               gACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAg       >  2:485454/1‑150 (MQ=255)

                                                                                                                                               gACACCACGACGGTGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAg       <  1:3464787/150‑1 (MQ=255)

                                                                                                                                                 cacCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa     <  1:171444/150‑1 (MQ=255)

                                                                                                                                                 cacCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa     <  1:172313/150‑1 (MQ=255)

                                                                                                                                                    cgcgACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAt  <  2:128351/150‑1 (MQ=255)

                                                                                                                                                    cgcgACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAt  <  2:2558214/150‑1 (MQ=255)

TAACGCGTAAAACGTTAGGGCTAACCGGGGAAGAGAAGATTGATATTGGCGATCTGCAAAACCTACAACCCGGCGCGACGGTTCCGGTGACGCTTACGCGCGCGGATGGTAGCCAGGAAGTCGTACCCTGCCGTTGTCGTATCGACACCGCGACGGAGTTGACCTACTACCAGAACGACGGCATTTTGCATTATGTCATTCGTAATATGTTGAAGTAACAACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAT  >  CP009273/1332546‑1332843

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