BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,282,439	A→T	6.8%	intergenic (+457/+83)	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,439	0	A	T	6.8%	546.2 / 3.9	222	intergenic (+457/+83)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base A (107/83); new base T (8/7); total (122/100)
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 226 total aligned reads displayed.

GGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG  >  CP009273/1282290‑1282587

ggCAGATTTTTTGATTGCGGATTCGTTAGAGAATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                                                                      <  2:2287687/150‑1 (MQ=255)

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

    gAATATTTGATTGCGGATTCGTTTGTGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCg                                                                                                                                                  <  2:5441259/150‑1 (MQ=255)

     aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                 <  1:1236474/150‑1 (MQ=255)

       tatTTGATTGCTGATTCGTTTTAGAAGTCCGGTTCTTTTGACAGTGATGGTGGTGTGGGAAGGATTACTCAGCTCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                                               <  2:2325716/150‑1 (MQ=14)

        atTTGATTGGGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTTGTCGAAGc                                                                                                                                              >  2:3722780/1‑150 (MQ=255)

        atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                                              >  1:5588166/1‑150 (MQ=255)

        atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                                              >  2:5382434/1‑150 (MQ=255)

        atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTGAGTCGAAGc                                                                                                                                              >  2:5382091/1‑150 (MQ=255)

        atTAGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                                              <  2:3889515/150‑1 (MQ=255)

             aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                                         >  1:1640189/1‑150 (MQ=12)

               ttcGGATTCGTTTGAGTATTCCGGGTCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                                                       <  2:3265107/148‑1 (MQ=255)

                 cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCGCTCTCGCTGACGCTCGAGCCGACCCTTAGGCGAAGCTTCTCATcc                                                                                                                                     >  2:253004/1‑150 (MQ=14)

                 cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                                     >  2:1062977/1‑150 (MQ=255)

                  ggATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                                                                                                                    <  1:3486667/150‑1 (MQ=255)

                  ggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTAGAGCCCAACCATAGTCGAAGCATATCAtaca                                                                                                                                    >  1:1533583/1‑147 (MQ=14)

                  ggATGCGTTTGAGCATTGCGTGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                                                                                                                    <  1:358892/150‑1 (MQ=14)

                  ggATGCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                                                                                                                    <  1:3064482/150‑1 (MQ=255)

                   gATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                                                                                                                   <  1:1367373/150‑1 (MQ=255)

                    aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCAcatccatc                                                                                                                                  >  1:1661009/1‑150 (MQ=255)

                       cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                                               >  2:1064551/1‑150 (MQ=255)

                        gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCg                                                                                                                                                  <  2:603272/130‑1 (MQ=9)

                        gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCg                                                                                                                              >  1:2797807/1‑150 (MQ=9)

                        gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTCGAGCCGAAACTTAGTAGAAGCTTATCATCCATCCCCg                                                                                                                              >  1:528650/1‑150 (MQ=14)

                        gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCAGCTCATCCAAccccc                                                                                                                              >  1:603439/1‑149 (MQ=14)

                           tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCaa                                                                                                                           >  1:3389508/1‑149 (MQ=255)

                            gagaTTTCGGGGTCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                                          <  2:4890828/150‑1 (MQ=255)

                            gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                                          <  2:2130769/150‑1 (MQ=255)

                                aTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                      >  1:2600125/1‑150 (MQ=255)

                                   ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACCAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCACaa                                                                                                                   >  1:4322071/1‑150 (MQ=255)

                                     ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGATCCGAACCGTAGTCGAAGCATCTAATCCTTCCCCGCATGGGCAGAata                                                                                                                 >  1:2253086/1‑150 (MQ=255)

                                      gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCAGCCCCGCATAGGCacgagta                                                                                                                >  2:5214382/1‑144 (MQ=11)

                                        gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                                  >  2:2668979/1‑146 (MQ=38)

                                         cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACATTTGTCGAAGCTTCTCATCCTACACCGCATGGGCAGAAAAttta                                                                                                             >  2:1824711/1‑149 (MQ=11)

                                                aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                                                      <  1:1901726/150‑1 (MQ=39)

                                                    gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                                  <  2:2341916/150‑1 (MQ=39)

                                                    gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                                  <  1:5427075/150‑1 (MQ=39)

                                                     aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                                                                 <  1:982411/150‑1 (MQ=39)

                                                        gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGCATATTTGATTGCGGagacggtg                                                                                              >  2:4691803/1‑143 (MQ=11)

                                                          ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCCGATCGGCAGAAAGTATGGTTGGGGAttagagtta                                                                                            >  2:2560382/1‑143 (MQ=11)

                                                            tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgtga                                                                                          >  1:4561181/1‑150 (MQ=255)

                                                             gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGagag                                                                                         <  2:641801/150‑1 (MQ=12)

                                                             gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGagag                                                                                         <  1:1612460/150‑1 (MQ=255)

                                                              ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTggaagt                                                                                        >  2:1689369/1‑145 (MQ=255)

                                                                   aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGTAGAGTTCAgg                                                                                   >  2:3653962/1‑150 (MQ=255)

                                                                   aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCACATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAAAGTTCAgg                                                                                   >  2:3653586/1‑150 (MQ=255)

                                                                        tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                                              >  2:3736270/1‑150 (MQ=255)

                                                                        tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                                              <  1:2405006/150‑1 (MQ=255)

                                                                         aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCGAATATTTGATTGCGGATTCGCTTGTGAGTTAAGGGACtgt                                                                             >  1:2831310/1‑148 (MQ=255)

                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                            >  2:3155157/1‑150 (MQ=255)

                                                                               cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                       <  1:5635930/150‑1 (MQ=255)

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

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

                                                                                  tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                    <  1:404250/150‑1 (MQ=255)

                                                                                   gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                                   <  2:5394606/150‑1 (MQ=255)

                                                                                   gcgcgcGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                                   <  2:2895770/150‑1 (MQ=255)

                                                                                     gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCTGGGACTTTTGAGAGTGAtgg                                                                 >  2:4258770/1‑150 (MQ=255)

                                                                                     gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGATTCGCTTGTGAGTTCAGGGCCTTTTGAAAGTGAtgg                                                                 >  2:5306476/1‑150 (MQ=255)

                                                                                      cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAGTGCGGATTCGCTGGGGCGTTCAGGGACTTTTGAAAGTGaggtg                                                                >  2:3002278/1‑146 (MQ=14)

                                                                                      cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATATTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGTGAGTTCAGGGACTTTTTAAAGTTAtggt                                                                >  2:2046076/1‑150 (MQ=14)

                                                                                      cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCGTTGTCGAAGCTTCTCAACCTTCCCCGCAAGGGCCGAATATTTGATTGCGGATTCGGTTGTGAGTTAAGGGACTTTTGAAACTGAtggt                                                                >  1:4927967/1‑150 (MQ=14)

                                                                                         ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                             <  1:1496747/150‑1 (MQ=255)

                                                                                         ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                             <  2:5588166/150‑1 (MQ=255)

                                                                                         ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                             <  1:4491794/150‑1 (MQ=255)

                                                                                         ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCACGAGCCGAACCTTAGTCGAAGCTTCTCATACTTCCCCGCAAGGGCAGAATATTTTATTTCGGATTCGCTTGTTAGTTAAGGGAATTTTGAAAGTGAtggtggt                                                             >  1:3339750/1‑150 (MQ=14)

                                                                                              ccTTCGGGTCGTTGCCTGCGGCAACTCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                        <  1:341257/150‑1 (MQ=255)

                                                                                              ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTATTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                        >  2:1419683/1‑150 (MQ=255)

                                                                                              ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                        <  2:2349972/150‑1 (MQ=255)

                                                                                              ccTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAGTCGCGGGTGCGGTCAGGGACGTTTGAAAGTGATGGTGGTggggg                                                        >  2:2058997/1‑150 (MQ=14)

                                                                                                 tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                                     <  2:2807658/150‑1 (MQ=255)

                                                                                                 tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                                     <  1:4176598/150‑1 (MQ=255)

                                                                                                 tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCCCAAGGGCCGAAAAATTGATTGCGGATTCGCTTGGGAGTTAAGGGACTTTTGAAAGTGTTGGGGGTGGGGGAAg                                                     >  1:4726526/1‑150 (MQ=14)

                                                                                                 tCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATACTTCCCCGCCAGGGCCGAATATTTTATTGCGGATTCGCTTTGTAGTTAAGGGACTTTTTAAATTGATGGTGGTGGGGGAAg                                                     >  1:4726366/1‑150 (MQ=14)

                                                                                                   gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGCAAGGGCAGAATATTTGATTTCGGATTCGCTTGTTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGggagga                                                   >  1:2722947/1‑150 (MQ=14)

                                                                                                     gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATCTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                                 >  1:5001481/1‑150 (MQ=255)

                                                                                                        ggtGCCCGCGGCAACGCTCCCTCGCAGACGCACGAGCCGAACCTTAGTCGAAGCCTCTCATCCTTCCCCGCATGGGCAGAATATCTGATTGCGGATTCGCCTGAGAGTTCCGGGACTTTTGAAAGCGATGGTGGTGGGGGAAGGATTACt                                              <  2:2025076/148‑1 (MQ=14)

                                                                                                            ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGTGTGTGGTGGAAGGATTACTCTgc                                          >  2:4119077/1‑150 (MQ=17)

                                                                                                             cTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                                         <  1:2141482/150‑1 (MQ=255)

                                                                                                               gCGGCAACGCTCGCTCGCTGACGCTCCATCCTAACCTTAGTCGCAGCTTCTCATCCGTCCCCGCATGGCCATATTTTTAGATTTGGTATTCGCTTGAGAGTTCAGTGACTTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCt                                       <  2:3582326/150‑1 (MQ=17)

                                                                                                                  gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTAGATAGTGAGTGAGGAGGGGGGAGGATTACTCAGCGTTgcg                                    >  2:1344786/1‑150 (MQ=17)

                                                                                                                  gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCATCGATgcg                                    >  2:5486056/1‑150 (MQ=255)

                                                                                                                      cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                <  2:1524128/150‑1 (MQ=255)

                                                                                                                      cGCACTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATAGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                <  2:5670225/150‑1 (MQ=255)

                                                                                                                       gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGCtt                               >  1:1936200/1‑150 (MQ=255)

                                                                                                                       gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAAGGACTTTTGAAAGAGATGGTGGTCGGGGAAGGATTACTCAGCGCTGCGCGCtt                               >  1:4425555/1‑150 (MQ=255)

                                                                                                                         tctctcGCTGTCGCTCGAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                             <  2:13278/150‑1 (MQ=255)

                                                                                                                             tcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                         >  1:5192822/1‑150 (MQ=255)

                                                                                                                              cGCTGGCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt                        <  1:3247879/150‑1 (MQ=255)

                                                                                                                                  gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg                    >  1:1992581/1‑150 (MQ=255)

                                                                                                                                    cGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                  >  1:2564819/1‑150 (MQ=255)

                                                                                                                                     gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                 <  2:2170689/150‑1 (MQ=255)

                                                                                                                                          aGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGAGTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGcc            >  2:3675266/1‑150 (MQ=255)

                                                                                                                                           gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGccg           >  2:2269622/1‑149 (MQ=255)

                                                                                                                                           gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt           <  1:4149018/150‑1 (MQ=255)

                                                                                                                                            ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg          <  2:2579681/150‑1 (MQ=255)

                                                                                                                                             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTACGGGTCGTTTccatc         >  1:699967/1‑147 (MQ=255)

                                                                                                                                               aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCgg       >  2:1467801/1‑150 (MQ=255)

                                                                                                                                                aCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTGAGCGCAGAGCGCTACGCACAGAGGGTAGGTGCCCGCGGc      >  2:1309257/1‑150 (MQ=14)

                                                                                                                                                   tgaGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc   <  2:4452434/148‑1 (MQ=255)

                                                                                                                                                    tAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACg  <  1:3736270/150‑1 (MQ=255)

GGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG  >  CP009273/1282290‑1282587

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