BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,872,672	A→C	5.9%	Q234P (CAG→CCG)	yeaP →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,872,672	0	A	C	5.9%	1058.6 / 50.2	323	Q234P (CAG→CCG)	yeaP	diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base A (169/135); new base C (13/6); total (182/141)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 325 total aligned reads displayed.

GGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAA  >  CP009273/1872525‑1872820

 gcgggcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTAAgg                                                                                                                                                    >  2:2498345/1‑149 (MQ=255)

   gggcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGtt                                                                                                                                                  <  2:2198684/149‑1 (MQ=255)

      cgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGc                                                                                                                                               <  1:3564814/149‑1 (MQ=255)

         tttttGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAg                                                                                                                                            >  2:873098/1‑149 (MQ=255)

           tttGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAgcg                                                                                                                                          >  2:409592/1‑149 (MQ=255)

           tttGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAgcg                                                                                                                                          <  1:1548068/149‑1 (MQ=255)

           tttGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAgcg                                                                                                                                          >  2:2751615/1‑149 (MQ=255)

              gAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATCGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCt                                                                                                                                       >  2:1540917/1‑149 (MQ=255)

              gAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCt                                                                                                                                       <  1:2240158/149‑1 (MQ=255)

              gAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCt                                                                                                                                       <  1:415279/149‑1 (MQ=255)

                  aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAAt                                                                                                                                   <  1:2020183/149‑1 (MQ=255)

                  aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAAt                                                                                                                                   >  2:1285096/1‑149 (MQ=255)

                   tCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATa                                                                                                                                  <  1:7750/149‑1 (MQ=255)

                      gacgacACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGc                                                                                                                               >  2:223326/1‑149 (MQ=255)

                      gacgacACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGc                                                                                                                               >  2:1332605/1‑149 (MQ=255)

                       acgacACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCt                                                                                                                              >  2:3308636/1‑149 (MQ=255)

                          acacTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCa                                                                                                                           >  1:1249168/1‑149 (MQ=255)

                                ttttCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCaaaa                                                                                                                     <  1:148394/149‑1 (MQ=255)

                                      cTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGa                                                                                                               >  1:132639/1‑149 (MQ=255)

                                          cccGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGtt                                                                                                           <  1:400723/149‑1 (MQ=255)

                                            cGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTAt                                                                                                         <  1:2997705/149‑1 (MQ=255)

                                              gCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTg                                                                                                       >  1:3300087/1‑149 (MQ=255)

                                               cATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTgg                                                                                                      >  1:4102105/1‑149 (MQ=255)

                                                aTCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgt                                                                                                     >  2:553246/1‑149 (MQ=255)

                                                aTCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgt                                                                                                     <  2:215610/149‑1 (MQ=255)

                                                  cTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAAAGGCGAAGTTATTGgtcg                                                                                                   >  1:297462/1‑149 (MQ=255)

                                                   ttAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcgt                                                                                                  >  1:1044322/1‑149 (MQ=255)

                                                    tAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcgtc                                                                                                 <  1:4600711/149‑1 (MQ=255)

                                                    tAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcgtc                                                                                                 <  1:698211/149‑1 (MQ=255)

                                                     aaCCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATTATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGtct                                                                                                <  1:755875/149‑1 (MQ=255)

                                                      aCCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGtctc                                                                                               >  1:4681210/1‑149 (MQ=255)

                                                      aCCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGtctc                                                                                               <  1:4394656/149‑1 (MQ=255)

                                                        cATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTcgg                                                                                             >  2:331809/1‑149 (MQ=255)

                                                          tAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTcggcg                                                                                           >  2:3703680/1‑149 (MQ=255)

                                                          tAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTcggcg                                                                                           >  1:261040/1‑149 (MQ=255)

                                                            aGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGt                                                                                         <  2:4448178/149‑1 (MQ=255)

                                                             gATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtg                                                                                        <  1:2511537/149‑1 (MQ=255)

                                                              aTAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtga                                                                                       <  1:2210677/149‑1 (MQ=255)

                                                                aaTGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtgatg                                                                                     >  2:1354136/1‑149 (MQ=255)

                                                                 aTGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATACTCGTGGCGATCTGTTTCTCATTCCGGTTGGCGAGCGCCTTAATACGCTCCAGCACAATGGCGAAGTTAGTGGTCGTCTCGGCGGtgatga                                                                                    >  2:1057558/1‑149 (MQ=255)

                                                                  tGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAg                                                                                   >  1:2410331/1‑149 (MQ=255)

                                                                  tGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAg                                                                                   <  1:2872098/149‑1 (MQ=255)

                                                                  tGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAg                                                                                   <  1:4041638/149‑1 (MQ=255)

                                                                    aTCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGtt                                                                                 <  2:3171350/149‑1 (MQ=255)

                                                                      cgcgTTTATCGATCTGGATAACTTCAGATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGtttt                                                                               >  2:1878069/1‑149 (MQ=255)

                                                                        cgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTa                                                                             >  2:3667151/1‑149 (MQ=255)

                                                                         gTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAg                                                                            <  2:2637497/149‑1 (MQ=255)

                                                                          tttATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAgt                                                                           <  1:1901523/149‑1 (MQ=255)

                                                                           ttATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAgtt                                                                          >  1:1707811/1‑149 (MQ=255)

                                                                            tATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAgttg                                                                         >  2:1987889/1‑149 (MQ=255)

                                                                              tCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAgttgtt                                                                       <  1:3780924/149‑1 (MQ=255)

                                                                                gATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTc                                                                     >  1:4448402/1‑149 (MQ=255)

                                                                                 aTCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCa                                                                    >  2:3145131/1‑149 (MQ=255)

                                                                                      gATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTaaa                                                               >  1:2168022/1‑149 (MQ=255)

                                                                                      gATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTaaa                                                               >  1:3121536/1‑149 (MQ=255)

                                                                                      gATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTaaa                                                               <  2:769908/149‑1 (MQ=255)

                                                                                      gATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTaaa                                                               <  1:1644947/149‑1 (MQ=255)

                                                                                       aTAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAaac                                                              <  1:3471984/149‑1 (MQ=255)

                                                                                        tAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAaaca                                                             <  1:3732906/149‑1 (MQ=255)

                                                                                         aaCTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAaacaa                                                            >  2:1100602/1‑149 (MQ=255)

                                                                                          aCTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAaacaac                                                           >  2:4408176/1‑149 (MQ=255)

                                                                                           cTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAaacaaag                                                          >  1:4309705/1‑147 (MQ=255)

                                                                                           cTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACg                                                          >  2:3087851/1‑149 (MQ=255)

                                                                                            ttCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACga                                                         <  1:991917/149‑1 (MQ=255)

                                                                                            ttCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACga                                                         >  2:3376570/1‑149 (MQ=255)

                                                                                             tCAAATTACTCAATGATCGTTTTGGTCATAAGAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACgag                                                        <  2:542427/149‑1 (MQ=255)

                                                                                                   tAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCg                                                  >  1:247328/1‑149 (MQ=255)

                                                                                                    aaTCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCgg                                                 <  1:54247/149‑1 (MQ=255)

                                                                                                     aTCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGa                                                >  1:3201622/1‑149 (MQ=255)

                                                                                                       cAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGata                                              <  1:2750351/149‑1 (MQ=255)

                                                                                                        actGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGatat                                             <  2:4272882/147‑1 (MQ=255)

                                                                                                         aTGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATAtt                                            <  1:2283621/149‑1 (MQ=255)

                                                                                                         aTGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATAtt                                            <  1:2533474/149‑1 (MQ=255)

                                                                                                          tGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATAttt                                           <  1:3244336/149‑1 (MQ=255)

                                                                                                          tGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATAttt                                           >  1:2705556/1‑149 (MQ=255)

                                                                                                          tGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATAttt                                           <  2:1667126/149‑1 (MQ=255)

                                                                                                             tCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTtcgt                                        <  1:4187043/149‑1 (MQ=255)

                                                                                                                ttttGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGc                                     <  1:2905075/149‑1 (MQ=255)

                                                                                                                ttttGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGc                                     >  1:1417497/1‑149 (MQ=255)

                                                                                                                   tGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGc                                  >  1:3973611/1‑149 (MQ=255)

                                                                                                                   tGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGGTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGACTGCGGCTATTTCGTCGCTGc                                  >  2:2217109/1‑149 (MQ=255)

                                                                                                                      tCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCgag                               >  1:1702973/1‑149 (MQ=255)

                                                                                                                       cATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCgaga                              <  2:680654/149‑1 (MQ=255)

                                                                                                                       cATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCgaga                              >  2:414077/1‑149 (MQ=255)

                                                                                                                        ataataGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGaa                             <  1:3660224/149‑1 (MQ=255)

                                                                                                                          aataGTGGCGATCTGTTTCTCATTCCGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACg                           >  1:4358221/1‑149 (MQ=255)

                                                                                                                           ataGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGc                          >  1:2096947/1‑149 (MQ=255)

                                                                                                                              gTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCAtt                       <  1:2374408/149‑1 (MQ=255)

                                                                                                                                ggCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTca                     >  2:279193/1‑149 (MQ=255)

                                                                                                                                  cGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTcagc                   <  2:3055592/149‑1 (MQ=255)

                                                                                                                                   gATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTcagca                  >  1:111634/1‑149 (MQ=255)

                                                                                                                                   gATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTcagca                  <  2:400826/149‑1 (MQ=255)

                                                                                                                                    aTCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTcagcag                 >  2:2609998/1‑149 (MQ=255)

                                                                                                                                       tGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCaa              >  1:1375801/1‑149 (MQ=255)

                                                                                                                                         tttCTCATTCCGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAAt            >  1:4201233/1‑149 (MQ=255)

                                                                                                                                              cATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTg       >  2:3795995/1‑149 (MQ=255)

                                                                                                                                              cATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTg       >  2:2381153/1‑149 (MQ=255)

                                                                                                                                                ttCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGga     >  1:4104745/1‑149 (MQ=255)

                                                                                                                                                ttCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGga     <  1:2751790/149‑1 (MQ=255)

                                                                                                                                                ttCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGga     >  2:3454851/1‑149 (MQ=255)

GGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAA  >  CP009273/1872525‑1872820

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