BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,140,030	T→A	7.3%	N23Y (AAT→TAT)	yghW ←	DUF2623 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,140,030	0	T	A	7.3%	835.9 / 41.6	271	N23Y (AAT→TAT)	yghW	DUF2623 family protein
Reads supporting (aligned to +/- strand): ref base T (116/135); new base A (5/15); total (121/150)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 274 total aligned reads displayed.

CATTACCATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGC  >  CP009273/3139881‑3140179

 aTTACCATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTa                                                                                                                                                      <  2:1847228/150‑1 (MQ=255)

  ttACCATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTaa                                                                                                                                                     >  1:3644354/1‑150 (MQ=255)

    aCCATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAAcc                                                                                                                                                   >  2:2311557/1‑150 (MQ=255)

    aCCATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATAGGCACAGAATTTTGTCACATTAAcc                                                                                                                                                   <  1:3475582/150‑1 (MQ=255)

     ccATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCg                                                                                                                                                  >  2:5035225/1‑150 (MQ=255)

      cATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGc                                                                                                                                                 >  1:3148645/1‑150 (MQ=255)

         ctctTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACt                                                                                                                                              <  1:384861/150‑1 (MQ=255)

          tctTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTg                                                                                                                                             >  1:1859734/1‑150 (MQ=255)

          tctTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTg                                                                                                                                             >  1:2090648/1‑150 (MQ=255)

          tctTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTg                                                                                                                                             <  2:1092530/150‑1 (MQ=255)

              tATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCgg                                                                                                                                         >  2:1535770/1‑150 (MQ=255)

               aTCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGc                                                                                                                                        >  2:2918735/1‑150 (MQ=255)

                 ccAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCAt                                                                                                                                      >  2:3870636/1‑150 (MQ=255)

                 ccAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCAt                                                                                                                                      <  1:3034192/150‑1 (MQ=255)

                 ccAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCAt                                                                                                                                      <  2:4159253/150‑1 (MQ=255)

                 ccAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCAt                                                                                                                                      <  1:3806649/150‑1 (MQ=255)

                  cAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATg                                                                                                                                     <  2:803670/150‑1 (MQ=255)

                       cATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTg                                                                                                                                >  2:520405/1‑150 (MQ=255)

                       cATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTg                                                                                                                                >  2:512977/1‑150 (MQ=255)

                          aGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCtt                                                                                                                             <  1:3378628/150‑1 (MQ=255)

                          aGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCtt                                                                                                                             <  1:1484007/150‑1 (MQ=255)

                           gcggcgCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCttt                                                                                                                            >  2:2639778/1‑150 (MQ=255)

                            cggcgCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCtttt                                                                                                                           >  2:4719356/1‑150 (MQ=255)

                            cggcgCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCtttt                                                                                                                           <  2:4969025/150‑1 (MQ=255)

                               cgcgTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCATGCGTTGCTTTTAAt                                                                                                                        <  1:559130/150‑1 (MQ=255)

                               cgcgTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCATGCGTTGCTTTTAAt                                                                                                                        <  1:3146194/150‑1 (MQ=255)

                                gcgTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATc                                                                                                                       <  2:1320881/150‑1 (MQ=255)

                                  gTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATccc                                                                                                                     >  1:2148567/1‑150 (MQ=255)

                                  gTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATccc                                                                                                                     >  1:4293784/1‑150 (MQ=255)

                                    cAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGc                                                                                                                   <  2:3435034/150‑1 (MQ=255)

                                     aGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGcc                                                                                                                  >  2:3825237/1‑150 (MQ=255)

                                     aGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGcc                                                                                                                  <  2:4385125/150‑1 (MQ=255)

                                     aGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGcc                                                                                                                  <  1:4046754/150‑1 (MQ=255)

                                       aaTACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCAt                                                                                                                <  1:734813/150‑1 (MQ=255)

                                        aTACCCGCTTCCCAGGGGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCAtt                                                                                                               <  1:1247544/150‑1 (MQ=255)

                                        aTACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCAtt                                                                                                               <  1:4275686/150‑1 (MQ=255)

                                        aTACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCAtt                                                                                                               >  1:4344169/1‑150 (MQ=255)

                                        aTACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCAtt                                                                                                               >  1:4470988/1‑150 (MQ=255)

                                        aTACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCAtt                                                                                                               <  1:2716586/150‑1 (MQ=255)

                                         tACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTa                                                                                                              >  1:4998969/1‑150 (MQ=255)

                                          aCCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTaa                                                                                                             <  1:519558/150‑1 (MQ=255)

                                           cccGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAg                                                                                                            >  2:383009/1‑150 (MQ=255)

                                            ccGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGc                                                                                                           >  1:3750729/1‑150 (MQ=255)

                                              gCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCt                                                                                                         <  2:2006674/150‑1 (MQ=255)

                                              gCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCt                                                                                                         >  2:967944/1‑150 (MQ=255)

                                               cTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCtt                                                                                                        <  1:1658489/150‑1 (MQ=255)

                                               cTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCtt                                                                                                        >  1:1381964/1‑150 (MQ=255)

                                                ttCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCttt                                                                                                       >  2:3101361/1‑150 (MQ=255)

                                                 tCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTa                                                                                                      >  1:4599656/1‑150 (MQ=255)

                                                  cccAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTAc                                                                                                     <  2:4613410/150‑1 (MQ=255)

                                                     aGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCaa                                                                                                  >  2:2275647/1‑150 (MQ=255)

                                                     aGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCaa                                                                                                  >  1:2999071/1‑150 (MQ=255)

                                                      ggCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCaaa                                                                                                 <  1:4719356/150‑1 (MQ=255)

                                                      ggCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCaaa                                                                                                 <  2:1834612/150‑1 (MQ=255)

                                                        cgcTTATCTGGCTATCTCCTGTCTTTTCGTACATCCGGTTTGAGTAGCCTAATACAAAACCGCGTTTTTAATCGGCACAGAATTTTGTCACATTAACCGATCTGTCGTCATGAGTTGCTTTTAATCCCGCCATTTAGCCTTTACCAAAAt                                                                                               <  1:1449937/150‑1 (MQ=255)

                                                         gcTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATg                                                                                              <  2:2967376/150‑1 (MQ=255)

                                                             aGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTa                                                                                          <  2:2749153/150‑1 (MQ=255)

                                                              gCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTAt                                                                                         <  2:3794662/150‑1 (MQ=255)

                                                                  gCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCAt                                                                                     <  2:5069398/150‑1 (MQ=255)

                                                                   cGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATc                                                                                    >  2:4699572/1‑150 (MQ=255)

                                                                   cGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATc                                                                                    <  1:1140272/150‑1 (MQ=255)

                                                                   cGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATc                                                                                    <  2:1919340/150‑1 (MQ=255)

                                                                      tctcCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGtt                                                                                 <  2:3363810/150‑1 (MQ=255)

                                                                      tctcCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATAAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGtt                                                                                 <  2:3520142/150‑1 (MQ=255)

                                                                       ctcCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTc                                                                                <  1:4117542/150‑1 (MQ=255)

                                                                          cGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCAcc                                                                             >  2:3994680/1‑150 (MQ=255)

                                                                           ggTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCt                                                                            <  1:3823314/150‑1 (MQ=255)

                                                                            gTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTc                                                                           <  2:254308/150‑1 (MQ=255)

                                                                             tCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCtaa                                                                          <  1:2041815/150‑3 (MQ=255)

                                                                              cTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCAt                                                                         <  2:2096544/150‑1 (MQ=255)

                                                                               ttttCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCAtt                                                                        <  1:848099/150‑1 (MQ=255)

                                                                                    gTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCta                                                                   >  2:1198104/1‑150 (MQ=255)

                                                                                     tACAACCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATATTGTCACATTAACCGAACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCtat                                                                  >  1:2599088/1‑150 (MQ=255)

                                                                                       cATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTTATTCGCtatat                                                                <  1:48467/150‑1 (MQ=255)

                                                                                       cATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCtatat                                                                >  1:1934322/1‑150 (MQ=255)

                                                                                              tgtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCa                                                         <  1:663824/150‑1 (MQ=255)

                                                                                                tgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCata                                                       <  1:2554808/150‑1 (MQ=255)

                                                                                                 gAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCatat                                                      >  1:3925482/1‑150 (MQ=255)

                                                                                                          aaTACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGata                                             <  2:2343075/150‑1 (MQ=255)

                                                                                                             aCAAAACCGCGTTGATCATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAttt                                          <  2:1382091/150‑1 (MQ=255)

                                                                                                               aaaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTa                                        >  2:3248240/1‑150 (MQ=255)

                                                                                                                aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAg                                       <  1:3257525/150‑1 (MQ=255)

                                                                                                                 aaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                      <  1:1025880/150‑1 (MQ=255)

                                                                                                                 aaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                      >  2:1624192/1‑150 (MQ=255)

                                                                                                                  aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGACGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                     >  1:904074/1‑150 (MQ=255)

                                                                                                                   ccGCGTTTATAATCGGCACAGAAATTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACAAAAATGGTTATTCATCGTTCACCTCATTAGCTATATATTGTAATATATAGCGATATTTTAGGCa                                    >  1:4088605/1‑150 (MQ=255)

                                                                                                                            taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGAc                           >  2:3204069/1‑150 (MQ=255)

                                                                                                                              aTCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACtt                         >  2:580438/1‑150 (MQ=255)

                                                                                                                                cGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGACCGACtttt                       >  1:4886730/1‑150 (MQ=255)

                                                                                                                                cGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACtttt                       >  1:1278790/1‑150 (MQ=255)

                                                                                                                                 ggCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                      <  1:440894/150‑1 (MQ=255)

                                                                                                                                  gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgc                     <  2:4318363/150‑1 (MQ=255)

                                                                                                                                      aGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCAc                 >  2:573852/1‑150 (MQ=255)

                                                                                                                                         aTTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCt              <  1:933700/150‑1 (MQ=255)

                                                                                                                                            ttGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAtt           <  2:3023095/150‑1 (MQ=255)

                                                                                                                                             tGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAttt          >  2:2094299/1‑150 (MQ=255)

                                                                                                                                              gTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAtttt         <  2:4740953/150‑1 (MQ=255)

                                                                                                                                               tCAGATTTATCGCACTTTCGGCTTGCGTTGCTTTTAATCCCTTCTTTAATCCTTTACCAAAATGGTTATTATTCGTTCACCTCATTCGCTATATATTTTCATTTATAGCGATTTTTTAGGAATTGAACGACTTTTGCGCACGCTAttttt        <  1:5042680/150‑1 (MQ=255)

                                                                                                                                               tCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAttttt        <  2:72487/150‑1 (MQ=255)

                                                                                                                                                     tAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGc  >  1:224482/1‑150 (MQ=255)

CATTACCATCTCTTTATCCAGTCCATAGCGGCGCGTCAGAATACCCGCTTCCCAGGCGCTAAGCTGGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGC  >  CP009273/3139881‑3140179

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: `‑`