BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,282,421	A→G	11.0%	intergenic (+439/+101)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,282,421	0	A	G	11.0%	399.0 / 29.8	186	intergenic (+439/+101)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base A (69/79); new base G (13/4); total (88/98)
Fisher's exact test for biased strand distribution p-value = 2.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGG  >  CP009273/1282274‑1282569

cATCCGTCCCCGCTTGGGCAGCCTTATTGATTGGGGATTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                    <  1:3673971/150‑1 (MQ=14)

  tCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                  <  2:2497425/150‑1 (MQ=255)

        cccGCTTGGGCAGAATTTTTGATTGCGGATTCGTTTGAGAATTGCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                            <  2:177779/150‑1 (MQ=255)

         ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                           <  2:1527837/150‑1 (MQ=255)

          cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGa                                                                                                                                          >  2:1644740/1‑150 (MQ=255)

          cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                          >  2:1645191/1‑150 (MQ=255)

            cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAAc                                                                                                                                        >  2:1066713/1‑150 (MQ=255)

             gtGGGCAGAATATAAGAGTGCGGATTCGTTTGAGTATTCGGGGGCTTTTGAAAGTGATGGTGGTGGGGGTAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                                       <  2:2499263/149‑1 (MQ=14)

               gggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt                                                                                                                                     <  1:3101512/150‑1 (MQ=255)

               gggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTCGAGGCGAACCtt                                                                                                                                     >  2:1984023/1‑150 (MQ=255)

                ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCttg                                                                                                                                    >  2:4139200/1‑149 (MQ=18)

                   aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTc                                                                                                                                 >  2:4031061/1‑150 (MQ=12)

                   aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                 <  2:810489/150‑1 (MQ=255)

                   aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                 <  1:3212009/150‑1 (MQ=255)

                   aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCTCGCGGACGCTCGAGCCGAACCTTAGTc                                                                                                                                 >  1:3657149/1‑150 (MQ=255)

                   aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACAAAAGCTGAAGCGCGAGGCGAACCGTAGta                                                                                                                                 >  1:423760/1‑149 (MQ=14)

                    gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCg                                                                                                                                >  2:3191491/1‑150 (MQ=12)

                     aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACACGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                               >  2:3330539/1‑150 (MQ=255)

                       tatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGGGGCCAAGCGATCACGCAGACGCACGACACGAACCAAAGTAGaaa                                                                                                                             >  1:746341/1‑149 (MQ=14)

                       tatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                             >  2:2860974/1‑150 (MQ=255)

                       tatTTGATTGCGGATTAGTTTGAGAATTGCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCATCGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAg                                                                                                                             <  1:3913337/150‑1 (MQ=255)

                        atTTGATTGCGTATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                            <  1:3888445/150‑1 (MQ=255)

                         tttGATTGCGGATTCGTTTGAGTATTACTGGTCTTTAGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                                                                                                           <  2:2171733/150‑1 (MQ=14)

                         cttGATTGGGGATTAGTTTGAGTATTCCTGGGTTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                                                                                                           <  2:4171764/149‑1 (MQ=14)

                          ttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCGTAGTCGAAGCtt                                                                                                                          >  2:2241719/1‑150 (MQ=255)

                              ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGATGCGAACCGTAGTCGAAGCTTCTCa                                                                                                                      >  2:3354565/1‑150 (MQ=14)

                                   gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCtt                                                                                                                 >  2:2297430/1‑150 (MQ=255)

                                       cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                             <  1:2172495/150‑1 (MQ=255)

                                       cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACGCTCGAGGCGAACCTTTGTCGAAGCTTCTCATCCTTcccc                                                                                                             >  2:1280036/1‑150 (MQ=255)

                                        gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCACGCTGACGCTCGAGCCGTACCATAGTCGAAGCTTCTCATCCTTACCCg                                                                                                            >  1:1013704/1‑150 (MQ=14)

                                           tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCATCCCCCCaa                                                                                                         >  2:3257513/1‑149 (MQ=255)

                                           tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGGCTTTGCCGGCGGCAACGCACTCGCGCTGACGCTCTCTCCGAACCCTAGTCGACGCTAATCAACCTTCCCCGCag                                                                                                         >  2:229246/1‑149 (MQ=14)

                                             agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                       <  1:842098/150‑1 (MQ=255)

                                             agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCGCGCTGGCGCTCGAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                       >  2:1021429/1‑150 (MQ=255)

                                             agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCTCTCTCGCGGACGCGAGAGCCGGACCGTAGTCGTAGCTTCTAATCCTTGCCCGCAggg                                                                                                       >  1:1802201/1‑150 (MQ=14)

                                              gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                      >  2:489562/1‑150 (MQ=255)

                                              gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                      >  2:2421959/1‑150 (MQ=255)

                                               taTTCCGGGTTTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                     <  1:586539/149‑1 (MQ=255)

                                                aTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCATCTCATCCTTACCCGCATGGGCa                                                                                                    >  2:3754617/1‑150 (MQ=255)

                                                aTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGAATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCa                                                                                                    >  2:2326451/1‑150 (MQ=255)

                                                  tCCGGGGCTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                                  <  1:2749470/150‑1 (MQ=255)

                                                   ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCTACCAGAAGCACGAGCCCAAACATAGTCGAAGCTTCACAAACATACCCCCAAGGGCCGaa                                                                                                 >  1:4034278/1‑150 (MQ=14)

                                                    cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTCGAGCCGAACCGTATTAGAAGATTCTAATCCTTACACGCAATGGCAGAat                                                                                                >  1:473880/1‑150 (MQ=14)

                                                      gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                                              >  2:4098870/1‑150 (MQ=255)

                                                       ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                             <  1:1263860/150‑1 (MQ=38)

                                                        ctttttGAAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                            <  2:1375819/148‑1 (MQ=37)

                                                         tttaaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCATCACATCCGTCCCCGCATGGGCAGAATATTTg                                                                                           >  2:3883699/8‑150 (MQ=11)

                                                          ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGACAGAATAATTGa                                                                                          >  2:3814157/1‑150 (MQ=255)

                                                           tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACACGCTGACCCTCGAGCCGAACCGTAGTCGAAGCTTATCATCCTTCCCCGCATGGGCCGAATAATTGAt                                                                                         >  1:4136680/1‑150 (MQ=14)

                                                              gtaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGc                                                                                      <  1:634815/148‑1 (MQ=39)

                                                                  gcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                  <  1:1180554/148‑1 (MQ=39)

                                                                  gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCATCGCCCTTCGGGTCCTGTCCGGACGCAAAGCCCTCTCGCGAACCCACGCGCCGAACCTTATTCGAAGCATATCATACGTCCCCCCATGCGCAGAATATTTTAATTCGGAt                                                                                  >  2:3866850/1‑150 (MQ=2)

                                                                      tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                 <  1:435184/83‑1 (MQ=33)

                                                                      tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                 >  2:506234/1‑83 (MQ=33)

                                                                      tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                              <  1:2964523/150‑1 (MQ=40)

                                                                       ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                             <  2:2681127/150‑1 (MQ=255)

                                                                        gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                            <  1:2406104/150‑1 (MQ=255)

                                                                        gtggtggGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                            <  1:3872602/150‑1 (MQ=255)

                                                                               gggAAGGATTTCTCTGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCACCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                     <  1:4163372/150‑1 (MQ=255)

                                                                               gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGtt                                                                     >  2:3703355/1‑150 (MQ=255)

                                                                               gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                     <  2:1185690/150‑1 (MQ=255)

                                                                                 gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                   <  1:1578394/150‑1 (MQ=255)

                                                                                     gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGATTAGCTTGAGAGTTCAgggc                                                               >  2:3146450/1‑149 (MQ=255)

                                                                                      aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGCATCTTTGATTGCGGATTCGCGTGAGAGTTCAGGGAc                                                              >  2:204820/1‑150 (MQ=255)

                                                                                      aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAc                                                              <  1:318421/150‑1 (MQ=255)

                                                                                      aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCGGACGCACGCGCCGAACCATGGTCGAAGCGTCTCATCCGTACCCCCCAGGGCCGGATAATTGATTGCGGATTCGGTTGGGAGTTCAGGGAc                                                              >  1:3503436/1‑150 (MQ=11)

                                                                                         aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTGATTGCGGATTGGGGTGTGAGTTCAGGGAcggt                                                           >  2:748044/1‑147 (MQ=14)

                                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                          <  1:3540813/150‑1 (MQ=255)

                                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                          <  2:4064681/150‑1 (MQ=255)

                                                                                             aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                       <  2:1410180/150‑1 (MQ=255)

                                                                                             aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                       <  1:1046592/150‑1 (MQ=255)

                                                                                             aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTGGCGGACGCGCGTGCCGAACCGTAGTCGAAGCTTCTCAACCTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGACGTTTgga                                                       >  1:1681702/1‑148 (MQ=14)

                                                                                                gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                    <  2:3080008/150‑1 (MQ=255)

                                                                                                gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                    <  2:556373/150‑1 (MQ=255)

                                                                                                 cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                                   <  1:436675/150‑1 (MQ=255)

                                                                                                  tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGA                                                  <  2:881578‑M1/150‑1 (MQ=255)

                                                                                                  tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGA                                                  <  1:2666919‑M1/150‑1 (MQ=255)

                                                                                                        cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCAACCTTCCCCGCCTGGGCAGAATATTTTATTGCGGATTCGGTTGGTAGTTCAGGGACGTGTGAAAGGGCTGGTGG                                            >  2:2856830‑M1/1‑150 (MQ=255)

                                                                                                          tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTG                                          <  2:3419296‑M1/150‑1 (MQ=255)

                                                                                                          tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTG                                          <  1:327259‑M1/150‑1 (MQ=255)

                                                                                                            gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                        >  2:3033572‑M1/1‑150 (MQ=255)

                                                                                                            gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                        <  2:1585787‑M1/150‑1 (MQ=255)

                                                                                                               cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGA                                     <  2:3955333‑M1/150‑1 (MQ=255)

                                                                                                               cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGA                                     <  2:3955091‑M1/150‑1 (MQ=255)

                                                                                                                  cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATACTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGG                                  >  2:3059950‑M1/1‑150 (MQ=255)

                                                                                                                  cGGGTCGTTGCCTGCGGCAACCCTCTATCGCTGAAGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCATCCCCCCAAGGGCACAATATTTTATTGAGGAATCTCGTGATAGTTCTGGGGCTGTTGTGTGTGATGGTGGGGGGGGAAgg                                  >  2:1491513/1‑150 (MQ=9)

                                                                                                                      tCGGTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                              <  2:1205981/150‑1 (MQ=255)

                                                                                                                          tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGATTAATCa                          >  1:1758372/1‑150 (MQ=255)

                                                                                                                           gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg                         >  2:3846910/1‑150 (MQ=255)

                                                                                                                            ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATCTTTAATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTagtaagc                        >  2:1880780‑M1/1‑143 (MQ=255)

                                                                                                                            ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc                        <  2:1251882/150‑1 (MQ=255)

                                                                                                                            ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAga                        >  2:2668426/1‑149 (MQ=255)

                                                                                                                             cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAgat                       >  1:742231/1‑148 (MQ=14)

                                                                                                                               gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCt                     <  1:2860974/150‑1 (MQ=255)

                                                                                                                                   cAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGCGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagt                 <  2:423760‑M1/150‑15 (MQ=255)

                                                                                                                                       gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGTGTTCAGGGACTTTTGAAAGTGATTGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt             >  1:4058375/1‑150 (MQ=255)

                                                                                                                                        ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc            >  1:1913577/1‑150 (MQ=255)

                                                                                                                                            ctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTACCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGGTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTagaaccatcgtagaccatgacgg        >  1:1007804‑M1/1‑127 (MQ=255)

                                                                                                                                                 tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCgg   >  1:2497363/1‑150 (MQ=255)

                                                                                                                                                 tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTctg   >  1:917404/1‑148 (MQ=255)

CATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGG  >  CP009273/1282274‑1282569

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